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Abstract 

A comprehensive study of exclusive hadronic i?-meson decays into final states con- 
taining two pseudoscalar mesons (PP) or a pseudoscalar and a vector meson (PV) 
is presented. The decay amplitudes are calculated at leading power in Aqcd /nih 
and at next-to-leading order in Og using the QCD factorization approach. The 
calculation of the relevant hard-scattering kernels is completed. Important classes 
of power corrections, including "chirally-enhanced" terms and weak annihilation 
contributions, are estimated and included in the phenomenological analysis. Pre- 
dictions are presented for the branching ratios of the complete set of the 96 decays 
of B~ , B^, and Bg mesons into PP and PV final states, and for most of the 
corresponding CP asymmetries. Several decays and observables of particular phe- 
nomenological interest are discussed in detail, including the magnitudes of the 
penguin amplitudes in PP and PV final states, an analysis of the np system, and 
the time-dependent CP asymmetry in the Kcp and Kr]' final states. 



1 Introduction 



As the B factories [1, 2] continue to accumulate large data samples, an increasing number 
of different i?-decay modes becomes accessible to investigation. Many of these modes 
carry interesting information on CP-violating interactions or hadronic flavor-changing 
neutral currents, but except for a number of decay channels considered "theoretically 
clean" a theoretical framework is required to correct for the effects of the strong inter- 
action. 

For a long time, exclusive two-body 5-decay amplitudes have been estimated in the 
"naive" factorization approach (see, e.g., [3, 4, 5, 6, 7] and references therein) or modifi- 
cations thereof. In many cases this approach provides the correct order of magnitude of 
branching fractions, but it cannot predict direct CP asymmetries due to the assumption 
of no strong rescattering. It is therefore no longer adequate for a detailed phenomeno- 
logical analysis of 5-factory data. Naive factorization has now been superseded by QCD 
factorization [8, 9]. Although not yet proved rigorously, this scheme provides the means 
to compute two-body decay amplitudes from first principles. Its accuracy is limited 
only by power corrections to the heavy-quark limit and the uncertainties of theoretical 
inputs such as quark masses, form factors, and light-cone distribution amplitudes. For 
the charmless decays considered in this paper the lowest-order approximation in QCD 
factorization coincides with "naive" factorization. 

Among the charmless hadronic B decays, the modes S — > tttt and B — > ttK have been 
studied first and most extensively within QCD factorization [8, 10, 11, 12, 13], because 
they are the simplest decays for which a significant interference of tree and penguin 
amplitudes is expected. Other specific final states that have been investigated include 
those with vector mesons and exotic mesons [14, 15, 16, 17], as well as r] or rj' along with a 
pseudoscalar or vector kaon [18]. In the present paper, we complete the phenomenology 
of hadronic B decays into two light pseudoscalar mesons or one pseudoscalar and one 
vector meson, all from the ground-state nonet. We consider the complete set of 96 
decay modes, including decays of Bg mesons, and decays into mesons with fiavor-singlet 
components. We also summarize in a compact notation all the required decay coefficients 
and hard-scattering kernels up to the next-to-leading order in Ug. In a series of recent 
papers [19, 20, 21] Du et al. have considered a subset of these decay modes, including 
decays of Bg mesons. In [22, 23] fits to the data on pseudoscalar-vector meson final 
states based on the generalization of the hard-scattering kernels of [10] to these final 
states have been performed. While we agree qualitatively with the conclusions reached 
from these fits, we beheve that it is currently more useful to investigate in more detail the 
dynamical origin of agreements and discrepancies with experimental data. In this respect 
the analysis presented in this paper is more complete than the fits performed previously 
as far as decay channels, error estimates, and analysis of observables are concerned. 

The motivation for our analysis is threefold: First, the larger set of decay channels 
provides additional information on the CKM phase 7 and on hadronic flavor-changing 
currents, complementary to the tttt and ttK final states. Second, the different strong- 
interaction dynamics underlying specific decay modes provides valuable tests of the QCD 
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factorization framework. Being able to compute correctly a large number of decays, 
including "uninteresting" ones, increases our confidence in the reliability of the approach 
where it is necessary to extract fundamental parameters. An example of this concerns 
the role of so-called scalar penguin operators, which are an important contribution to 
the penguin amplitude in decays to pions and kaons, but which are expected to be 
suppressed for decays into vector mesons [15, 24]. Finally, our analysis provides up-to- 
date expectations for the branching fractions of decay modes yet to be discovered. 

The outline of the paper is as follows: In the rather technical Section 2 we briefly 
review the QCD factorization approach and then collect the notation, definitions, and 
analytic formulae required to compute the complete set of ^ PP and B PV decay 
amplitudes at next-to-leading order. This includes a new notation for the parameteriza- 
tion of flavor amplitudes, which we find more convenient than the conventional notation 
in terms of parameters Oj. Sections 3 and 4 contain a summary of the input parameters 
entering our analysis and an outline of the adopted analysis strategy for the large set of 
final states. Also included in Section 4 is a discussion of the tt^tt" tree amplitude. 

We then analyze separately penguin- dominated AS = 1 and AD = 1 decays, and 
tree-dominated AD = 1 decays in Sections 5 and 6, respectively. We begin with an 
investigation of the magnitude of the penguin amplitude in AS = 1 decays, which dom- 
inates the overall decay rate. The factorization approach makes distinct predictions for 
this amplitude depending on whether the final state is PP, PV, or VP. The comparison 
with existing data provides important information on the usefulness of next-to-leading 
order calculations in QCD factorization and simultaneously on the magnitude of pen- 
guin weak annihilation. We proceed to discuss ratios of branching fractions, and a few 
observations that we find difficult to accommodate in the Standard Model: we quantify 
the correction to the measurement of sin 2/3 in the final states (f)Ks and rj'Ks, and we 
define a new ratio involving the tt^K^ final state that may suggest an anomaly in the 
electroweak penguin sector. We also consider penguin-dominated AD — 1 decays, which 
have small branching fractions, and discuss empirical constraints on the size of weak 
annihilation contributions. 

Section 6 on tree-dominated decays focuses on the information that results from 
the measurement of irp final states. In addition to CP-averaged branching fractions, 
direct CP asymmetries, and certain ratios of branching fractions, we consider the five 
asymmetries that can be defined in the time-dependent study of the final states vr^p^. We 
find that the asymmetry S is well-suited to constrain 7. S together with the analogous 
quantity in B ^ tt'^ti^ decay imply 7 = 70° with an error of about 10° at la. In this 
analysis as in the rest of the paper we adopt the quasi two-body assumption, which 
considers the vector mesons as stable, neglecting interference effects related to the non- 
neghgible widths of the decaying vector mesons. 

In Section 7 we summarize the results for final states containing 77 or 77', which 
complete those of [18], where an additional K or K* in the final state was assumed. We 
briefiy discuss the decays of Bg mesons in Section 8, and conclude in Section 9. In three 
appendices we collect: (i) the complete set of decay amplitudes expressed in terms of 
the amplitude parameters ccj (formerly Oj) and including weak annihilation; (ii) the 
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convolutions of the new hard-scattering kernels with the distribution amplitudes of vector 
mesons up to the second Gegenbauer moment; and (iii) the results of the BaBar, Belle, 
and CLEO experiments, from which we computed the experimental averages used in the 
paper. 

Guide to reading the paper 

Since, in the course of time, this paper has grown to an undue volume, we have structured 
it such that most sections can be read independently of each other. Most of the material 
in Section 2 and Appendices A and B provides an up-to-date summary of the technical 
work necessary to implement QCD factorization at next-to-leading order. The reader 
familiar with the basic ideas of the theoretical framework may read only Section 2.1. The 
reader interested only in specific phenomenological applications should consult Sections 3 
and 4 (omitting perhaps Section 4.3) and the introduction to Section 5, where we specify 
the theoretical input, the general strategy of the analysis, and the definition of specific 
analysis scenarios. The analysis sections 4.3 (tree amplitude in S — > tttt decays), 5 
(penguin-dominated decays), 6 (tree-dominated decays), 7 (final states with rj or rj'), 
and 8 {Bs decays) could then be read independently. 

2 Two-body decay amplitudes 

In this section we detail the theoretical framework within which our results for the decay 
amplitudes are computed. We begin with a brief summary of the QCD factorization 
method. We then introduce a new notation for the basic transition operator T, which 
allows us to describe pseudoscalar and vector mesons in the final two-body state in terms 
of a single expression. The decay amplitude for B — > M1M2 is proportional to the matrix 
element {MiM2\T\B). The transition operator is decomposed into a complete basis of 
"fiavor operators" accounting for the different topologies of the various decay mechanisms 
(tree, penguin, annihilation, etc.). After a summary of light-cone distribution amplitudes 
we collect the results for the coefficients of the different terms in the transition operator. 
This discussion is necessarily rather technical, as one of our goals is to provide a reference 
for the decay amplitudes and hard-scattering kernels in QCD factorization at next-to- 
leading order. Some readers may wish to read only Section 2.1 and then continue with 
Section 3. 

2.1 The QCD factorization approach 

A detailed discussion of the QCD factorization approach can be found in [8, 9, 10]. 
Here we recapitulate the basic formulae to set up the notation. The effective weak 
Hamiltonian for charmless hadronic B decays consists of a sum of local operators Qi 
multiplied by short-distance coefficients Ci and products of elements of the quark mixing 
matrix, Xp = VpbVpj^, where D = d,s can be a down or strange quark depending 
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on the decay mode under consideration, and p — u,c, t. Using the unitarity relation 
Xi^^ + + aJ^^ = we write 

^eff = ;i E ^p""^ ( ^1 ^1 + ^2 Q2 +Y1 + Q77 + Csg Qsyj + h.C. , (l) 

p=u,c ^ j=3 ' 

where Qi 2 are the left-handed current-current operators arising from 1^-boson exchange, 
Qs ... 6 and Q7,...,io are QCD and electroweak penguin operators, and Qy^ and Qsg are 
the electromagnetic and chromomagnetic dipole operators as given in [10]. The effective 
Hamiltonian describes the quark transitions b uuD, b ccD, b Dqq with q = 
u, d, s, c, b, and b Dg, b — > Dj, as appropriate for decay modes with interference 
of "tree" and "penguin" contributions. The Wilson coefficients are evaluated at next- 
to-leading order, consistent with the calculation of operator matrix elements described 
below. For the coefficients of electroweak penguin operators the evaluation described in 
[10], which incorporates some terms normally counted as next-to-next-to-leading order, 
is employed. 

The QCD factorization formalism allows us to compute systematically the matrix 
elements of the effective weak Hamiltonian in the heavy-quark limit for certain two- 
body final states M^Mg. In condensed notation, the matrix element of every operator in 
the effective Hamiltonian is evaluated as 

{M[M!,\Q,\B) = Yl ^r""" ^ij * f^-^^- 

{Mi,M2}e{M{,M^} 

+ * fs^B * fM[^M[ * fu^^M^ , (2) 

where p^^^'^ is an appropriate form factor, $m are leading-twist light-cone distribution 
amplitudes, and the star products imply an integration over the light-cone momentum 
fractions of the constituent quarks inside the mesons. A graphical representation of this 
result is shown in Figure 1. Whenever the spectator antiquark line goes from the B 
meson to one of the final-state mesons we call this meson Mi and the other one M2. 
The sum in the first term on the right-hand side of (2) accounts for the possibility that 
for particular final states the spectator antiquark can end up in either one of the two 
mesons. If the spectator antiquark is annihilated we use the convention that Mi is the 
meson that carries away the antiquark from the weak decay vertex. 

A justification of the factorization formula for final states in which Mi is a heavy 
meson (applicable to D mesons in the heavy-quark limit) can be foTind in [9, 25]. In this 
case the term in the second line of (2) is absent. For final states with two light mesons 
factorization has been proved at order as [8], but a complete proof has not yet been 
given (for recent developments in this direction, see [26]). The factorization formula 
reduces the complicated hadronic matrix elements of four-quark operators to simpler 
non-perturbative quantities and calculable hard-scattering kernels and T/^. In this 
paper we complete the calculation of all relevant kernels at next-to-leading order in ag- 
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Figure 1: Graphical representation of the factorization formula (2). Only one 
of the two form-factor terms is shown for simplicity. 

The light-cone expansion implies that only leading-twist distribution amplitudes are 
needed in the heavy-quark limit. There exist however a number of subleading quark- 
antiquark distribution amplitudes of twist 3, which have large normalization factors for 
pseudoscalar mesons, e.g. for the pion 



For realistic 6-quark masses these "chirally-enhanced" terms are not much suppressed 
numerically. We therefore include in our analysis all quark-ant iquark twist-3 ampli- 
tudes. (The quark-antiquark-gluon amplitude at twist-3 does not have an anomalously 
large normalization.) In order to perform the same analysis for all final states we also 
include the quark-antiquark twist-3 amplitudes for vector mesons, even though there is 
no particular enhancement in this case, being replaced by 2mv/mi, times a ratio of 
two decay constants (see below), with my the vector-meson mass. 

The inclusion of chirally-enhanced terms is important to account for the large branch- 
ing fractions of penguin-dominated decay modes with pseudoscalar final-state mesons, 
such as S — > ttK [10], but it also causes a number of conceptual problems. Factorization 
is not expected to hold at subleading order in Aqcd/"?.?, and, somewhat unfortunately, 
is indeed violated by some of the chirally-enhanced terms [8] . In contrast to the leading- 
twist distribution amplitudes, the twist-3 two-particle amplitudes do not vanish at the 
endpoints but rather approach constants. The kernels Tfj in the first term of the fac- 
torization formula also approach constants at the endpoints (modulo logarithms), and 
hence there is no difficulty with this term. These kernels include the important scalar 
penguin amplitude mentioned in the introduction, conventionally denoted by ag. How- 
ever, the second term in the factorization formula, which accounts for the interactions 
with the spectator quark, contains integrals that are dominated by the endpoint regions 
if the distribution amplitudes do not vanish at the endpoint. These integrals formally 
diverge logarithmically in a pcrturbative framework. This implies a non-factorizable soft 
interaction with the spectator quark, while Mi is formed in a highly asymmetric con- 
figuration, in which one quark carries almost all the momentum of the meson. Similar 
factorization-breaking effects occur in weak annihilation contributions, which are also 
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suppressed by at least one power of Aqcd/^6 in the heavy-quark expansion. 

In [10] we have parameterized the power corrections from chirally-enhanced and weak 
annihilation terms in an ad hoc way and included a variation of the corresponding pa- 
rameters in our error estimate. Based on the analysis of tctt and ttK final states we 
concluded that these factorization-breaking terms introduce a sometimes substantial un- 
certainty into the theory but do not render the framework unpredictive. In particular, 
the data so far give no indication that these corrections should be outside the range 
specified by our error estimate. Given this situation, we follow our previous analysis 
strategy as regards twist-3 and annihilation effects. The inclusion of twist-3 effects for 
vector mesons then allows us to estimate the impact of some power-suppressed effects, 
which turn out to be small. But we should stress that this is far from a complete account 
of such contributions. We also note that some other classes of power corrections have 
been analyzed using the method of light-cone QCD sum rules [27, 28]. 

The pseudoscalar mesons t] and rj', which contain a flavor- singlet component in their 
wave function, require special attention, because they can be formed from two gluons at 
leading twist. The new effects associated with this possibility have been worked out in 
the context of QCD factorization in [18]. For these mesons the factorization formula (2) 
holds provided it is extended by an additional non-local form-factor term. At the order 
we are working the corresponding complications do not appear for vector mesons, since 
they do not have a two-gluon component. There is, however, one novel effect for neutral 
vector mesons, namely that they can be produced via a virtual photon. Although this 
effect is very small, we shall discuss it, because it is of the same order as other electroweak 
contributions to the decay amplitudes. 

2.2 Flavor operators 

In the literature on QCD factorization the result of computing the hard-scattering kernels 
for the various operators in the effective weak Hamiltonian is usually presented in terms 
of "factorized operators" with coefficients ai{MiM2) ■ (The matrix element of a factorized 
operator is simply proportional to a form factor times a decay constant.) The reasons 
for this are largely historical. Besides containing some redundancy, such a notation leads 
to inconveniences for vector mesons since the factorized operator '^q{qb)s-p ® {Dq)s+p 
vanishes, while there exist non-vanishing terms at order as that one would naturally 
associate with the corresponding coefficient a^. Here we introduce a new notation that 
eliminates these inconveniences, and which at the same time is sufficiently general and 
explicit to faciUtate the construction of the decay amplitude for any particular decay 
with either pseudoscalar or vector mesons in the final state. The main point of the 
new notation is to keep track of only the flavor structure of an operator. Its Dirac 
structure, which distinguishes some of the aj(MiM2) coefficients in the conventional 
notation, becomes irrelevant, since operators with different Dirac structure (but identical 
flavor structure) always contribute in the same combination to a particular set of final 
states. 

We first introduce our phase convention for the flavor wave functions. We take the 
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light quark triplet as {u, d, s) and the antitriplet as {u, d, s). Consequently, 

B~ r^ub, Bd^db, B,-^sb. (4) 

Assuming isospin symmetry, the non-singlet members of the light nonet of pseudoscalar 
mesons are given by 

7r° ~ —= iuu — dd) , 7r~ ~ , tt"*" ~ , 

V2 (5) 

K^^ds, K^r^sd, K-r^us, K+r^su, 

and analogous expressions hold for the corresponding vector mesons. The treatment of 
mesons containing a flavor-singlet component in their wave function has been explained 
in detail in the dedicated paper [18], to which we refer the reader for all details. The 
meson rj (and similarly rj' and the vector mesons u and 0) can be treated as a coherent 
superposition of the flavor components 

rjg ~ —= {uu -\- dd) , rjs ss . (6) 

This should not be confused with the representation of the meson state as a superposition 
of flavor states, but simply means that these mesons have matrix elements of ss and uu 
(or dd, these two being equal due to the assumed isospin symmetry) operators that are 
a priori unrelated. Consequently, each meson is described by two decay constants; for 
instance, for the rj meson 
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Similarly, there are two transition form factors, two leading-twist quark-antiquark distri- 
bution amplitudes etc., corresponding to the two distinct flavor components. In addition, 
for the pseudoscalar case there exist contributions proportional to the two-gluon content 
of T] and T]' . Finally, certain contributions to the h — > Dgg amplitude are conveniently 
interpreted in terms of cc flavor components r^c, r]'^ for r] and 77'. The QCD factoriza- 
tion approach presented here is sufficiently general to allow for a discussion of states 
containing flavor-singlet contributions without assumption of a particular flavor-mixing 
scheme. However, in our phenomenological analysis we will adopt the Feldmann-KroU- 
Stech scheme for 77 and rj' mixing, in which it is sufficient to introduce a single mixing 
angle in the flavor basis [29]. The corresponding expressions for the hadronic matrix 
elements (decay constants, form factors, etc.) needed in the calculations of the decay 
amplitudes can be found in [18]. For the vector mesons uj and we will assume ideal 
mixing, so that a; = cUg ~ {uu + dd)/\/2 and = 0s ~ ss. 

We match the effective weak Hamiltonian onto a transition operator such that its 
matrix element is given by 

{M[M'^\H.s\B) = ) {M[M'^\n + n\B) ■ (8) 

p=u,c 
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The two terms account for the flavor topologies of the form-factor and hard-scattering 
terms in (2), respectively, corresponding to the two diagrams of Figure 1. Each term in 
the transition operator contains several operators labeled only by the flavor composition 
of the four-quark final state. Note that below we will include in the definition of Tq also 
some power-suppressed terms, such as weak annihilation contributions. 

Including electroweak penguin topologies there are six different flavor structures one 
can write down for the left diagram of Figure 1. We define (with D — d or s) 

ri = Spu ai(MiM2) A{[qsu][uD]) + Spu a2(MiM2) A{[qsD][uu]) 

where the sums extend over q — u,d, s, and qg denotes the spectator antiquark. The 
operators ^([^Mi?Mi][?M2?M2]) ^-Iso contain an implicit sum over = u,d,s to cover 
all possible S-meson initial states. The coefficients a^{MiM2) contain all dynamical 
information, while the arguments of A encode the flavor composition of the flnal state 
and hence determine the final state to which a given term can contribute. We define 

{M[M^\c^,{M,M,)A{[..][..])\Bj = ca^{M[M^)AM^M^ (10) 

whenever the quark flavors of the flrst (second) square bracket match those of M[ (M2). 
The constant c is a product of three factors of 1, ±l/-\/2 etc. from the flavor composition 
of the B meson and Mi 2 as specifled above. The quantity AmiM2 is given by 

^ r m|Fo^^^^ (0)/m2 ; if Ml = M2 = P, 

Am,m, =^^\ -2mvelj^-pBA^-'''{0)fM, ; if Mi = V, = P, (11) 
[ -2mv elj^ ■ pb F^^^^ (0)/m2 ; if Mi = P, = V. 

Here -F+,o and Aq denote pseudoscalar (P) and vector [V) meson form factors in the 
standard convention (see, e.g., [3]). The decay constants /mj are normalized according 
to 

(7r-(g)|(i7^75M|0) = -i/^g^, {p~{q)\dlt,u\0) = -ifpmpe*^. (12) 

For a vector meson, /y = fy always refers to the decay constant of a longitudinally 
polarized meson. We neglect corrections to the decay amplitudes quadratic in the light 
meson masses, so that all form factors are evaluated at g^ = 0. (At this kinematic point, 
the form factors F+ and Fq coincide.) Parameters referring to the B meson depend 
on whether the decaying meson is Bs or B^^a- This will be implicitly understood in 
the following. The above expression can be simplified by replacing 2'my e* ■ ps — ^ tti^, 
since the left-hand side squared and summed over the polarizations of the vector meson 
gives m% (neglecting again quadratic meson-mass corrections). As expected by angular 
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momentum conservation the final-state vector meson is longitudinally polarized. When 
the replacement above is done the polarization sum for vector mesons must be omitted, 
and the decay rate is simply given by 

16^ |(^(^2|Heff|^)r, (13) 

where 5 = 1/2 if M[ and are identical, and S = 1 otherwise. 

We do not discuss in this paper final states containing two vector mesons. The 
extension to this case is straightforward, since the leading amplitude in the heavy-quark 
limit is the one for two longitudinally polarized vector mesons, which in many ways 
behave like pseudoscalar mesons. 

In order to exemplify the notation consider the decay — > 7r°p°, for which Qs = d 
and D — d. The spectator quark can go to either one of the two mesons, so Mi can be 
7r° or p°. Hence, e.g. 

{n'p'\al{MiM,)J2Am[mBa)^-l [a?(7rV°) ^.o,o + (pV) A,o,o] . (14) 
g 

The order of the arguments of is relevant as will be seen from the explicit expressions 
given in Section 2.4 below. The factor c — 1/2 arises from the flavor wave functions of 
the mesons. On the other hand, 

{n'pV,iMiM,)Y,A{[dd][qq])\Bd)^0, (15) 
g 

since q = u,d contribute equally but with opposite sign for the mesons tt^ and p°. 

We now discuss the flavor structure of the hard-scattering term in (2), i.e., the second 
diagram in Figure 1. Since all six quarks participate in the hard scattering, the generic 
flavor operator in 7^ is of the form -B([^MiQ'Mi] [^M2 9m2] [^s^]) with possible sums over 
quarks from penguin transitions or flavor-singlet conversion g — > qq. We deflne the 
matrix element of a 5-operator as 

{MiM2\B{[...][...]][...])\Bg) = cBM,M,, with BM,M,^±i^fBjMjM,, (16) 

whenever the quark flavors of the three brackets match those of Mi, M2, and Bg. The 
constant c is the same as in (10). The upper sign in the definition of BmiM2 applies when 
both mesons are pseudoscalar, and the lower when one of the mesons is a vector meson. 
This matches the sign conventions for the quantities AmiM2 i^i (H)- The most important 
case of spectator scattering is when the spectator-quark line goes from the B meson to 
a final-state meson, which we then call Mi. These are the hard spectator interactions, 
which are the only terms of leading power in the heavy-quark limit. (We refer to all 
other contributions to as annihilation.) This special case implies qs = (Jmi, leaving 
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six different amplitudes tfiat are in one-to-one correspondence witfi the six A-operators 
defined above, because 

B{[qsqMMM2(lM2][Qsh]) = -^^iMl A{%qMMM2<lM2\) ■ (17) 
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We therefore absorb the spectator scattering contributions into the definition of the 
coefficients af(MiM2). 

With parts of absorbed into 7j, the transition operator 7j now contains all scat- 
tering mechanisms except weak annihilation. In particular, it contains all contributions 
of leading power in the heavy-quark limit. The remaining, power-suppressed annihilation 
part of 7^ is parameterized in its most general form as 

r| = 5puhi{M^M2) Y,B{W][qu][Dh]) + 5puh2{MrM^) ^ 5 ( [wg'] [g'D] ) 

q' q' 

Q,g' q,q' 
+ Spubsi{M,M2) J2B{[uu]W][Db]) + Sp^bs2{M,M2) J2B{[uD][q'q'][ub]) 

q' q' 

Q,q' 9,?' 

+ t^S3,EwiMiM,) B{[qD] [q'q'] [qb]) + lfs4,EwiMiM,) J] ^ e, B{[qq] [q'q'] [Db]) , 

9,9' 9,9' 

where the sums extend over q, q' — u, d, s. The sum over q' arises because a quark- 
antiquark pair must be created hj g ^ q'q' after the spectator quark is annihilated. The 
definitions of the first six coefficients coincide with the corresponding definitions in [10] 
for TTTT and nK final states. Note that the operator g> B{[qq'][q'D][qb]) is redundant, 
because it is equivalent to J2q ^i[qsq][qB]) (which includes an implicit sum over qg). We 
allow this redundancy to keep the parameterization of annihilation effects separate from 
the others. The six new coefficients with subscript '5" contribute only to final states 
containing fiavor-singlet mesons or neutral vector mesons and were not needed in our 
previous analysis. It will be convenient to use the notation 

/5f (M1M2) = 1^ 6f (M1M2) (19) 

AmiM2 

whenever AmiM2 does not vanish. However, for some pure annihilation decays such as 
Bg — > 7r+7r~, where Pf is not defined since Fq^^~*'^(0) = 0, we express the decay amphtude 
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in terms of the coefficients 6^. Tfie redundancy mentioned above implies tfiat a\{MiM2) 
and ^^{MiM^) always appear in the combination 

(M1M2) = al{M^M2) + PUM^M^) . (20) 

It is now straightforward to express the Bq — > M[M2 decay amphtude in terms of 
linear combinations of the coefficients «f and /9,f . The results are collected in Appendix A, 
where we also give a convenient master formula suitable for implementation in a computer 
program, which generates the entire set of 96 decay amplitudes by evaluating a single 
expression. It remains to express the flavor coefficients in terms of the hard-scattering 
kernels of the QCD factorization approach. 



2.3 Distribution amplitudes 

We now summarize the definitions of the light-cone distribution amplitudes for light 
pseudoscalar and vector mesons. The corresponding amplitudes for B mesons have 
been discussed in [30, 31]. While our treatment of the leading-twist distributions is 
completely general, at the level of twist-3 power corrections we work in the approximation 
of neglecting the qqg Fock state of the meson. The motivation for this approximation is 
that it retains all effects with large (chirally-enhanced) normalization factors but neglects 
"ordinary" power corrections of order Aqcb /n^'b- The following discussion relies on the 
analysis of distribution amplitudes in coordinate space presented in [32, 33], while the 
momentum-space projectors are taken from [31]. A more detailed discussion of the 
pseudoscalar case can be found in [10, 34]. 

We recall that in general the coUinear approximation for the parton momenta can be 
taken only after the light-cone projection has been applied. We therefore assign momenta 

U2 p 

k1 = xpi" + + ^ p", H = xp" -k^ + _ ^ _ p" (21) 

2xp-p 2xp-p 

to the quark and antiquark in a light meson with momentum p, where p is a light- 
like vector whose 3-components point into the opposite direction of and x = 1 — x. 
The light-cone projection operator of a light pseudoscalar meson in momentum space, 
including twist-3 two-particle contributions, then reads 



(22) 

where /xp is defined as with defined as in (3). The convention for the projec- 

tion is that one computes tr (M^A) if uAv is the scattering amphtude with an on-shell 
quark and antiquark. The derivative acts on the scattering amplitude A, and it is un- 
derstood that, after the derivative is taken, the momenta ki and k2 are set equal to xp 
and xp, respectively. The overall sign of (22) corresponds to defining the projector (in 
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coordinate space) through the matrix element {P{p)\qf3{z)qa{0)\0) rather than the oppo- 
site ordering of the quark fields. (For simplicity, we suppress the gauge string connecting 
the two fcrmion fields.) The meson t] (and similarly rj', u, and 0) is described by two 
quark-ant iquark distribution amplitudes corresponding to the two independent fiavor 
components rjq and rjs in (6). 

The leading-twist distribution amplitude is conventionally expanded in Gegenbauer 
polynomials, 



$p(a;, ii) = 6x{l — x) 



n=l 



(23) 



since the Gegenbauer moments a^in) are multiplicatively renormalized. When three- 
particle contributions are neglected, the twist-3 two-particle distribution amplitudes are 
determined completely by the equations of motion, which then require 



%{x) 



1 



6 ' 6 

It is not difficult to derive from this the simpler projector [35] 



XX . 



ifp 
4 



k2 ■ ki 



$p(a;) 



(24) 



(25) 



a/3 



The T) and rj' mesons also have a two-particle two-gluon distribution amplitude at leading 
twist, for which we adopt the convention given in [18]. 

The corresponding equations for vector mesons are very similar. In general, the light- 
cone projection operator in momentum space contains two terms, = + M^- 
However, the transverse projector does not contribute to the B PV decay amplitudes 
at leading and first subleading power in Aqcd/"^6. The longitudinal projector is given 

by 



V 



mvfl 



X) 



V 



p-p 



hf{x) 



dv 




(26) 



The polarization vector has been replaced by e* p^/niv, which is correct up to cor- 
rections quadratic in the meson mass. After this replacement no polarization sum must 
be taken after squaring the decay amplitude. The leading-twist distribution amplitude 
<l>y(x) is expanded in Gegenbauer polynomials exactly as in (23), but with Gegenbauer 
moments a^{fi). When three-particle amplitudes are neglected, the twist-3 amplitudes 
that appear in (26) can all be expressed in terms of the leading-twist amplitude ^±{x) 
of a transversely polarized vector meson [33]. We define 

$.(x)= fdv^^- fdv^^^Zy^alMPn+i{2x-l), (27) 

Jo Jx n=0 
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where aQ j_ — 1, and Pn{x) are the Legendre polynomials. The second equation is 
obtained by inserting the Gegenbauer expansion of ^±{x). The equations of motion 
then lead to the Wandzura-Wilczek relations [33] 



= h'^\x)^-(x-x)^^(x), 

(28) 



X 

dv 



—XX ^v{x) 



which allow us to express the twist-3 two-particle projection in terms of the single func- 
tion Comparing the vector meson projection in (26) with the pseudoscalar pro- 
jection in (22) we see that the function ^v{x) is analogous to ^p{x). In analogy with the 
pseudoscalar case we obtain the simpler projector 

= (^M.)-^^*.(.))^^. (29) 

valid in the approximation where one neglects three-particle contributions. Like $p(a;), 
the function $^(x) does not vanish at the endpoints x = 0, 1. We note, however, the 
vanishing of the integral 

'1 

dx^^{x) = 0. (30) 





With these remarks, we can obtain the hard-scattering kernels for vector mesons in 
the final state from those for two pseudoscalars given in [10] by performing the replace- 
ments $p(a;) — >• $v{x), ^p{x) ~^ ^v{x), fp ~^ fv, and by interpreting as in (32), 
(33) below. However, one must pay attention to certain sign changes arising from the 
absence of 75 in the vector-meson projector (29) and in the matrix elements that define 
vector-meson form factors. This leads to the sign alternations in the definitions of a^, 
^4, a^EW' '^4EW ("^1) below, and to sign alternations in the weak annihilation 
terms to be discussed later. 



2.4 Coefficients of the decay operators in QCD factorization 

In the following we give the expressions for the coefficients of the decay operator up to 
the next-to-leading order. First, some comments on the treatment of electromagnetic 
corrections and electroweak penguin effects, and on the treatment of weak annihilation, 
are in order. 



Electroweak penguin effects and electromagnetic corrections 

1) We neglect electromagnetic corrections to the QCD coefficients ai, 0:2, af, and 
a^, since these are much smaller than the next-to-lcading order QCD corrections. 
We also neglect corrections of order a^Cy-io to these coefficients, since the Wilson 
coefficients C7-10 are proportional to the electromagnetic coupling a. 
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2) In leading order the electroweak penguin coefficients and g^v involve only 
the Wilson coefficients Cy-io- We include the QCD corrections of order cCgCy-io 
to «3,EW a-nd a4,EW 

3) We include electromagnetic corrections to ojgg^ and ct-^^^ only when they are 
proportional to the large Wilson coefficients Ci,2,77, ^® include the corrections 
of order aCi^2,7'y, but neglect those proportional to aC^-iQ. 



Weak annihilation 

1) We neglect weak annihilation mechanisms involving photons (7 q'q'), which 
in fact would introduce more operators in (18), i.e., we drop annihilation terms of 
order a Ci, but keep terms of order Cj. The reason for this is that the former can 
never be CKM-enhanced when one of the large Wilson coefficients Ci,2 is involved. 

2) We use an approximation where all singlet annihilation coefficients are set to zero 
except for Ps3, see below and [18]. 

In the following we first give the generic results for the ctj and Pi coefficients applicable 
to all final states, and then discuss contributions specific to particular pseudoscalar or 
vector mesons. 



A-operators (generic results) 

The coefficients of the flavor operators can be expressed in terms of the coefficients 
af defined in [8, 10] as follows:^ 

ai{MiM2) = ai(MiM2), 
a2{MiM2) = a2{MiM2) , 



af(MiM2) 



a^(MiM2) - (MiAfs) ; if M1M2 = PP, VP , 
af (M1M2) + af (MiiVfa) ; if M1M2 = PV , 



' ^ a^(MiM2) - a^(MiM2) ; ifMiM2 = yP, ^ ^ 



<EW(^1^2) 
<EW(A^1^2) 



4(MiM2) - a?(MiM2) ; if M1M2 = PP, VP , 
al{MiM2) + a?(MiM2) ; if M1M2 = PV , 

a?o(MiM2) + rf 2 al{M,M2) ; if M1M2 = PP, PV , 
a?o(MiM2) - 2 a^(MiM2) ; if M1M2 = FP . 



-'^The numerical values of the coefScicnts ai{MiM2) also depend on the nature of the initial-state B 
meson. This dependence is not indicated explicitly by our notation. The same remark applies to the 
annihilation coefficients fof defined below. 
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Figure 2: Leading-order contribution to the coefficients af. The weak decay 
of the b quark through a four-fermion operator is represented by the black 
square. The 6-quark hne comes in from the left. The outgoing line to the right 
represents the quark in Mi. The two lines directed upward represent the meson 
M2. The spectator antiquark is not drawn, because it does not participate in 
the hard scattering. 



Note that the order of the arguments in af (M1M2) and of (M1M2) is relevant. For pions 

X 



and kaons, the ratios r^^ are defined as 



while their generalizations to rj and rj' can be found in [18] . All quark masses are running 
masses defined in the MS scheme, and nig denotes the average of the up and down quark 
masses. For vector mesons we have 



= (33) 

where the scale-dependent transverse decay constant is defined as 

e*)\qa^.q'\0) = f^ip.el - p.e;) . (34) 

Note that all the terms proportional to r^^ are formally suppressed by one power of 
in the heavy-quark limit. Numerically, however, these terms are not always 

small. 

The general form of the coefficients af at next-to-leading order in as is 



+ 



Nr 47r 



Vi{M2) + —Hi{MiM2) 



+ P!{M,) , (35) 



where the upper (lower) signs apply when i is odd (even). It is understood that the 
superscript 'p' is to be omitted for i = 1,2. The quantities Vi{M2) account for one-loop 
vertex corrections, Hi{MiM2) for hard spectator interactions, and Pf (M1M2) for penguin 
contractions. We now present the explicit expressions for these objects, extending the 
results of [10] to the case when Mi or M2 is a vector meson. 
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?V W \C Vi 



Figure 3: Next-to-leading order vertex contribution to the coefficients a". The 
meaning of the external lines is the same as in Figure 2. 



Ni{M2) 



(36) 



Leading order (Figure 2). The leading-order coefficient Ni{M2) is simply the normal- 
ization integral of the relevant light-cone distribution amplitude $p,y or It follows 
that 

' ; i = 6, 8 and M2 = y , 
1 ; all other cases. 

The special case for vector mesons arises because of the vanishing integral (30). This 
implies the absence of scalar penguin contributions for vector mesons at tree level, since 
the coefficients ag g originate from {V — A) ® {V + A) penguin operators in the weak 
Hamiltonian, which must be Fierz-transformed into the Form (—2) (5* — P)®{S + P) to 
match the quark flavors in the second [S + P) current with those of the meson M2. 



Vertex terms (Figure 3). The vertex corrections are given by 



Vi{M2) = { 



dx 



dx 



121n— -18 + ^(x) 

-121n — + 6-^(1 -x) 
A* 



dx^m^i.^) — Q + h{x) 



i = 1-4, 9, 10, 
i = 5,7, 
i = 6,8, 



(37) 



with 



g{x) — 3 ( In X — in 



+ 



1-x 



2Li2(a:) -ln'^x + 



2 In a; 
1 — X 



(3 -|- 2i7r) Inx — {x 1 — x) 



(38) 



h{x) — 2 Li2{x) — hi^x — (1 + 27ri) In x — (x <-> 1 — x) . 

The expression for g{x) has already been presented in [10], while the kernel h{x) is new. 
The constants —18, 6, —6 are scheme dependent and correspond to using the NDR 
scheme for 75. The light-cone distribution amplitude is one of the leading-twist 
amplitudes $p,v, depending on whether M2 is a pseudoscalar or vector meson, whereas 
^rn2 is one of the twist-3 amplitudes We recall that ^p{x) = 1, so / dx ^rn2{^) [—6-1- 
h{x)] = —6 for pseudoscalar mesons, which rcprodTiccs the result of [10]. On the other 
hand, because of (30) the scheme-dependent constant —6 does not contribute for vector 
mesons as it should be, since there is no leading-order contribution that could cause a 
scheme dependence at next-to-leading order. 



16 



Figure 4: Next-to-leading order penguin contribution to the coefficients aF^. 



Penguin terms (Figure 4)- At order ctg a correction from penguin contractions is present 
only for i = 4, 6. For i = 4 we obtain 



_ln_ + --G^^(.,) 



+ C3 



6 Id 6 



(39) 



where n/ = 5 is the number of light quark flavors, and s„ = 0, Sc = {rric/mbY are mass 
ratios involved in the evaluation of the penguin diagrams. Small electroweak corrections 
from Cy-io are neglected within the approximations discussed above. The function 
Gu^is) is given by 



Gm2{s)= / dxG{s — ie,l — x)^M2{^) ^ 
Jo 

G{s,x) — —A / duu{l — u)l'a[s — u{l — u)x] 
Jo 

2(12s + 5x-3x In s) 4V4s - x (2s + x) 



(40) 



9x 



arctan 



X 



As — X 



The interpretation of is the same as in the discussion of vertex corrections. For 
i — 6, the result for the penguin contribution is 



4: TTlb 2 ^ , . 



+ C3 



^ln^ + ^-GM.(0)-GM.(l) 



^ In ^ - (n/ - 2) GmM - GmA^c) - GmM 

O jJi 



-2C: 



eflf 



(41) 
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Figure 5: Hard spectator-scattering contribution to the coefficients a^. The 
meaning of the external hues is the same as in Figure 2, but the spectator-quark 
hne is now included in the drawing. 



if M2 is a pseudoscalar meson, and 



Ci Gm2{sp) + 



GmM + Gm,{1) 



+ {C^ + a) {rif - 2) GmM + GmAsc) + GmA'^) 



(42) 



if M2 is a vector meson. In analogy with (40), the function GM^i^) is defined as 

Gm2 {s) = / dxG{s-ie,l- x) {x) . (43) 

As mentioned above we take into account electromagnetic corrections only for 0:3 
and Q!4EW' ^^"^ they are proportional to the large Wilson coefficients Ci,2 and 

C^^. These corrections are present for i = 8, 10 and correspond to the penguin diagrams 
of Figure 4 with the gluon replaced by a photon. (An additional contribution for neutral 
vector mesons will be discussed separately below). For i = 10 we obtain 



a 



For i = 8 we find 

P|(M2) 



{C, + N,C2) 



4: rrib 2 ' 
_ln_ + --G^^(.^) 







a 



9nNr 



(Ci + iVeCs) 



if M2 is a pseudoscalar meson, and 

Pi{M2) = - 



a 



3ln- + --GM.(.p) 



(Ci + A^cC2)Gm.(Sp) 



l-x 



(44) 
(45) 

(46) 



if M2 is a vector meson. 
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Hard spectator terms (Figure 5). The correction from hard gluon exchange between M2 
and the spectator quark is given by 



H.{M,M,) fdx fdy 



^M2{x)^mM ^ M^^M2{x)^raM 

+ = 



xy 



xy 



(47) 



for i = 1-4,9,10, 



BmiM2 '^b 



iMiM2 



1 .1 

dx I dy 
Jo 



xy 



xy 



(48) 



for i = 5,7, and Hi{MiM2) = for i = 6, 8. In these results Xb is defined by [8] 



(49) 



with ^b{0 of ^'^o hght-cone distribution amphtudes of the B meson. We recall 
that the term involving is suppressed by a factor of Aqcd / "t.^ in heavy-quark power 
counting. Since the twist-3 distribution amplitude $mi(l/) docs not vanish at y 1, the 
power-suppressed term is divergent. We extract this divergence by defining a parameter 



X^^' through 



Jo y Jo y Jo y 

' dy 



$mi(y)-<^'mi(i; 



Jo [y\+ 



(50) 



The remaining integral is finite (it vanishes for pseudoscalar mesons since $p(y) = 1), but 
X^^ is an unknown parameter representing a soft-gluon interaction with the spectator 
quark. Since the divergence that appears in an attempt to compute this soft interaction 
pcrturbativcly is regulated by a physical scale of order Aqcd (i-e. at y ~ Aqqb /n">-b), 
we expect that X^ ~ In(m6/AQCD), however with a potentially complex coefficient, 
since multiple soft scattering can introduce a strong-interaction phase. A consequence 
of this is that power corrections to the spectator interaction, including chirally-enhanced 
ones for pseudoscalar mesons, are unavoidably model dependent. As in [10], our model 
consists of varying X^^'^ within a certain range (specified later) and to treat the resulting 
variation of the coefficients as an uncertainty. We also assume that X^^^ is universal, 
i.e., that it does not depend on Mi and on the index i of Hi{MiM2). 

In [10] the convolution integrals relevant to S — > PP decays were evaluated exphcitly 
up to second non-trivial order in the Gcgcnbaucr expansion. In Appendix B we present 
the corresponding results for convolution integrals involving the function $^(a;), which 
appear only in 5 — > PV decays. 
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Figure 6: Additional contributions (representative diagrams shown) to the 
decay amphtude when M2 = 77 or 77'. Left: Charm-loop contribution to the b 
Dgg amplitude. Right: Soft spectator-scattering diagram, which contributes 
to Aa^{Mir]q^s) in (52). The gluon to the left is soft and part of the operator 
whose matrix element defines the non-local form factor T^^^'^. The gluon 
that originates at the weak vertex is semi-hard. 



A-operators (peirticulEir results) 

In addition to the generic expressions given above, there exist additional contributions 
to the decay amplitudes for final states containing an 77 or 77' meson (related to their 
two-gluon content), and for final states containing one of the neutral vector mesons 
a;, (j) (related to their coupling to the photon). 

The contributions specific to t;*^'^ mesons have been discussed in [18]. There are two 
effects relevant to the A-operators. First, the leading contribution of the h — > Dgg 
amplitude can be interpreted in terms of a "charm decay constant" of the 77^'^ meson, 
which is calculable in a Aqcd/"^c expansion (left side of Figure 6). To account for this 
effect the term 

[5pc«2(Mi7^W) + a?(Mi7^W)] A{\qsD][cc]) (51) 

is added to in (9). Since this contribution is very small, we will only need the 
leading-order expressions for the new ctj coefficients, which coincide with the general 
leading-order result. We will also neglect the corresponding electroweak penguin effect, 
which is even smaller. Second, there exists an additional form-factor type contribution 
Aaf (Mi77^'s) to the flavor-singlet coefficient a^{Mi'qq}s) ■, which is given by 

A«3(iWir?,, J - g^^^ ^8. {t^h) Udx + . . . I ^^^^^^ , (52) 



where /^^ = v^m^A^ with = 0.5 GcV serves as a typical scale for the semi-hard 
subprocess (right side of Figure 6), and the dots stand for a contribution from the 
leading- twist two-gluon distribution amplitude of 77 or 7/', which vanishes for asymptotic 
distribution amplitudes. The new feature of (52) is that this contribution is proportional 
to a non-local "form factor" J^^~'^^(0) (rather than the usual local form factors F — Fq, 
Ao) defined in terms of the B Mi matrix element of a bilocal quark-antiquark-gluon 
operator [18]. The ratio of the two form factors in (52) is of order unity in the heavy- 
quark limit and has been estimated to be close to 1; however, there is a large uncertainty 
associated with this estimate. 
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1 

Figure 7: Additional contributions (representative diagrams shown) to the 
decay amphtude when M2 — p^,uj, 4>- Wavy hues denote a photon. 

For the vector mesons a;, the particular contributions are related to the electro- 
magnetic penguin diagrams shown in Figure 7. These contributions are very small, but 
it is interesting to discuss them from a conceptual point of view. In both diagrams the 
photon propagator is canceled. The left diagram where the photon originates from the 
electric dipole operator Q^^ is therefore effectively a local vertex. The other diagram in 
the figure can be calculated explicitly in an expansion in Aqcd/'^tiq when the quark in 
the loop is heavy (with mass mg). In the case of a charm quark this produces logarithms 
of the form ln(m6/mc), which could be summed by introducing a cc decay constant and 
distribution amplitude for the vector mesons, similar to the treatment of the b Dgg 
amplitude for 77 and 1]'. The quark loop contains an ultraviolet divergence, which must 
be subtracted in accordance with the scheme used to define the Wilson coefficients. The 
scale and scheme dependence after subtraction is required to cancel the scale and scheme 
dependence of the electroweak penguin coefficients Cy-io- When the quark in the loop is 
massless (light quarks), the loop integral is infrared and ultraviolet divergent. The ultra- 
violet divergence is subtracted as before. The infrared divergence must be factorized into 
the longitudinal decay constant of the vector meson. As a consequence, parameters such 
as fpo become scale-dependent.^ The reason for this is that the operator q^y^q has a non- 
vanishing anomalous dimension in the presence of electromagnetic interactions. With 
these remarks it is straightforward to compute the additional contribution Aa^ -^^(MiV) 
to the electromagnetic penguin amplitude a^-^,^{MiV) when V ~ p^,u},(j). Within our 
approximation of keeping only the aCi^2,77 corrections we find 

Aal^^{M,V) = ^ \4{C\ + N,C2) (I - 6,, In ^ - In ^"j - 9Cf\ . (53) 

Here refers to the renormalization scale of the Wilson coefficients Ci{ij), and v to the 
scale of the decay constant fv{i^)- The v dependence of (53) cancels the v dependence 
of the leading-order decay amplitude to order a. In our analysis we set u — /i. 

S-operators/weak annihilation (generic results) 

We now turn to a discussion of power-suppressed weak annihilation effects. Except for 
certain sign changes for final states with vector mesons the "non-singlet" annihilation 

^Strictly speaking, we should then distinguish fpO [v) from fp+ , which is scale-independent. However, 
given the smallness of the effect we set the two decay constants equal at f = /U. 
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Figure 8: Weak annihilation contributions. 



coefficients can be taken from [10]. We consider 6-quark decay and use the convention 
that Ml contains an antiquark from the weak vertex with longitudinal momentum frac- 
tion y. For non-singlet annihilation M2 then contains a quark from the weak vertex with 
momentum fraction x. The basic building blocks when both mesons are pseudoscalar 
are given by (omitting the argument M1M2 for brevity) 







4 


= 0, 






4 


= 0, 


^3 




4 





/■ 

^0 



1 



y{l - xy) x^y^ 



+ 



Jo 



= nos I dxdy <{ r^^ ^M^i^) ^miiv) - 



x(l — xy) xy'^ 

2y 



xy 



(54) 



xy{l- xy) 
2{l + x) 



2x 



Ai = Tras Cdxdy |rf ^ $m.(x) ^mM ^ti±^ + rf ^mM ^m,{x) 

JO K X y 



xy{l - xy) 

2(1 + ?/) 
xy'^ 



When Ml is a vector meson and M2 a pseudoscalar, one has to change the sign of the 
second (twist-4) term in A\, the first (twist-2) term in A\, and the second term in A^ 
and A3. When M2 is a vector meson and Mi a pseudoscalar, one only has to change the 
overall sign of A\. 

In (54) the superscripts 'i' and '/' refer to gluon emission from the initial and final- 
state quarks, respectively (see Figure 8). The subscript 'A;' on A^'-'' refers to one of the 
three possible Dirac structures Fi ® which arise when the four-quark operators in 
the effective weak Hamiltonian are Fierz-transformed into the form (^ife)ri(^2Q'3)r2) such 
that the quarks in the first bracket refer to the constituents of the B meson. Specifically, 
we have k ^ 1 ior {V - A) {V - A), k ^ 2 for (V - A) ® {V + A), and A; = 3 for 
{—2){S — P) <S> {S + P). The power suppression of weak annihilation terms compared 
to the leading spectator interaction via gluon exchange is evident from the fact that 
annihilation terms arc proportional to J'b rather than fs^B/^B- 

In terms of these building blocks the non-singlet annihilation coefficients are given 
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by 



If 



If 



7V2 



C,A\ + CiAl, + Ai) + N.CeAi 



CiA\ + CqA'^ 



C,A\ + Cj{Al + Ai) + N,CsAi 



C^oA{ + CsAi, 



(55) 



omitting again the argument M1M2. These coefficients correspond to CTirrcnt-currcnt 
annihilation (&i,&2)) penguin annihilation (63,64), and electroweak penguin annihilation 
(6f^,6f^), where within each pair the two coefficients correspond to different flavor 
structures as defined in (18). 

The weak annihilation kernels exhibit endpoint divergences, which we treat in the 
same manner as the power corrections to the hard spectator scattering. The divergent 
subtractions are interpreted as 



^0 



dy 

y 



fdy 
Jo 



Iny 

y 



A J 



(56) 



and similarly for M2 with y ^ x. The treatment of weak annihilation is model-dependent 
in the QCD factorization approach, and the explicit results of this subsection are useful 
mainly to keep track of overall factors from Wilson coefficients and color. We treat 
as an unknown complex number of order ln(mfe/AQCD) and make the simplifying 
assumption that this number is independent of the identity of the meson M and the 
weak decay vertex. (The first assumption will be relaxed in a specific scenario, where 
we allow different Xa for the three cases PP, PV , and VP.) Since the treatment of 
annihilation is model-dependent anyway, we further simplify our results by evaluating 
the convolution integrals with asymptotic distribution amplitudes $(x) = ^11(3:^) = Qxx, 
^p{x) = 1, and ^v{x) 3(x — x). We then find the simple expressions 



^ 67ro;, (rf 



rf )(X1-2X^ + 



TT 

"3 



(57) 



4 



67ra, (rf + rf ^) {2X\ - Xa) , 



.M2 
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Figure 9: Examples of flavor-singlet annihilation graphs. Diagrams where 

both gluons attach to the constituents of the B meson belong to the B rj^'^ 
form factors. The fourth diagram and similar ones with one of the two gluons 
attached to the constituents of the B meson do not contribute to /J^g. 




and A{ = A2 = Q when both final state mesons are pseudoscalar, whereas 



TT 



3(X^-4+V)+rfvf (Xl-2X^) 



Al 67ras 



(58) 



Ai ^ 67ra, [3rf ^ {2Xa - 1)(2 - Xa) - rf (2X1 - ^a)] , 

and A{ — AI — when Mi is a vector meson and M2 a pseudoscalar. For the opposite 
case of a pseudoscalar Mi and a vector M2, one exchanges in the previous 

equations and changes the sign of A^. 



S-operators/weak annihilation (particular results) 

The calculation of the singlet weak annihilation coefficients bsi is even more uncertain. 
Some of the diagrams that can contribute to these quantities are shown in Figure 9. 
Since we consider only diagrams formally proportional to ag, the second meson M2 
must have a two-gluon content, hence the effect vanishes for vector mesons. We shall 
also make the approximation adopted in [18], in which the small electroweak singlet 
annihilation coefficients fe^a.Ew and &54,ew are neglected, and in which also 651 and 652 
are not computed and only 653 is kept. The reason for neglecting 651 and 652 is not 
that these terms are smaller than 653, but that in decay amplitudes they always appear 
together with the large tree coefficients cti and a2- On the other hand, 653 appears in 
conjunction with the singlet penguin amplitude . The singlet weak annihilation effect 
is then confined to final states with an r] or rj' meson. 

With these approximations only the first three diagrams shown in the figure (and the 
corresponding crossed diagrams) have to be calculated. The result is [18] 

(.S3(M,P„.) = ^ (a + N^C,) 5fi I'dy !£i|lM ^ . (59) 

where ^pg{x) is the leading- twist two-gluon distribution amplitude of P, whose Gegen- 
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bauer expansion reads 



^Pg{x) = 5B^^{n) x^x'^ {x-x) + ... . (60) 

The coefficient as well as all higher Gegenbauer moments vanish for — > oo. 

Expression (59) is endpoint- divergent as any other weak annihilation term. Introducing 
X^'^ as before, and truncating the Gegenbauer expansion of ^pg{x) after the leading 
term, we obtain the estimate 

hsz{M^P,,s) = %i.C^ + ^cCe) h^:asB^'{^xn)X^' , (61) 

c 

which we will use in our numerical analysis below. The process 77* r]^''' can in principle 
be used to constrain i?^^(/i). In our normalization convention the analysis performed in 
[36] gives -B^^(l GeV) = 2 ± 3. Note, however, that the second Gegenbauer moment of 
the singlet quark-antiquark amplitude is much smaller, of order —0.1. 



3 Input parameters 

The predictions obtained using the QCD factorization approach depend on many input 
parameters. First there are Standard Model parameters such as the elements of the CKM 
matrix, quark masses, and the strong coupling constant. Of those our results are most 
sensitive to the strange- quark mass, which sets the scale for the ratios for pseudoscalar 
mesons defined in (32). (We work with a fixed ratio rriq/ms, hence implicitly depends 
on rUs-) Some branching ratios, and in particular the direct CP asymmetries, are also 
very sensitive to the value of \Vub\ and the weak phase 7 = aig{V*ij). Next there are 
hadronic parameters that can, in principle, be determined from experiment, such as 
meson decay constants and transition form factors. In practice, information about these 
quantities often comes from theoretical calculations such as hght-cone QCD sum rules or 
lattice calculations. The corresponding uncertainties in the form factors are substantial 
and often have a large impact on our results. Finally, we need predictions for a variety of 
light-cone distribution amplitudes, which we parameterize by the first two Gegenbauer 
coefficients in their moment expansion. Experimental information can at best provide 
indirect constraints on these parameters. Fortunately, it turns out that the sensitivity of 
our predictions to the Gegenbauer coefficients is usually small. The notable exception is 
the color- suppressed tree amplitude «2(^i^2), which shows a considerable dependence 
on the first inverse moment of the i?-meson distribution amplitude (A^) and the second 
Gegenbauer moment of the light mesons through the hard-spectator interaction. 

A summary of the input parameters entering our numerical analysis is given in Ta- 
ble 1. Some additional parameters related to 77 and rj' mesons, such as their decay 
constants, form factors, and the rj^' mixing angle in the quark-flavor basis, can be 
found in [18]. Not given in this reference are the values of the second Gegenbauer mo- 
ments of the quark-antiquark twist-2 distribution amplitudes for the flavor components 

Tjg^ and r)s \ for which we take a^"'^ = ± 0.3. 
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Table 1: Summary of theoretical input parameters. All scale-dependent quan- 
tities refer to /i — 2 GeV unless indicated otherwise. Parameters related to rj 
and rj' are given in [18]. 



QCD scale and running quark masses [GeV] 




mb{mb) 


rricimb) 


m,(2GeV) 




0.225 


4.2 


1.3 ±0.2 


0.090 ±0.020 


0.0413 


CKM parameters and 5-meson lifetimes [ps] 


\Vcb\ 




7 


r(fi-) 






0.041 ±0.002 


0.09 ±0.02 


(70 ± 20)° 


1.67 


1.54 


1.46 


Pseudoscalar-meson decay constants [MeV] 




fK 


fB 


fBs 


131 


160 


200 ± 30 


230 ± 30 


Vector-meson decay constants [MeV] 


fp 








209 ± 1 


218 ±4 


187 ±3 


221 ±3 


Transverse vector-meson decay constants [MeV] 


J p 






J 


150 ± 25 


175 ± 25 


150 ±25 


175 ± 25 


Form factors for pseudoscalar mesons 


:at g2 = 0) 


^0 


i^B^K 
^0 


T^B.^K 
^0 


0.28 ±0.05 


0.34 ± 0.05 


0.31 ± 0.05 


Form factors for vector mesons (at = 0) 


.B^p 
^0 


aB^K* 
^0 


AB~*u) 
^0 


aB.^K" 
^0 


aBs-^4> 

^0 


0.37 ±0.06 


0.45 ±0.07 


0.33 ±0.05 


0.29 ±0.05 


0.34 ± 0.05 


Parameters of pseudoscalar-meson distribution amplitudes 






a§ 


\b [MeV] 


Ab. [MeV] 


0.1 ±0.3 


0.2 ±0.2 


0.1 ±0.3 


350 ± 150 


350 ± 150 


Parameters of vector-meson distribution amplitudes 


p p 


af, afl 


K* K* 




^21 ^2,± 


0.1 ±0.3 


0.2 ±0.2 


0.1 ±0.3 


0±0.3 


0±0.3 
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A few additional comments are in order. The values for the decay constants of 
pseudoscalar mesons and longitudinally polarized vector mesons can be determined with 

good accuracy from experimental data on the leptonic decays 7r~ — > /xz/^, K~ — > /xz/^, the 
semileptonic decay r ^ p~ Ur, and the electromagnetic decays V e^e~ with V = p°, 
LU, or (p. We have updated the values obtained in [4] by using the most recent results 
for the various decay rates. We will neglect the small uncertainties on these parameters 
in our numerical analysis. The values we take for the decay constants of the B and 
Bg mesons are in the ball park of many theoretical calculations using QCD sum rules 
and lattice gauge theory. The values for the heavy-to-light form factors are close to the 
results of light-cone QCD sum rules where available [37, 38, 39]. In other cases we base 
our values on a crude estimate of SU(3) flavor symmetry breaking effects. The B — >■ r/^'^ 
form factors receive an unknown two-gluon contribution. We therefore parameterize the 
form factor as [18] 

-^0 - -^0 -T~ + ^2 7=- • (D^j 

/t V3/^ 

(For Bg decay replace F^'^'^ f^^,) by Fq^~*^ /*(,).) Information on F2 can in principle 

be obtained from semileptonic B ^ rjlu decay at = 0. At present, however, the 
parameter F2 is completely undetermined, and for lack of better knowledge we adopt 
the value F2 = 0, for which Fq^^ = 0.23 and Fq~'^ = 0.19. The modes with rj' in the 
final state are rather sensitive to this choice. This introduces an additional theoretical 
uncertainty not taken into account in the error ranges given below (see [18] for the 
dependence of the B K^*^ri^'^ modes on the choice of F2). The values for the transverse 
decay constants and Gegenbauer moments of vector mesons are rounded numbers taken 
from [38], however we have inflated the small errors quoted there. ^ 

The quark masses are running masses in the MS scheme. Note that the value of 
the charm-quark mass is given at // = m;,. The ratio Sc — {mc/mbY needed for the 
calculation of the penguin contributions is scale independent. The values of the light 
quark masses are such that = r^. Finally, the value of the QCD scale parameter 
corresponds to ^^(M^) = 0.118 for the two-loop running coupling in the MS scheme. 
The corresponding results for the Wilson coefficients Cj are tabulated in [10]. 

As discussed in detail in [10], there are large theoretical uncertainties related to 
the modeling of power corrections corresponding to weak annihilation effects and the 
chirally-enhanced power corrections to hard spectator scattering. As in our earlier work 
we parameterize these effects in terms of the divergent integrals Xh (hard spectator 
scattering) and Xa (weak annihilation) introduced in (50) and (56). We model these 
quantities by using the parameterization 

XA^il + QAc"'^)!!!^; QA<1, A^ = 0.5GeV, (63) 



^To facilitate the comparison with the results of [10, 18], where they overlap, we note the following 
changes of input parameters relative to those papers: nis was (110 ±25) McV in [10] and (100 ± 25) MeY 
in [18]; IKb/V^bl was 0.085±0.017, fs was (180±40)MeV, and F^^^^ was 0.9/if//^F(f ^'^ in [10]; af 
was 0.3 ± 0.3, t{B~) was 1.65 ps, and T{Bd) was 1.56 ps in [10, 18]. 
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and similarly for Xh- Here (pA is an arbitrary strong-interaction phase, which may be 
caused by soft rescattering. In other words, we assign a 100% uncertainty to the "default 
value" Xa = In^mMf-^h) ~ 2.4. Unless otherwise stated wc assume that Xa and Xh arc 
universal for all decay processes. Finally, in the evaluation of the hard-scattering and 
annihilation terms we evaluate the running coupling constant and the Wilson coefficients 
at an intermediate scale Hh ~ (AqcD^ft)^^^ rather than fj, ~ mf,. Specifically, we use 



4 Analysis strategy for B and decays 

Wc are now in a position to discuss the phcnomenological implications of our results, and 
to compare them to the experimental data compiled in Appendix C. Unless otherwise 
stated, all results for branching fractions and decay rates refer to an average over CP- 
conjugate modes, i.e., a result for Bt{B — > /) actually refers to 

^ [Br{B f) + Br{B /)] . (64) 

Our convention for the direct CP asymmetry follows the standard minus B" conven- 
tion 

Br(E° ^ /) - Br(E° ^ /) 
^ = Br(SO ^ /) + Br(SO ^ /) ' ^^^^ 

which is opposite to the sign convention employed in our previous paper [10]. 



4.1 Outline 

In the following sections we will compare the data with our theoretical results. Besides 
presenting our default values for the branching fractions and CP asymmetries, along with 
a detailed estimation of the various sources of theoretical uncertainty, we will consider a 
series of scenarios of specific parameter sets which elucidate correlations between different 
quantities and their sensitivity to hadronic parameters. Sections 5 to 7 deal with decays 
of B~ and B^ mesons (and their CP conjugates), while Bg decays are treated in Section 8. 

The decay modes are categorized according to which flavor topology (tree, penguin, 
or annihilation) gives the dominant contribution to the amplitude. The AS = 1 decay 
modes are always penguin dominated and arc discussed first, in Section 5. They typically 
have "large" branching ratios of order few times 10~^ to few times 10~^. The theoretical 
predictions for these modes often suffer from large uncertainties due to the strange-quark 
mass and due to power corrections contributing to the weak annihilation coefficient 
which is part of the dominant amplitude — al + (3^. These two sources of uncertainty 
are almost fully correlated between the various decay channels, since they always con- 
tribute to the penguin coefficients a^. We illustrate methods that allow one to extract 
from data the magnitudes and strong phases of the dominant penguin contributions in 
the three topologies for which M1M2 — PP, PV, and VP. Once the measurements 
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become more precise, it will be possible to reduce significantly the corresponding corre- 
lated theoretical uncertainties. For the time being we consider different scenarios for the 
unknown annihilation contributions as representatives for modifications of the penguin 
amplitude a1 away from our default parameter choices. There are also some AD = 1 
decays which are penguin or annihilation dominated. In the last part of Section 5 we 
show how they can be used to derive constraints on annihilation parameters. 

Most of the AD — 1 decay modes are dominated by tree amphtudes, which do not 
suffer from a large sensitivity to light quark masses or chirally-enhanced power correc- 
tions. These decays are studied in Section 6. In many cases the dominant theoretical 
uncertainty arises from the variation of CKM parameters or form factors and decay con- 
stants. The latter source of uncertainty can be reduced once better data on semileptonic 
or leptonic B decays become available. The sensitivity to CKM parameters is not a 
theoretical limitation but rather provides access to \Vub\ and 7, which is important for 
CKM fits. 

Wc omit from our discussion in Sections 5 and 6 the decays with an rj or rj' meson 
in the final state. As discussed in [18] these modes are characterized by a complicated 
interplay of many different fiavor topologies and suffer from additional large theoretical 
uncertainties due to 77-77' mixing, the two-gluon component in the r]^'^ wave functions, 
and an annihilation contribution to the B i]'^'^ semileptonic form factors. We will, 
however, consider modes with an cj or meson in the final states, taking the fiavor wave 
functions tu = Uq ^ {uu + dd) / \/2 and (j) = (j)s ^ ss corresponding to ideal mixing, and 
neglecting singlet annihilation contributions (which for vector mesons vanish within our 
approximations). The large number of final states with r) or rj' mesons is then analyzed 
briefiy in Section 7. The decays with additional pseudoscalar or vector kaons have already 
been discussed in [18]. Here we give the results for the complete set, concentrating on 
B^ — > Tc^T]^'^ and B~ — > p~r]^'^ decay, which have sizable branching fractions. 

4.2 Simplified expressions for the decay amplitudes 

Appendix A contains the exact expressions for all 96 decay amplitudes in terms of the 
fiavor parameters ai and Pi. While these results (along with the expressions for the 
flavor operators collected in Section 2.4) are used in our numerical evaluations, for a 
phenomenological analysis is will be very useful to have simpler, approximate expressions 
at hand, which capture the dominant contributions to the amplitudes. These can be 
obtained by making the following approximations: 

• AS — 1 Decays: We neglect annihflation contributions proportional to Xu \ since 
they are strongly CKM suppressed with respect to the corresponding terms pro- 
portional to Ac . This amounts to setting (32, (3^, and [32 to zero in Ap- 
pendix A. Of the electroweak penguin contributions we only keep ag since all 
other electroweak penguin terms are strongly suppressed. We neglect in particular 
all electroweak annihilation contributions. 
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• AD = 1 Decays: We neglect all electroweak contributions, since they are never 
CKM enhanced and formally of order a, hence most likely smaller than unknown 
QCD corrections of order a^. 

Within these approximations the various contributions to the decay amplitudes can be 
classified as tree topologies (cki and q;2), penguin topologies (af and a^), an electroweak 
penguin topology (03 ew)' annihilation topologies The penguin annihilation 

contribution is always combined with the penguin contribution into the combi- 
nation 64 = 0:4 + The electroweak penguin contribution is kept only for = 1 
decays. 



4.3 The B n tree amplitude 

Before discussing the various penguin and tree-dominated i?-meson decay modes it is 
instructive to consider the process B~ — > 7r~7r°, which among the charmless modes 
discussed here is the single pure tree decay (within the approximations mentioned above) . 
For the corresponding CP-averaged branching ratio we obtain 



10^Br(S- ^ tt-tt") = {G.ltl-}) x 



IK 



ub\ 



-^0 



-(0) 



0.0037 0.28 



(66) 



where the largest sources of uncertainty come from the parameter of the 5-meson 
distribution amplitude, the second Gegenbauer moment of the pion distribution ampli- 
tude, and the quantity Xh- This value can be compared with the experimental result 
Br(S~ — > 7r~7r°) = (5.3 ± 0.8) • 10~^. The good agreement of the central values indi- 
cates that the magnitude of the tree amplitude is obtained naturally with the default 
set of parameters. However, given the substantial uncertainty in the overall normal- 
ization due to \Vu})\ and Fq^^''(0), a similarly good agreement could also result if the 
amplitude coefficients Q;i,2(7r7r) were to take values rather different from our expecta- 
tions |Q;i(7r7r)| = 0.991^;^^ and |Q;2(7r7r)| = 0.201^;}]'. 

The comparison could be made independent of the values of \Vuh\ and F(f^'^(0) if 
the semileptonic B ^ ttIu rate were measured near = 0. One could then perform a 
direct measurement of the tree coefficients via the ratio [40, 41] 



r(B- 



TT TT 



dr{BO^n+l-u)/dq\,^^ 



= Sn \Vud\ |«i(vr7r) + a2{mT) 



(67) 



In QCD factorization we find |ai(7r7r) -|- Q;2(7r7r)| = 1.17_q;q7, which yields the value 
(0.66 to,ll) GeV^ for the above ratio. 

The CLEO collaboration [42] has measured the semileptonic decay spectrum in three 

bins. Using their result for the lowest bin. 



I 



^^•^^^ dBr 

2 /■ dO 



dq^ 



(S" ^ TT+ri/) = (4.31 ± 1.06) • 10" 



(68) 
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and making the assumption (fulfilled in all form- factor models) F^^^{q'^) > F^^''(0) 
for < 8GeV^, we obtain the upper bound 

\Vub\ ^"(0) < V6.05 • 10-5 ( T{Bd) ^ dq^ \p^\^ j = 1.22 • 10"=^ (69) 

at 90% confidence level, which is already close to our default input value (1.03 ± 0.30) • 
10"^. (The central experimental value would give the upper bound 1.03 • 10"^.) With 
Br{B- 7r-7r°) > 4.0 • 10"^ at 90% CL we find 

|ai(7r7r) + a2{7T7r)\ > 0.87 (90% CL) , (70) 

or > 1.10 for central values. Although the bound at 90% CL is not yet in an interesting 
range, the limit obtained using central values shows the potential of the ratio (67) for 
constraining Q;i,2(7r7r) directly. This underlines the importance of the semileptonic decay 
spectrum for understanding the pattern of non-leptonic tttt final states. 

Significantly stronger limits can already be obtained if one assumes a model for the 
dependence of the form factor, which may be guided by QCD sum-rule calculations 
or lattice data at large q^. A fit of such form factor models to the CLEO decay spectrum 
and lattice data has recently been performed [43] , resulting in 

\Vub\ F+^''(0) = (0.83 ± 0.16) • 10-3 (71) 

and a small form factor F:^'~*'^(0) = 0.23 ± 0.04. This would require a sizable value of 
|Q;i(7r7r) -|- Q;2(7r7r)| to account for the B~ — > n~n^ branching fraction. 

A scenario where a2{'K7r) is large may be of interest, since the current experimental 
value of the ratio 

—-— ^ = 0.86 ±0.15 72 

Bt{B- TT-TfO) ^ ^ 

is not in good agreement with the QCD factorization result 1.47lo;43 (with 7 = 70° and 
|^&/K;6| — 0.09 fixed). This is often taken as an indication for destructive interference 
of tree and penguin contributions to B^ n~7r~^, implying a large value of 7 or a large 
strong phase, in contradiction with factorization. However, another possibility could be 
that S° 7r~7r° is enhanced relative to B^ — > 7r~7r"'" because the color-suppressed tree 
amplitude a2{nn) is large, of order 0.4 — 0.5. To keep the B^ — > 7r~7r° branching ratio 
in agreement with its experimental value we must then assume that the B ^ n form 
factor F^^^{0) is about 0.25, on the lower edge of the range of model predictions, but 
in agreement with the trend indicated by the semileptonic decay spectrum. Below, when 
we define a set of parameter scenarios different from the default parameters to exhibit 
correlations among the decay modes, we shall therefore consider this "large-Q;2" scenario 
as one of our options. 

This scenario is not as contrived as it may seem, since it can be realized by simply 
taking the the second Gegenbauer moment of the pion to be a2 = 0.4 and Xb = 200 MeV, 
both at the boundary of our parameter region. Indeed, there is no rigorous information 
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available on the parameter of the 5-meson distribution amplitude, and although the 
light-cone distribution amplitude of the pion is now believed to be close to the asymptotic 
form already at small scales (disfavoring a large Gegenbauer moment 0^2), the calculations 
of the low-energy hadronic processes from which this information is extracted may have 
significant theoretical uncertainties. Furthermore, the heavy-to-light form factors at 
— are usually taken from QCD sum-rule calculations, but it now appears that these 
form factors exhibit a substantially more complex dynamics than B ^ D form factors 
or pion transition form factors. It remains to be investigated whether the dynamics of 
these transitions at large momentum transfer is adequately represented in the sum-rule 
framework. 

5 Penguin-dominated decays 

All B~ and B*^ two-body decays with kaons in the final state are dominated by penguin 
amphtudes. The modes with a single kaon are based on A^" = 1 transitions and have 
"large" branching ratios of order few times 10~^ to few times 10~^. Modes with two 
kaons proceed through AD = 1 transitions and therefore are expected to have smaller 
branching fractions. In this section we focus on decays without r} or rj' mesons in the 
final state. These modes will be studied in Section 7. 

A successful prediction for the branching ratios and CP asymmetries in rare B decays 
requires theoretical control over the magnitudes and relative strong-interaction phases 
of tree topologies, penguin topologies, electroweak penguin topologies, and annihilation 
topologies. The fact that generically a given decay amplitudes receives several interfering 
contributions complicates the interpretation of the experimental data in terms of flavor 
topologies. We expect that our predictions for annihilation effects (as modeled by the 
complex parameters qa and ^pa) and for the strong phases of the various contributions 
are afflicted by the largest theoretical uncertainties. The former effects are incalculable 
within QCD factorization and so can only be estimated using a simple model. Strong 
phases are predicted to vanish in the heavy- quark limit. The leading contributions to 
these phases are then of order ag and calculable, or of order A/m^ and incalculable. 
Since in practice the two expansion parameters and A/m^ are not very different, we 
expect that our perturbative analysis at one-loop order can only give a rough estimate of 
the values of strong phases. While we do trust the generic prediction that these phases 
are small, we expect signiflcant deviations from the values obtained at order ag from 
higher-order perturbative and power corrections. This expectation is indeed supported 
by an analysis of hadronic B decays using the renormalon calculus [44] . 

Given these limitations, it is instructive to test the reliability of our predictions for the 
tree, penguin, and annihilation topologies in decays without interference of the different 
topologies. This will probe the magnitude of these topologies but not their strong- 
interaction phases. Note that the penguin and annihilation terms cannot be completely 
separated phenomenologically, since and (3^ always enter in the combination 6^ — 
•^4 + Ps- However, some other annihilation contributions can be probed directly. 
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A discussion of the tree contribution to the decay B" — > 7r~7r° has aheady been 
given in Section 4.3. In the following we present a detailed study of penguin and anni- 
hilation terms, using penguin-dominated decay processes based on AS' = 1 and AD = 1 
transitions. For each class of decays we give explicit expressions for the decay am- 
plitudes in terms of flavor parameters, adopting the approximations described earlier. 
We then present numerical predictions for the CP-averaged branching fractions and CP 
asymmetries and compare them with the world-average experimental data compiled in 
Appendix C. These results are always obtained using the exact representations of the 
decay amplitudes in terms of flavor parameters as given in Appendix A. 

In the tables we flrst present our "default results" (column labeled "Theory") along 
with detailed error estimates corresponding to the different types of theoretical uncer- 
tainties detailed in Section 3. The first error shown corresponds to the variation of the 
CKM parameters |Vc6|, iKbl, and 7 ("CKM"), the second error refers to the variation of 
the renormalization scale, quark masses, decay constants (except for transverse ones), 
form factors, and (in later sections) the ri-r]' mixing angle ("hadronic 1"). The third er- 
ror corresponds to the uncertainty due to the Gegenbauer moments in the expansion of 
the light-cone distribution amplitudes, and also includes the scale-dependent transverse 
decay constants for vector mesons ("hadronic 2"). Finally, the last error reflects our 
estimate of power corrections parameterized by the quantities Xa and Xh ("power"), 
for which we adopt the form (63) with gA,H < 1 and arbitrary strong phases ^Pa,h- 

In order to illustrate correlations between errors we explore four parameter scenarios 
in which certain parameters are changed within their error ranges while all others take 
their default values. Specifically, these scenarios are defined as follows: 

• Scenario SI ("large 7"): 

To study the dependence of the various observables on the CP-violating phase 7 
we use the default parameter values but set 7 = 110° instead of 70°. 

• Scenario S2 ( "large 02" ) : 

As discussed in Section 4.3, the experimental branching ratios for S — > tttt modes 
can be reproduced in QCD factorization cither by using a large value of 7 > 90°, 
or by enhancing the ratio a2(7r7r)/ai(7r7r) with respect to its default value of about 
0.2. This can be done using the form factors F(f ^'^(0) = 0.25 and F^^^{0) = 0.31, 
the Standard Model parameters iKib/Kbl = 0.08 and = 70MeV, as well as the 
Gegenbauer moments = 200 MeV and a2 — 0.4, all of which are within the 
error ranges specified earlier. 

• Scenario S3 ("universal annihilation"): 

The dominant penguin coefficient a^, which includes the annihilation contribution 
/^f , is rather sensitive to the parameter Xa- While B decays into two pscudoscalar 
mesons are well described by a small annihilation contribution, certain B PV 
decay amplitudes favor a penguin coefficient that is larger than its default value. 
Here we adopt the simplest model in which a moderate enhancement of 0:4 is 
attributed to weak annihilation. Specifically, we treat Xa as a universal parameter 
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obtained by using Qa — ^ and a strong phase (p^ = —45°. The sign of this phase 
is not predicted, but is chosen such that the sign of the direct CP asymmetry 
Acp{7r'^K~) agrees with the data. 

• Scenario S4 ("combined"): 

Here we consider a combination of S2 and S3, but with the more moderate pa- 
rameter values rUs = 80MeV and = 0.3, as well as non-universal annihilation 
phases cpA = -55° (PP), ipA = -20° (PV), and cpA = -70° (VP). The signs 
of these phases are not predicted. As a result, our predictions for the signs of 
CP asymmetries in this scenario must be taken with caution. While each of the 
previous scenarios fails to describe well some classes of decay modes, scenario S4 is 
an attempt to combine certain parameter choices (all within our theoretical error 
ranges) in such as way as to obtain a good description of all currently available 
experimental data. In view of this, this is our currently favored scenario. 

In Section 5.4 we will also study a modification of scenario S3 in which we use the large 
value qa — with a universal phase (pA — —60°. 



5.1 Penguin-dominated AS = 1 decays 

We start by presenting in Tables 2 and 3 our results for the CP-averaged branching 
fractions and direct CP asymmetries for the decays B nK, B nK*, B Kp, as 
well as for the modes B — > Ku and B — > K(j). The B — > ttK modes have the largest 
branching fractions, of order (l-2)-10~^. The data show that the corresponding rates for 
the PV modes B — > nK* and B Kp are smaller by about a factor of two, indicating 
a sizable suppression of the penguin amplitudes in the cases with a final-state vector 
meson. QCD factorization predicts small direct CP asymmetries for most decay modes, 
however with a few exceptions. All predictions for CP asymmetries agree with the data 
within errors. 

Adopting the approximations described in the previous section, the B — > ttK decay 
amplitudes are given by 



V2Ab --*%0K A, 

Ago^TT+K- = A^K 



Spu Oil + Oil 



+ A 



Kit 



Spu Oi2 + Spc fois^EW 



(73) 



V2A 



BO 



A 



ttK 



Ol'a 



Spu OL2 + Spc 2'^3,EW 



where it is understood that each term must be multiplied with Ap and summed over 
p = u,c. The order of the arguments of the coefficients a^{MiM2) and f3f{MiM2) is 
determined by the order of the arguments of the AmiM2 pref actors. The expressions for 
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Tabic 2: CP-averaged branching ratios (in units of 10^^) of penguin-dominated 
B — > PP decays (top) and B — > PV decays (bottom) with AS = 1. The 
theoretical errors shown in the first column correspond (in this order) to the 
uncertainties referred to as "CKM" , "hadronic 1" , "hadronic 2" , and "power" 
in the text. The numbers shown in the remaining columns correspond to 
different scenarios of hadronic parameters explained in the text. 



Mode 



Theory 



SI S2 S3 



S4 



Experiment 



B- 
B- 
B^ 
B^ 



TT+K' 



19.3 



-1.9+11.3+1.9+13.2 
-1.9- 7.8-2.1- 5.6 
1 1 1 +1.8+5.8+0.9+6.9 
-'--'-•-'--1.7-4.0-1.0-3.0 
1 a q+2.6+9.6+1.4+11.4 
-'-"•'-'-2.3-6.5-1.4- 4.8 
7 n +0.7+4.7+0.7+5.4 
' -0.7-3.2-0.7-2.3 



18.8 

14.0 
20.3 
6.5 



20.7 
11.9 
18.8 
8.3 



24.8 
14.0 
21.0 
9.3 



20.3 
11.7 
18.4 
8.0 



20.6 ± 1.3 
12.8 ± 1.1 
18.2 ±0.8 
11.2 ± 1.4 



B- 
B- 



TT+K*- 



3.6 



+0.4+1.5+1.2+7.7 
' -0.3-1.4-1.2-2.3 
q q +1.1 +1.0+0.6 +4.4 
'-'•'-' -1.0-0.9-0.6-1.4 

q q+1.4+1.3+0.8+6.2 

'^•'J -1.2-1.2-0.8-1.6 
n 7+0.1+0.5+0.3+2.6 
"-"• ' -0.1-0.4-0.3-0.5 



3.4 

5.5 
5.9 
0.6 



2.2 
2.6 
2.4 
0.4 



7.3 
5.4 
6.6 
2.1 



8.4 
6.5 
8.1 
2.5 



13.0 ±3.0 

< 31 
15.3 ±3.8 

< 3.6 



B- 
B- 
B' 
B^ 



K-p+ 
K^p'> 



tr o +0.6+7.0+1.5+10.3 
'-'•° -0.6 -3.3 -1.3- 3.2 
9 +0.9+3.1+0.8+4.3 
^•" -0.9-1.4-0.6-1.2 
7 A +1.8+7.1+1.2+10.7 
'•^-1.9-3.6-1.1- 3.5 
A e +0.5+4.0+0.7+6.1 
^•"-0.5-2.1-0.7-2.1 



5.6 
1.3 
4.3 
5.0 



13.6 
6.0 
13.9 
8.4 



10.8 
4.7 
12.5 
7.5 



9.7 
4.3 
10.1 
6.2 



< 48 

< 6.2 
^.9 ±2.2 

< 12 



B- 
B^ 
B- 
5° 



► K-uj 

> K-(t) 



3.5 
2.3 

4.5: 



-1.0+3.3+1.4+4.7 



, +0.3+2.8+1.3+4.3 
-0.3-1.3-0.8-1.3 
+0.5+1.8+1.9+11.8 
-0.4-1.7-2.1- 3.3 
^ +0.4+1.7+1.8+10.6 



-0.4-1.6-1.9- 



3.0 



1.9 
1.9 

4.4 
4.0 



7.9 
6.6 
2.5 
2.3 



5.8 
4.5 
10.1 
9.1 



5.9 
4.9 
11.6 
10.5 



5.3 ±0.8 

5.1 ± 1.1 

9.2 ± 1.0 
7.7 ±1.1 



the B — > T^K* and B — > Kp amphtudes are obtained by replacing (tt, K) — > (tt, K*) and 
(p, K) , respectively. The remaining amplitudes take the form 



V2Ab-^k-u, 
Ab--.k 



Aku> 



a2 + 2a^ + Spc gC^a ew 



+ Au,R 

+ A^R 6?4^ , 



AbO- 



>K°(t> — AkcI> 



a- 



pc 2'^3,EW 



Q;| + Spc 2 Ois^EW 



(74) 



An important feature of all these modes is that, within our approximations, annihilation 
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Tabic 3: Direct CP asymmetries (in units of 10^^) of penguin-dominated B 
PP decays (top) and B — > PV decays (bottom) with = 1. See text for 
explanations. 



Mode 



Theory 



SI 



S2 



S3 



S4 



Experiment 



B- 
B- 
B^ 
5° 



n+K- 



A Q +0.2 +0.3 +0.1 +0.6 
"-"•^ -0.3-0.3-0.1 -0.5 
7 1 +1.7+2.0+0.8+9.0 
-1.8-2.0-0.6-9.7 
A r +1.1+2.2+0.5+8.7 
^■"J-1. 1-2. 5-0.6-9.5 
_o o+l. 0+1. 3+0.5+3. 4 
"^■"^-0.8-1.6-1.0-3.3 



0.9 
5.7 
3.6 
-3.5 



0.8 
6.3 
3.0 
-3.6 



0.4 
-1.3 
-3.6 
-1.2 



0.3 
-3.6 
-4.1 

0.8 



-2 ±9 
1± 12 
-9 ±4 
3 ±37 



B- 
B- 

5° 



n+K*- 



1 c +0.4+0.6+0.5+2.5 
-'-■"-0.5-0.5-0.4-1.0 
o 7+2.1+5.0+2.9+41.7 
"■ ' -2.6-4.3-3.4-44.2 



9 -1 +0.6+8.2+5.1 
^■^ -0.7-7.9-5.8 



62.5 
64.2 

19 o +4.0+4.7+2.7+31.7 
-3.2 -7.0 -4.0 -35.3 



1.7 
5.2 
1.2 

-15.5 



1.6 
8.7 
1.7 
-17.7 



0.8 
-20.4 
-33.8 

1.5 



0.8 
-6.5 
-12.1 

1.0 



26 ±35 



B- 
B- 
B^ 
B^ 



K~p+ 
K^p'> 



0.3 



+0.1+0.3+0.2+1.6 



1 Q c +4.5 +6.9 +3.7 +62.7 
-'-"^■'J -5.7 -4.4 -3.1 -55.4 
_o Q +1.3+4.4+1.9+34.5 
"^■° -1.4 -2.7 -1.6 -32.7 
7 r +1.7+2.3+0.7+8.8 
' ■^ -2.1-2.0-0.4-8.7 



0.3 
-27.3 
-6.6 

6.9 



0.4 
-9.3 
-2.6 

5.5 



0.4 
26.6 
18.3 

1.4 



0.8 
31.7 
20.0 
-2.8 



26± 15 



B- 
B- 



■ K-uj 

■ K-(j) 



_7 Q +2.6+5.9+2.4+39.8 
' ■o -3.0 -3.6 -1.9 -38.0 

_Q 1 +2.5+3.0+1.7+11.8 
o^-*- -2.0-3.3-1.4-12.9 
1 R +0.4 +0.6 +0.5 +3.0 
-■-■"-0.5-0.5-0.3-1.2 
1 7+0.4+0.6+0.5+1.4 
-■-■' -0.5-0.5-0.3-0.8 



-14.5 
-9.6 

1.6 
1.7 



-5.5 
-4.3 
1.7 
1.9 



18.4 
0.0 
0.6 
0.9 



19.3 
3.7 
0.7 
0.8 



0± 12 

3±7 
19 ±68 



effects enter the decay amplitudes only through the effective penguin coefficients 0.% — 
O-t + 13^- As a result, these effects can be readily constrained using data. This explains 
why the B — > t:K decay modes have been employed frequently (by many authors and 
using various analysis strategies) as a way of constraining the phase 7. 

Consider first our default results labeled "Theory" . While the experimental results 
for the B — > t:K modes are well reproduced within errors, it appears that in the case 
of the modes B nK* and B — * K(f) the central values for the branching fractions 
obtained from QCD factorization are consistently lower than the data by a factor 2 to 
3. No significant discrepancy (at the present level of precision) is seen in the modes 
B — > Kp and B — > Kcu, in which the dominant penguin coefficient is al{VP) instead 
of al{PV) (we may call these modes B — > VP). We will come back to this observation 
below. 

In all cases the theoretical predictions suffer from large uncertainties due to the 
strange-quark mass (largest contribution to the second error), and due to power correc- 
tions contributing to the weak annihilation coefficient /Sf, which is part of the dominant 
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penguin amplitude (largest contribution to the last error). We stress, however, that 
these two sources of uncertainty are almost fully correlated between the various decay 
channels, since they always contribute to the penguin coefficients a^. Below we will dis- 
cuss how these coefficients can be directly determined from the data. Once the measure- 
ments become more precise it will be possible to reduce significantly the corresponding 
correlated theoretical uncertainties. 

As a final comment, let us add that for some observables in the decays under con- 
sideration here (especially B K(f) and B — > ttK) the present experimental values are 
inconsistent with Standard Model predictions at the level of 2 to 3 standard deviations. 
We will discuss these discrepancies at length later. While we will not pursue New Physics 
explanations of these discrepancies in the present work, the reader should bear in mind 
that, unless these discrepancies disappear with more precise data, one might expect large 
deviations also in other observables in these decays, such as their CP-averaged branching 
fractions. Therefore, deviations of our predictions from experiment do not necessarily 
imply a failure of the QCD factorization approach. 



5.1.1 Magnitudes and phases of penguin coefficients 

The decay amplitudes for the three sets of decays B — > ttK, B — > ttK*, and B — > pK 
shown in (73) are particularly simple. In the B — > PP modes various ratios of CP- 
averaged decay rates have been used to derive information about the weak phase 7 as 
well as the strong phase of the ratio al/ai [10, 45, 46, 47, 48, 49]. 

The amplitudes for the decays B~ 7r~K^, B~ — * p~ ■, and B~ tt~K*^ are 
determined in terms of the coefficients a\. QCD factorization predicts that ^ a^io a. 
good approximation. Given that |Au''/Ai''''| ~ 0.02, it then follows that the contribution 
with p — u can be neglected at the present level of accuracy of the experimental data. 
(This assertion would be invalidated if a significant direct CP asymmetry was found 
in any of these modes.) It follows that r{B~ tt~K^) oc |Q;4(7rfC)p and Aq^{B~ — > 
Ti^ K^) to a very good approximation, and similar relations hold for the other two 
decays. These decays therefore provide a direct test of the coefficient 64, which contains 
the penguin annihilation parameter /^g. These penguin annihilation contributions, as 
well as potential other non-perturbative contributions to the penguin coefficient a% are 
sometimes referred to as "charming penguins" [50]. 

Uncertainties related to heavy-to-light form factors can be eliminated by normalizing 
the pure-penguin rates to the pure-tree rate for the decay B~ 7r~7r° discussed in 
Section 4.3. Specifically, the magnitude of the penguin coefficient q.1{tiK) in the decays 
B — > "kK can be probed by considering the ratio 



a, 



[ttK) 





Vuh 








fK . 



T{B- 



71- 



2T{B- 



TT TT 



0^ 



1/2 



= 0.103 ±0.008 ±0.023, (75) 



where the numerical value uses the experimental measurements of the decay rates. The 
first error is experimental while the second reflects the uncertainty in \Vub\- This ratio. 



which is inversely proportional to the quantity called £exp in 



provides a crucial 
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Figure 10: Complex ratios al{MiM2)/{ai + a2){'n"n') for the cases M1M2 = ttK 
{PP) and 71 K* (PV) (horizontal axis = real part, vertical axis = imaginary 
part). The rings show the experimental values of the magnitudes of the two 
ratios with (light) and without (dark) including the uncertainty from iKtb/Kftl- 
The three contour lines in each plot correspond to = 1 (solid), 2 (dashed), 
and 3 (dashed-dotted). The triangular-shaped area in the first plot shows the 
constraint on the strong phase derived from the direct CP asymmetry in 
7r~K^ decays. The various symbols indicate theoretical points corresponding 
to different parameter scenarios: the dots on the solid contours correspond 
to S3, the squares correspond to S4, and the circles on the dashed contours 
correspond to the "large annihilation scenario" with qa = 2 and ifA = —60°. 



test of the QCD factorization approach, since it directly probes the magnitude of the 
dominant penguin-to-tree ratio in B ^ PP decays. With default parameters, QCD 
factorization predicts the values |a;i(7r7r) + a2(vr7r)| = 1.17 and \al{7iK) + pi^{nK)\ = 
I (-0.099 - O.OUi) - 0.0091 = 0.108, yielding 0.093 for the ratio of coefficients on the 
left-hand side, which is in excellent agreement with the data. Note that without the 
weak annihilation contribution from we would obtain 0.085 for this ratio, which is 
still in good agreement with the data. For comparison, this ratio equals 0.061 in naive 
factorization. It follows that the dominant penguin amplitude in i? ^ PP decays is 
correctly predicted in QCD factorization without any adjustment of parameters, and 
that the (incalculable) penguin annihilation contribution to this amplitude is a small 
effect, which in fact is consistent with the magnitude expected for a power correction. 
Specifically, we find that |/33/q;4| = 0.091q;o9. At least for the case of B ^ nK modes 
there is thus very little room for large non-perturbative contributions to the dominant 
penguin amplitude besides those already included in QCD factorization. 

This observation can be used to place a constraint on the quantity Xa in (63), which 
we use to parameterize our model estimate for the annihilation terms. In the first plot in 
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Figure 10 we compare our theoretical results for the ratio Q;4(7rAr)/(Q;i + a2)(7r7r) with the 
experimental value of its magnitude obtained from (75). The point with error bars gives 
our central value with errors from parameter variations (including Xa and Xh)- The 
solid "onion-shaped" curve indicates the variation allowed with = ^ and arbitrary 
annihilation phase (pA- This is the region allowed by the error analysis of [10]. The 
larger onion-shaped curves correspond to values of the annihilation parameters outside 
our usual error analysis, namely qa — '2 (dashed curve) and Qa — ^ (dashed-dotted 
curve). They are included here to indicate that for such large values of qa agreement 
with the data requires fine-tuning of the annihilation phase cpA- We will see below that 
such large qa values are all but excluded by current data. In fact, the B ttK decays 
alone {PP case) do not allow for values of qa above about 2.5 once data from direct CP 
asymmetries are included (triangular shape, see below). 

The penguin coefficients in the B PV modes B ttK* and B Kp can 
be probed via relations analogous to (75). For the case of B~ ti~K*^ data exist, 
and the corresponding result is shown in the second plot in Figiirc 10. In that case 
the ratio f-^j fx in (75) must be replaced by f-^/ fK*-, but there is still no dependence 
on hadronic form factors. We observe that the central value for the magnitude of the 
penguin coefficient 04 (tt^*) obtained in QCD factorization is significantly smaller than 
the experimental value, and the two barely overlap when theoretical uncertainties are 
taken into account. Unfortunately, the branching fraction for the decay B^ p'^ 
has not yet been measured. We therefore cannot test whether the penguin coefficient in 
the B — > VP modes B — > pK is correctly predicted in QCD factorization. 

It is also instructive to determine ratios of penguin coefficients that are independent 
of iKib/K&l) using the relations 



fK_ 

Ik* 



T(B- 



TT 



-X*0) 



T{B- TT-^o) 



1/2 



0.58 ± 0.07 (exp.) , 



-(0) 



A, 



(0) 



V(B- 



p-K') 



1/2 



(76) 



(not yet observed) , 



which can be used to relate the dominant penguin coefficients in the PV modes 
B — > ttK* and B — > pK to the corresponding coefficient in the PP mode B ttK. 
While unfortunately no data are available yet that would allow us to deduce the relative 
magnitude of the penguin coefficients in the B — > PV and B — > VP modes, the first ratio 
above provides a first clue about the magnitude of the penguin coefficient in S — > PV 
modes, which is seen to be significantly smaller than the corresponding coefficient in 
B — > PP decays. Qualitatively, this reduction can be understood in terms of the fact 
that the quantity Og in the relation = a^+r^ Og in (31) vanishes at tree level for the case 
where M2 is a vector meson. In fact, using default parameters we predict an even more 
drastic reduction of the penguin coefficient, \al{7[K*) + f3i^{']TK*)\ — |(— 0.030 — 0.002z) — 
0.005| = 0.035, where the power-suppressed twist-3 term in the projector (29) contributes 
about 10% to the result for al. This leads to a ratio \al{7:K*) /al{7iK)\ = 0.32. Note 
that the annihilation contribution is now potentially much bigger, l/^g/a^l ~ 0-16lo.'i4) 
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and also much more uncertain. The main reason for the enhancement of this ratio is 
the reduction of compared with the case of S — > PP decays. We expect a similar 
situation to hold in the case of VP modes, for which we expect \al{pK) + f3^(pK)\ — 
I (0.037 + 0.003i) + 0.0071 = 0.044 and \(3i/a2\ = O.lStoH 

Additional information can be obtained by combining the observables for the decays 
S° — > n'^K^ and B~ 7r~^°, as well as for the corresponding B — > PV modes [46, 48]. 
Neglecting the small parameter (which is justified at the level of a few percent) and 
introducing the tree-to-penguin ratio rpM = 
obtain 

r(fiO ^ 71+K-) , ,2 

-Kfm = 7^775 = 1-2 C0S7 RerpM + Ffm , 

r{B ^TT K^) (77) 

Rfm ■ ^cp(-B° tt^K') = —2 sin7 ImrpM ■ 

In the B TiK sector the experimental values of these ratios are -Rfm = 0.96 ± 0.07 
and RpM ■ ^cp{R^ ~^ 7i^K~) = —(8.2 ± 3.8)%. The second relation in (77) provides 
useful information on the relative strong- interaction phase of the ratio ai{nK) / al{nK) . 
Assuming that ai{TrK) has about the same phase as ai{Tnr) (both are expected to be 
almost real and close to 1), we can then extract an estimate for the strong phase of the 
penguin coefficient aKnK). The result is illustrated by the triangular-shaped area in 
the first plot in Figure 10. 

We stress at this point that, in the future, the analogous ratios in the B — > nK* and 
B — > pK systems will provide interesting information about tree-penguin interference 
in PV and VP modes. The corresponding ratios tfm in these two systems are expected 
to be larger than rFMi^^K) by roughly a factor of 2, thereby enhancing significantly the 
interference terms in (77) and the sensitivity to 7. Note that the sign of the interference 
term is expected to be different in the pK modes (Re(rFM) < 0) as compared with the 
7rX(*) modes (Re(rFM) > 0). 

The results for the penguin coefficients illustrated in Figure 10 motivated our sce- 
narios S3 to S4 explained earlier. The values of the penguin-to-tree ratios obtained in 
these scenarios are indicated by the dots on the solid contours (S3) and by the squares 
(S4). From Tables 2 and 3 we observe that these scenarios are rather consistent with the 
data within errors. The best match is obtained for scenario S4. A notable exception is 
the decay B^ — > 7r°K°, whose measured branching fraction exceeds the values obtained 
using QCD factorization. We will discuss below that this fact cannot be attributed to 
uncertainties of the factorization approach. 



5.1.2 Other decays 

Within the approximations described earlier the amplitudes for the decays B^ K^(f) 
and B^ — > K'^(f) coincide and are determined in terms of the combination ct^+ctl — ■^ct^^EW ^ 
where once again the terms with p — u can be neglected to a good approximation. The 
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experimental value of the ratio 



T(B- 



K- 



1.10 ±0.18, 



(78) 



is indeed consistent with unity. In QCD factorization we find the central values «q 



0.002, (i^ 



-0.038, and 



'2'-^3,EW ~ 



0.004, suggesting that the amplitude is dominated 



by the penguin coefficient a^- Under this assumption, and using the average of these 
two decay rates, we then obtain 



TpB^K 
^0 



(0) 



Fo^--(O) 



al{K(t)) 



K* 



T(B 



nl/2 



= 0.82 ±0.10. 



(79) 



Our theoretical value l.SS^Qigg for the left-hand side of this equation is in marginal 
agreement with the data. Leaving aside the possibility that the B ^ K form factor is 
much smaller than assumed in our analysis, this might indicate a more modest penguin 
enhancement in the PV modes B — > K(j) than in i? — > t:K* . 



5.2 Hints of departures from the Standard Model 

Some measurements related to the penguin-dominated AS* = 1 transitions are difficult to 
explain theoretically. We will now study the corresponding observables. The quantities 
we will focus on are the mixing-induced CP asymmetries in the decays B — > (pKs and 
B — > ri'Ks, and the CP-averaged B^ — > 7r°X° branching ratio. 



5.2.1 Mixing-induced CP asymmetries in B ^ 4>Ks and B rj'Ks decays 

In neutral S-meson decays into a CP eigenstate /, one can define a time-dependent CP 
asymmetry 

where Atub > denotes the mass difference of the two neutral B eigenstates. The 
notation B^{t) refers to a state that was B^ at time t — 0, and B^{t) refers to a state 
that was B^ at time t — 0. The identity of these states at i = is determined by tagging 
the "other B" , using the fact that and are produced in a coherent state at the 
T{AS). The quantity Sf is referred to as "mixing- induced CP asymmetry", while — C/ 
is the direct CP asymmetry. 

If all subdominant amplitudes are neglected, the mixing- induced CP asymmetries Sf 
for the three final states J/'4>Ks, (t>Ks, and r]'Ks all equal sin^^, where 0^ = 2/3 is the 
BcrBd mixing phase. As is well-known the correction to this result is smallest for the 
final state J /ijjKs- since the subleading amplitude is CKM-suppressed and suppressed by 
a penguin-to-tree ratio. The correction is somewhat larger for the other two final states, 
where one can rely only on CKM suppression. The measurements of Sf for / = (i)Ks 
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and r]'Ks will become more accurate during the coming years, so that it is interesting to 
quantify the correction. 

We can parameterize the Bd decay amplitude such that it is proportional to 1 + 
e~^^ df e^^f , where df > includes the small ratio | KidKibl / | VcdVctl of CKM elements, and 
6*/ is the strong- interaction phase difference between the two amplitudes with different 
weak phases. The CP-asymmetry Sf is then given by 

^ sin + 2df cosdf sin(0rf + ^y) + dj sm{(f)d + 27) 
^ l + 2dfCos6fCos^ + d'f 

— sin (f)d + 2df cos 9f cos (f)^ sin 7 -|- O(rfj) . (81) 
For the two modes of interest we obtain 

2d f COS Of = { 82 
' 1.8±0.4±0.4^?;^^°:^%; f = v'Ks. 

where the errors have the same meaning as in previous sections. The various scenarios 
give very similar values, ranging from 3.2-3.5% for (pKs and 0.5-1.7% for rj'Ks- 

As expected the correction is very small, since the QCD factorization approach does 
not contain any large enhancement mechanism that could compensate the CKM sup- 
pression. The theoretical uncertainties are under control in both cases, including those 
due to the power corrections parameterized by the quantities Xa and Xh- Using that 
cos(/)<iSin7 = 0.64 12 ?2 for 7 = (70 ± 20)°, we conclude that in the Standard Model 

SsKs - Smks = 0.025 ± 0.012 ± 0.010 , 

(83) 

S^'Ks - Sj/^Ks = 0.011 ± 0.009 ± 0.010 , 

where the second error estimates the uncertainty in Sj/^Ks- theoretical results 

obtained using QCD factorization agree with simple model estimates obtained in [51]. 
The correction for rj'Ks is smaller than for (pKs since the color-suppressed tree and 
penguin amplitudes partially cancel each other. The precision of the predictions (83) is 
higher than in a scheme where only SU(3) flavor symmetry is employed [52], in which 
case the two differences above can only be bound to be less than about 30%. The present 
experimental values of the two differences are 

ScliKs ~ ^J/ipKs — —1-11 ± 0.41 , 

(84) 

5^'i^,.-^j/^K, = -0.40 ±0.34. 

In the first case there is a 2.8o" discrepancy between the data and the theoretical pre- 
diction obtained in the Standard Model. If confirmed with more data, this would be a 
clean signal of physics beyond the Standard Model. 
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Figure 11: Theoretical prediction for the ratio i?oo as a function of 7. The 
central value is shown by the solid line. The inner (dark) band corresponds 
to the variation of the theory input parameters, while the outer (light) band 
includes in addition the uncertainty from weak annihilation and twist-3 hard- 
spectator contributions parameterized by Xa and Xh- The gray bands show 
the experimental result for Rqq with its la (dark) and 2a (light) errors. 



5.2.2 The large B° ir^K^ decay rate 

The experimental data for the tt^K^ decay rate are significantly larger than 

predictions obtained using QCD factorization. Sometimes this is interpreted as evidence 
for large rescattering effects. The purpose of this section is to point out that rescattering 
(or, more generally, any other source of hadronic uncertainty) cannot be invoked to 
explain the data. 

In many analyses the —>■ tt^K^ decay rate is normalized to the rate for B^ —>■ 7i^K~ 
in order to eliminate the dependence on the form factor [10, 48]. The resulting ratio is a 
strong function of 7, and in fact its high value can be explained by choosing a low value 
of 7. Here we propose to consider instead the ratio 



00 



T{B- 



TT 



-,exp 



(85) 



1.18 ±0.17. The 



0.79±0.02±0.06l[J:^?± 



of CP-averaged rates. The experimental value of this ratio is R, 
theoretical prediction obtained using QCD factorization is -Rqq 
0.04, essentially independent of 7. The largest uncertainties are due to the strange-quark 
mass, the S — >■ vr and B ^ K form factors, and Xa- The very weak dependence of this 
result on 7 is illustrated in Figure 11. The results obtained in the various scenarios are: 
i?th ^ Q 75 (^g^^^ Q_g7 (^g2), 0.81 (S3), and 0.86 (S4). Even the largest theoretical values 
are about 2 standard deviations away from the central experimental result. 
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The theoretical interpretation of the ratio it!oo in the Standard Model is very clean. 
Based on the amplitude parameterizations in (73) we define the ratios 



— 7: RttK 



0.12 -O.OH, 



fC — —RttK 



0.03- 0.02i, 



rp = 



0.02 



O.Oli, 



(86) 



where Rt^k = {.I-kI Ik) ■ {Fq^^/Fq^'') ~ 1- The numerical values correspond to the cen- 
tral results obtained in QCD factorization and are meant only to illustrate the smallness 
of the various ratios. The ratio tew determines the relative magnitude of electroweak to 
QCD penguin contributions. In the Standard Model, the theoretical prediction for this 
quantity is to a large extent free of hadronic uncertainties. To see this, note that tew 
can be written as the product of the two ratios | R^k Q;3,ew/ ('^1 + "^2) and (ai -\-a2)/6t'i. 
Using Fierz identities and top-quark dominance in the electroweak penguin diagrams the 
first of these ratios can be calculated in a model-independent way up to corrections that 
break V-spin (s ^ u) symmetry [47]. The second ratio can be determined experimen- 
tally using the strategies described in Section 5.1.1. Next, the quantity rc determines 
the ratio of the color- suppressed tree amplitude to the dominant penguin amplitude. 
This ratio is doubly CKM suppressed. The magnitude of the penguin coefficient al^ijiK) 
is determined by data, and the magnitude of the tree coefficient 0^2 is expected to be of 
order 0.2 (though it is larger in the scenarios S2 and S4). It then follows that rc cannot 
be larger than a few percent in magnitude. Finally, the ratio rp determines the rela- 
tive magnitude of up and charm-penguin contributions weighted by their CKM factors. 
Being again doubly CKM suppressed, this ratio is bound to be of order a few percent. 
Moreover, to first order this ratio cancels between numerator and denominator in the 
ratio i?oo- 

It is then justified to work to first order in the small ratios rc and rp and to neglect 
cross terms of order tew^c and tew^p- Parametrically, this amounts to neglecting 
terms of 0{\*, a\^, a^), where A ~ 0.22 is the Cabibbo angle and a the electromagnetic 
coupling. This gives 



where the neglected terms are safely below 1% in magnitude. For 7 70° the COS7- 
dependent term is about -1-0.02, which is negligible in view of the present experimental 

error m 

i?oo- The fact that R^^q < 1 can be understood in terms of a negative interference 
between QCD and electroweak penguin contributions. The value tew ~ 0.12 predicted 
by QCD factorization indeed leads to i?oo ~ 0-^9 in (87), which explains our numerical 
result. In the Standard Model the result i?oo < 1 is a practically model- independent 



Rqq = |1 - tewT + 2cos7Rerc + . . . , 



(87) 
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prediction. In this regard the theoretical analysis of the observable i?oo is almost as clean 
as that of the CP asymmetry S^Ks discussed in the previous section. If the experimental 
finding of i?oo > 1 is confirmed, this would be a clean signal of physics beyond the 
Standard Model. ^ 

The concerned reader may ask whether the conclusion just obtained could be affected 
by contributions to the decay amplitudes neglected in our approximation scheme. Using 
the exact expressions for the amplitudes collected in Appendix A but dropping again 



terms of 0(A , a\ 



a 



we find that the above analysis remains valid provided the ratios 



Tew and rc are redefined as 



rc ^ - 



3 R^K (^1,^w{t^K) + (5l■^^^{'KK) 
2 6il{T:k) 

R^K a2{nK) + P2{tiK) 



Am 



(88) 



The additional annihilation and electrowcak annihilation contributions are so small that 
they cannot affect the analysis presented above. 

Since (87) contains a term linear in the parameter tew, one might try to increase the 
value of Rqo (without changing much the decay rates for other B ttK modes) by adding 
an 0(1) New Physics contribution to the electroweak penguin coefficient al^^T^{7rK). 
Such large non-standard effects with the flavor structure of electroweak penguins can 
indeed arise in a large class of extensions of the Standard Model [54]. However, there 
is an alternative way of displaying the "vr'^i^^ anomaly" which points in a different 
direction. Using the isospin properties of the effective weak Hamiltonian and of the 
(ttK) final states it can be shown that the ratio 



Rt 



T{B- 7r-i?o) + r(BO ^ TT+K-) 



(89) 



equals 1 up to corrections that are at least quadratic in small amplitude ratios [55, 56]. 
In terms of the amplitude ratio tew defined in (86) and a new ratio 



Am 



is) 



0.18-0.02i, 



we find 



Rl = 1 + |rEwP - cos7Re(rTrEw) + ■ ■ ■ 



(90) 



(91) 



where we have neglected smaller second-order terms involving the ratios rc and vp. 
Because the corrections arc of second order in small ratios the theoretical expectation 



''The conclusion that the current wK data may point to an electroweak penguin contribution at 
variance with the Standard Model value of has been reached independently in [53] through an 

analysis of a larger set of ratios. 
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for Rl lies very close to 1. We find Rf = 1.01 ± 0.01 ± 0.01 1^^^ ± 0.01 using our default 
error analysis, while the different scenarios predict: R^^ — 1.03 (SI), 1.02 (S2), 1.01 
(S3), and 1.02 (S4). The current experimental result, R^^^ = 1.24 ± 0.10, deviates from 
the theoretical expectation by about 2o", as recently emphasized in [57]. To account 
for the experimental value requires that the magnitude of tew be many times larger 
than in the Standard Model (since rx is restricted by data to lie close to its theoretical 
value). Specifically, from (87) and (91) it follows that in order to reproduce the central 
experimental data one needs tew ~ ±0.5i with a large (weak or strong) phase. 

Before proceeding, we stress that ratios analogous to -Rqo and Rl can also be consid- 
ered for vector-pseudoscalar final states, for which the final states ttK are replaced 
by pK or nK*. Let us briefiy consider the case of i?oo in more detail. To lead- 
ing order in small quantities relation (87) holds also in these cases once the ratio 
R^j^ ^ 1.0 is replaced by corresponding ratios Rp^ = {fp/ fx) ■ {F^^^ / A^^'^) fa 1.2 
and Rt,k* = {fn/ Ik*) ■ {A^^^* / F'^^'^) ^ 1.0. Due to the suppression of the penguin 
coefficient a1 the vahies of the ratios Tew and are larger for the PV modes than for 
the PP modes. In addition, the sign of tew in the pK modes is negative and opposite 
to the cases of t^K and t^K*. This qualitatively explains our results R^^{pK) — 1.73 ^q'S 
and R^^{'kK*) — 0.411^q;26- Obviously the uncertainties are significantly larger than in 
the TiK case, but nevertheless it would be interesting to probe the qualitative features 
of these predictions. Unfortunately, the neutral B decays into p^K^ and ti^K*^ have not 
been seen yet experimentally. However, in the case of t:K* we may derive the experi- 
mental upper bound R^^^{'kK*) < 0.9 at 90% confidence level, which agrees with the 
prediction that this ratio should be smaller than 1. 



5.3 Penguin and annihilation-dominated AD = 1 decays 

The last set of decays from which interesting information on penguin and annihilation 

terms can be deduced is i? ^ KK and the corresponding PV modes. The decay ampli- 
tudes for 5° — > K^K^, K^K*^, and K*^K^ are particularly interesting, 
since they are the only decay modes that receive only weak annihilation contributions. 
The corresponding simplified amplitudes are given by 



Af 



*K-K° 



^A 



KK 



5pu P2 + &4 



A 



B0^K-K+ 



A 



KK 



A 



BO^KOKO 



KK 



(92) 



where it is understood that each term must be multiplied with Ap"^^ and summed over 
p = u, c. The expressions for the B KK* and B K*K amplitudes are obtained 
by replacing (KK) (KK*) and {K*K), respectively. We will also consider the modes 
B — > 7r0, whose amplitudes are extremely simple and given by 



Af 



■V2A 



A^s oc 



3 ■ 



(93) 
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Our default theoretical predictions for these modes are shown in Tables 4, 5 and 6. 

In principle, the modes — > K~K'^ (and the corresponding PV modes) can provide 
interesting information about the weak annihilation amplitudes. Unfortunately, however, 
their branching fractions are too small to be observed in the near future. The other decay 
modes are penguin dominated. Using this information we predict 

T(B- K-K^) T(B- K-K*^) T(B- K^K*-) , , 

The relevant branching fractions are found to be of order few times 10~^, which should 
be within the long-term reach of the B factories. If these equalities could be established, 
this would be another indication that annihilation contributions are suppressed with 
respect to the dominant penguin amplitudes governed by a^. 

In the approximations described above, the decays B~ — > 7r~(/) and B^ — > 7r°0 are 
determined by the singlet penguin coefficients af, which are predicted to be highly 
suppressed in the QCD factorization approach. In addition these modes are CKM sup- 
pressed. An experimental upper bounds exists for B^ 7r^0, which however lies two 
orders of magnitude above the theoretical predictions. It will therefore not be possible 
to extract useful information from these modes is the foreseeable future. 



5.4 Bounds on large annihilation contributions 

Our discussion in this section was confined to the treatment of annihilation effects sug- 
gested in [10], in which the value of the model parameter qa is limited to qa < 1, 
corresponding to a 100% uncertainty on the default value for the quantity Xa- While we 
still believe that this is a reasonable assumption, which so far has not been invalidated 
by the data, one might ask what would happen if larger values of qa were considered. 
Note that the annihilation kernels in (57) and (58) are quadratically dependent on the 
quantity X^, so that increasing the value of qa can have a dramatic effect on the relative 
strength of annihilation terms as compared with the leading penguin amplitude. 

It is apparent from Figure 10 that for values of qa significantly larger than 1 a fine- 
tuning of the phase (fA is required so as not to be in confiict with the experimental 
constraints on the magnitudes of the parameters in the tiK and ■nK* systems. Let us 
assume for a moment that this fine-tuning is indeed realized in nature. To be specific, 
we set = 2 and adjust the phase ipA — —60° such that the resulting 0.% values fall in 
the center of the dark rings in Figure 10 (see the dots on the dashed curves). All other 
parameters are set to their default values. We might then ask whether such a scenario 
is consistent with the experimental data. 

Interestingly, we find that for several decay channels the results obtained with such 
large annihilation contributions are already in confiict with the data. In the upper part 
of Table 7 we collect results for branching ratios where significant discrepancies appear. 
In the lower portion of the table we give results for K^*^K^*^ modes, some of which 
are also very sensitive to annihilation effects. We conclude that the data support the 
assumption that a reasonable estimate of weak annihilation effects can be obtained by 
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Table 4: CP-averaged branching ratios (in units of 10 ^) of penguin- dominated 
B — > PP decays (top) and B — > PV decays (bottom) with AD = 1. 



B- 
B^ 



Mode 



Theory 



1.36 
1.35 



hO.45 +0.72 +0.14 +0.91 
-0.39-0.49-0.15-0.40 
h0.41 +0.71 +0.13 +1.09 
-0.36-0.48-0.15-0.45 



SI 

2.20 
2.12 



S2 
1.52 
1.53 



S3 
1.71 
1.83 



S4 

1.46 
1.58 



Experiment 

< 2.2 

< 1.6 



B- 
B^ 
B- 
B° 



. K-K*^ 
. K^K*- 



0.30^2: 



0.26 
0.30 
0.29 



+0 
-0, 
+0, 
-0, 
+0, 
-0, 



11+0.12+0.09+0.57 

09-0.10-0.09-0.19 
08+0.10+0.08+0.46 
07-0.09-0.08-0.15 
08+0.41+0.08+0.58 
07-0.18-0.07-0.17 
10+0.39 +0.08 +0.60 
09-0.17-0.07-0.17 



0.50 0.21 0.54 0.66 

0.41 0.18 0.45 0.56 

0.45 0.75 0.61 0.55 

0.47 0.71 0.61 0.54 



< 5.3 



B- 



► TT (p 



0.005 
0.002 



0.008 0.010 0.005 0.009 
0.004 0.005 0.002 0.004 



< 0.4 

< 5.0 



Table 5: Direct CP asymmetries (in units of 10^^) of penguin-dominated B 
PP decays (top) and B — > PV decays (bottom) with AD = 1. We only 
consider modes with branching fractions larger than 10~^. 



Mode 


Theory 


SI S2 S3 S4 


Experiment 


B- K-K° 


ifi Q+4.7+5.0+1.6+11.3 

-•-"■"J -3.7-5.7-1.7-13.3 

Ifi 7+4.7+4.5+1.5+4.6 
J-O-' _3.7_5.i_i.7_3.6 


-10.0 -14.7 -6.5 -4.3 
-10.6 -15.0 -12.7 -11.5 




B- K-K*^ 
B- K^K*- 


OQ cr +6.9+7.8+5.5+25.2 
^•-•■^ -5.7-9.0-6.5-36.8 
OK 7+7.4+7.2+5.7+10.9 
-5.7-9.0-6.9-13.4 
-.0 4+3.7+7.8+4.2+27.4 
-•-"^■^-3.0-3.5-4.7-36.7 
-1 q 1 +3.8+5.4+4.5+5.8 
-3.0-2.9-5.2-7.4 


-13.9 -22.5 -7.7 -9.6 
-16.9 -26.0 -16.8 -13.6 
-8.9 -12.7 -13.9 -21.1 
-8.0 -12.2 -8.8 -10.0 





Table 6: CP-averaged branching ratios (in units of 10 ®) of annihilation- 
dominated B PP decays (top) and B PV decays (bottom) with AD = 1. 



Mode 


Theory 


SI S2 S3 S4 


Experiment 


S° ^ K-K+ 


A niQ +0.005+0.008+0.000+0.087 
U.UlO _o.o05 -0.005 -0.000 -0.011 


0.007 0.014 0.079 0.070 


< 0.6 


S° ^ K-K*+ 
B° K+K*- 


n ni4+0-007 +0.010+0.000 +0.106 

U.Ul^ -0.006 -0.006 -0.000 -0.012 

n ni/i +0.007+0.010+0.000+0.106 

U.Ul^ -0.006 -0.006 -0.000 -0.012 


0.016 0.011 0.095 0.094 
0.016 0.011 0.095 0.056 
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Table 7: CP-averaged branching ratios (in units of 10~®) of some decays that 
imply bounds on the weak annihilation parameter X^. The theoretical results 
refer to the choice = 2 with strong phase ipA = —60°. 



Mode 



Default Large Annihilation 



Experiment 



5° 

B- 

B- 
5° 



K-p+ 
► K-(f) 



0.7 
2.6 
7.4 
4.5 
4.1 



6.0 

9.0 

19.3 

22.4 

20.2 



< 3.6 

< 6.2 
8.5 ±2.1 
9.2 ±0.7 
7.7 ± 1.1 



B- 
B- 
B- 

5° 



• K*-K^ 
■ K-K*° 
K-K+ 
K*-K+ 
K-K*+ 



1.36 
0.30 
0.30 
0.01 
0.01 
0.01 



1.65 
1.14 
0.95 
0.21 
0.26 
0.26 



< 2.2 

< 5.3 

< 0.6 



varying the parameter qa between and 1, allowing for an arbitrary strong phase ipA- 
Values of Qa larger than 1 are strongly disfavored by the data, and values > 2 can 
already be excluded, at least if universal annihilation is a reasonable approximation. It 
has occasionally been argued that one should allow for much larger values of < 8 to 
account for the theoretical uncertainties related to power corrections [58] . Table 7 shows 
that such large weak annihilation effects cannot be tolerated by the data. 



6 Tree-dominated decays 

Many of the decays with AD = 1 are dominated by tree amplitudes and do not suffer 
from a large sensitivity to light quark masses or chirally-enhanced power corrections. We 
now analyze this class of decays, including the time evolution of the 7rp final states. 

We first summarize the decay amplitudes simplified according to the discussion of 
Section 4.2. There are three independent amplitudes for the decays S — > 7rp given by 



-2 ABO-^nOpO 



A 



■Kp 



Spu (Qi2 - P2) - Oil 



A 



pn 



B°->-TT+p- — -^np 



'■TTp 



^pu 

5pu (a2 - A) - 6^ - 2/3f] + Ap, {a^ - A) - - 2/3^] . 



(95) 
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Table 8: Magnitudes of the amplitude parameters for tttt and up final states 
in the different scenarios. "Default" refers to the default input parameters. 
Scenarios not shown agree with the default. A minus sign indicates that a 
parameter is negative in naive factorization or (in case of the annihilation 
parameters) ii qa — ^- Theoretical errors are added in quadrature. 





Ueiault 


Co OA 




Deiault 


Co Co OA 

bZ b6 d4 


ai{Txp) 
ai{pn) 


n qq+O-04 
n qq+0.04 


0.84 0.88 
0.89 0.90 
0.94 0.95 


(3i{7rp) 
fSiiprc) 


025+0-046 
020+0-037 
n 024+0-045 


0.032 0.062 0.071 
0.023 0.051 0.063 
0.024 0.061 0.053 


a2(7r7r) 
tt2(7rp) 
a2{p7i) 


o.2ol°:lI 

U--LO_0.09 


0.57 0.48 
0.45 0.41 
0.31 0.30 


-/32(vr7r) 
-/32(vrp) 
-/32(p7r) 


o.oiotHS? 

0.008^HS^ 
009+0-018 


0.012 0.024 0.027 
0.009 0.019 0.024 
0.009 0.023 0.020 


-a2{7ip) 
al{p7r) 


n 1 02+0-020 

U.UOZ_Q QQg 

n 035+0-020 
U.UOd_o.oi3 


0.110 0.102 
0.027 0.028 
0.062 0.049 


-/?3(7r7r) 


009+0-032 
5+0-024 

U.UUO_o.004 

007+0-031 


0.013 0.035 0.041 
0.007 0.025 0.034 
0.008 0.032 0.028 



These expressions must be multiplied with X^p^ and summed over p = u,c. The order of 
the arguments of the coefficients (M1M2) and /3f (M1M2) is determined by the order of 
the arguments of the AmiM2 prefactors. The amplitudes for B~ — > 7rOp~ and ^o 7r~p+ 
are obtained from the first two expressions by interchanging n ^ p everywhere, including 
the arguments of the amplitude parameters. The expressions for the B tttt amplitudes 
are obtained by setting p — > tt, in which case they simplify to 



V2A 



B-- 
AbO- 



A 



TTTT '-'pu ' 

^pu 



(96) 



-A 



BO- 



rrO^O — At, 



Bose symmetry implies that the two-pion final state must have isospin 7 = 0,2 but not 
I — 1, which is the reason why the S — > tttt decay amplitudes are simpler than those for 
B — > np. Finally, 



V2Af 



-2 A 



A^, 



OLa 



Spu ("2 + P2) + 2a| 
Spu (o;2 - /3i) + 2q;| + al 



A^ 



+ A 



Spu (ttl + P2) + 



a'. 



Spu {-(X2 - Pi) + Ol\ 



(97) 



To understand the pattern of branching fractions and asymmetries it is useful to bear 
in mind the magnitudes of the amplitude parameters as given in Table 8. The table also 
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shows the values of the parameters in the four scenarios defined in Section 5. Generically 
a2 is large in scenarios S2 and S4 (and is somewhat reduced), while the annihilation 
coefficients are increased in S3 and S4. Because of the cancellation of vector and scalar 
penguin contributions, the combination a\{pn) = a^{p7i) — r'^aKpn) is very sensitive to 
the strange-quark mass. Except for the purely neutral final states the most important 
interference of amplitudes with different weak phases occurs between the tree coefficient 
CKi and the effective penguin amplitude a^. Comparing the Trp final states to tttt we 
note two important differences: first, the penguin amplitudes are smaller for irp. This 
implies reduced sensitivity to the angle 7 in CP-averaged branching fractions and smaller 
direct CP violation, but it also implies a reduced "penguin pollution" in time-dependent 
studies of CP violation. Second, while al{n7i) and Q;i(7r7r) have opposite signs (always 
assuming small relative phases) , which implies constructive tree-penguin interference for 
7 < 90°, both relative signs occur for np. 

6.1 Branching fractions and direct CP asymmetries 

Our results for the tttt, irp, and ttcu branching fractions and direct CP asymmetries are 
shown in Tables 9 and 10. In many cases the dominant theoretical uncertainty arises 
from the variation of CKM parameters or form factors. The latter source of uncertainty 
can be reduced once better data on semileptonic or leptonic B decays become available. 
The sensitivity to CKM parameters is not a theoretical hmitation but rather provides 
access to \Vub\ and 7. 

Prom the numbers shown in the tables one observes reasonable global agreement be- 
tween theory and data within their respective error ranges. Since many of the theoretical 
errors are correlated among the various final states we shall consider below certain ratios 
of observables, which provide more sensitive probes of the underlying hadronic ampli- 
tudes. With default parameters the decay rates to the charged final states tt+tt" and 
TT^p"^ are predicted to be significantly larger than the data. The discrepancy disappears 
for the two-pion final state if 7 is large (scenario SI) or a2 is large (scenario S2). How- 
ever, the table shows that the up modes do not favor the large-7 scenario. Both S2 and 
S4, which combines elements of the large-Q;2 scenarios with a moderate increase of weak 
annihilation, describe the branching fractions very well, but S4 is singled out if we add 
the information from penguin-dominated decays discussed in the previous section. 

The data on direct CP asymmetries are still too uncertain to draw conclusions from 
the comparison with theory. We note a 2-3 standard deviation discrepancy for the 7r"'"7r~ 
final state. The current central value 0.51 of the CP asymmetry is certainly too large 
to be understood in the QCD factorization framework, and appears even more puzzling 
in view of the small asymmetries for the ttK final states. Unless we allow for very 
large SU(3)-breaking effects the only difference between tt+tt" and 7r~^K~ occurs in the 
annihilation term Spuf^i + 2/54 , but it is hard to see how this could cause such a dramatic 
effect on the CP asymmetry, since (3^ is always very small in our framework and never 
exceeds a quarter (and often less) of the other penguin annihilation coefficient /3f . In 
the future this combination of annihilation amplitudes can be constrained independently 
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Table 9: CP-averaged branching ratios (in units of 10 ^) of tree-dominated 
B ^ PP decays (top) and B ^ PV decays (bottom) with AD = 1. The 
errors and scenarios have the same meaning as explained in Section 5. 



Mode 



Theory 



SI S2 S3 



S4 



Experiment 



B — >• TT TT' 







5° 



^0 ^ ^0^0 



fy n +3.0+2.1 +1.0+0.4 
U.U -2.4-1.8-0.5-0.4 
o Q+4.0+3.6+0.6+1.2 
-3.4-3.0-1.0-0.8 
n o +0.2+0.2+0.3+0.2 
'-'■'^-0.2-0.1-0.1-0.1 



5.8 

6.0 
0.7 



5.5 
4.6 
0.9 



6.0 
9.5 
0.4 



5.1 
5.2 
0.7 



5.3 ±0.8 

4.6 ±0.4 

1.6 ±0.7 
(< 3.6) 



B~ 
B~ 

B^ 



7rOp° 



-6.3+3.6+2.5+1.3 
5.0-3.1-1.2-1.1 
-6.5+5.1+1.0+0.8 
-5.5-4.3-0.6-0.7 
10.3+8.7+1.3+2.0 
8.4-7.2-2.3-1.6 



11.9 
14.0 

21.2: 

1 p; /I +8.0+5.5+0.7+1.9 
-'-'-'•^-6.4-4.7-1.3-1.3 

36.5 
0.4 



14.2 12.6 
10.7 10.4 



+18.2+10.3+2.0+3.9 
-14.7- 8.6-3.5-2.9 
0.2+0.2+0.9+0.5 
0.2-0.1-0.3-0.3 



18.6 
17.5 

36.1 
0.3 



11.0 
10.8 
21.8 
1.7 



12.2 
14.2 
22.2 
16.4 
38.6 
0.3 



12.3 
10.3 
11.8 
11.8 
23.6 
1.1 



9.1 ± 1.1 
11.0 ±2.7 
13.9 ±2.7 
8.9 ±2.5 
24.0 ±2.5 
< 2.5 



B- 
B' 



I +4.4+2.6+1.8+0.8 
' -3.5-2.2-0.9-0.9 



ni+ooo- 

-0.00- 



-0.02- 
-0.00- 



-0.02- 
-0.00- 



-0.03 
-0.00 



8.6 9.1 8.4 8.4 
0.01 0.07 0.01 0.01 



5.9 ±1.0 
< 1.9 



Table 10: Direct CP asymmetries (in units of 10~^) of tree-dominated B PP 
decays (top) and B —>■ PV decays (bottom) with AD — 1. We only consider 
modes with branching fractions larger than 10~^. 



Mode 



Theory 



SI 



S2 



S3 



S4 



Experiment 



B- 
B^ 



77 77 



77 ' 77 



77O77O 







-0.02- 



hO.Ol +0.05 +0.00 +0.01 
-0.01-0.05-0.00-0.01 
—6 5+2-1-^3.0+0.1+13.2 



-2.1-2.8-0.3-12.8 



45.1 



-18.4+15.1- 
-12.8-13.8 



4.3- 
14.1- 



-46.5 
61.6 



-0.02 -0.02 -0.02 -0.02 
-9.6 -9.1 5.6 10.3 
23.0 21.7 5.6 -19.0 



-7± 14 
51 ±23 



B- 
B- 
B' 
B' 
B' 



► 77- p" 

>77V 

77~p+ 
77V° 



A 1 +1.3+2.2+0.6+19.0 
-0.9-2.0-0.7-18.8 



-4.0 
-1.5 



+1.2+1.8+0.4+17.5 
-1.2-2.2-0.4-17.7 
0.4+1.2+0.2+8.5 
0.4-1.3-0.3-8.4 
n +0.2+1. 3 +0.1+11. 5 
"-'•^-0.1-1.6-0.1-11.7 
1 c 7+4.8+12.3+11.0+19.8 
-10.1 _4.7-i4.0-12.9-25.8 



3.4 
-5.3 

-1.7 
0.5 
-20.9 



4.6 
-6.3 
-1.8 

1.7 
-9.5 



-13.3 
12.2 
6.6 
-10.3 
-10.6 



-11.0 

9.9 

3.9 
-12.9 

10.7 



-17± 11 
23 ±17 
-11 ± 17 
-62 ± 27 



B- 



77 cu 



_i 0+0.5+2.7+0.8+2.1 
-0.5 -3.3 -0.7 -2.2 



0.6 



-2.1 -6.0 



9 ±21 
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using — > K~K'^ decays (up to SU(3)-breaking effects). 

Witfi tfie exception of the neutral final states the direct CP asymmetries are all 
predicted to be small, although with large uncertainties due to weak annihilation that 
affect even the signs of the asymmetries. However, the analysis of Section 5 has shown 
that for the largest values of weak annihilation allowed by our error definition not all 
values of the annihilation phase are compatible with data on branching fractions and CP 
asymmetries of penguin-dominated decays (see Figure 10). This might favor the signs 
displayed in scenario S3 and S4, although this conclusion must be regarded with great 
caution since it is mainly based on the direct CP asymmetry in the n^K"^ mode. 

It is worth emphasizing that the correlations between various asymmetries are pre- 
dicted more reliably in the factorization framework than the asymmetries themselves, 
at least if the annihilation parameter qa is not too different for the PP, PV, and VP 
amplitudes, or if annihilation is a small effect altogether. Table 10 shows that this is the 
case for all our scenarios. Since in the case of a small penguin-to-tree ratio the direct 
CP asymmetry is approximately given by 2sin7 \\^c'' /y^u^\ Ini(a5/ai), the relative mag- 
nitudes of the CP asymmetries in tt+tt^, 7r^p°, 7r°p^, 7r^p+, and vr+p^ provide direct 
access to the magnitudes and signs (phases) of the penguin amplitudes al for tttt, 7rp, 
and pTT, which are predicted to be rather different (see Table 8). We note specifically the 
case of the -k' and B^ vr^p" amplitudes, where according to (95) the PV 

and VP penguin amplitudes appear as — At^p al but with opposite sign relative 

to ai, and the charged final states, which probe the magnitude and sign of the PV and 
VP penguin amplitudes independently. A verification of these sign patterns alone would 
provide an impressive confirmation of the relevance of factorization to the calculation of 
direct CP asymmetries. 



6.2 Ratios of decay rates 

We now discuss a number of ratios of CP-averaged tttt and 7rp decay rates, which are 
either sensitive to the CKM angle 7 or to particular aspects of the hadronic amplitudes. 
In addition to Rt^t^ = V{B^ — >• 7r"'"7r~)/(2r(i^~ — >• 7r~7r°)) we consider the ratios: 



1(5° ^ 






r(so ^ 


TT+TT' 


-)' 


r(50^ 




-) 


r(so ^ 


7r~p' 




2r{B- - 







r(so ^ TT+p-) + r(so ^ 7r-p+) 



2r(S0 ^ TT+TT-) 

^3 = ; 'i-'^J , (98) 

p - --V- r . . 2r(i?-^7rV-) 

r(so^7r-p+) ' r(so^7r+p-) 

The theoretical and experimental results for these ratios are summarized in Table 11. 
The ratio i?7r7r has already been discussed in part in Section 4.3 as one of the motivations 
for exploring a scenario with large a2- Since QCD factorization does not contain a 
mechanism to generate a large strong phase between the two amplitudes with different 
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Table 11: Results for the ratios R^^^ and R1-5. 



Theory 



SI S2 S3 S4 



Experiment 





0.80 


Ri 


2.39 


R2 


2.06 


R3 


1.38 


i?4 


0.42 


i?5 


0.22 



+0.12+0.05+0.13+0.14 
-0.12-0.06-0.20-0.11 
+0.31+0.04+0.1.5+0.05 
-0.25-0.08-0.12-0.11 
+0.40+0.53+0.12+0.03 
-0.30-0.36-0.09-0.06 
+0.18+0.82 +0.03+0.02 
-0.17-0.59-0.04-0.05 
+0.04 +0.15 +0.45 +0.23 
-0.04-0.11-0.21-0.20 
+0.07+0.08+0.23+0.14 
-0.08-0.06-0.12-0.12 



0.56 
3.11 
3.02 
1.06 
0.49 
0.06 



0.45 
2.39 
2.37 
1.01 
1.14 
0.75 



0.85 
2.33 
2.03 
1.35 
0.37 
0.18 



0.55 
2.28 
2.28 
0.99 
0.92 
0.61 



0.47 ±0.08 
3.02 ±0.64 
2.61 ±0.35 
1.56 ±0.53 
0.88 ±0.57 
0.46 ±0.46 



weak phases contributing to — > tt+tt", Rt^tt is described well only if 7 is significantly 
larger than 100°, or if a2 is large and 747r7r|Ai''^| is small. 

We find that the ratio i?i is theoretically rather clean and yet sensitive to 7, because it 
is almost independent of the 5-meson form factors, and it is independent of the uncertain 
color-suppressed tree coefficient 02 as seen from (95). In fact, if the tree coefficient ai 
dominated the amplitude and were universal, Ri would equal (/p/Ztt)^ = 2.55, which is 
not far from the complete result. In our framework the largest theoretical uncertainty 
comes from the pion and p-meson light-cone distribution amplitudes, since these cause 
the largest non- universality of the amplitude parameters. From Table 11 we see that Ri 
is nearly constant, except in the scenario where 7 = 110°. The sensitivity to 7 arises 
from the fact that while in both cases the penguin-tree interference is constructive for 
7 < 90°, the effect is more pronounced for tttt than for irp due to the larger penguin 
amphtude in the former case. We should note, however, that the conclusion that Ri is 
theoretically clean hinges on the assumption that the phase of the penguin annihilation 
amplitude is not very different for the PP and PV amplitudes, or, if it is different, that 
the magnitude of penguin annihilation is not too large. The current experimental value 
of Ri is somewhat higher than the theoretical prediction for 7 = 70°, but the difference 
is only one standard deviation. 

Since the sum of the 7r~^p~ and 7r~p~^ decay rates is measured more accurately than the 
individual rates, we also consider the second ratio R2. In this case there is a substantial 
uncertainty from the ratio of form factors, Aq^'' (0) / F^^"^ (0) , evidenced by the second 
error of the "Theory" column in Table 11. Perhaps the only thing that can be said at 
present is that theory and experiment agree within their respective errors, which can be 
taken as a qualitative argument in favor of factorization. 

The ratio R^ of n^p^ to B^ 7r~ p'^ is mainly sensitive to the form-factor 

ratio A^^''(0)/Fq^^""(0), and to a lesser extent to 7. Since the information provided 
by i?3 is largely equivalent to the one from the parameter AC in the time-dependent 
analysis of the n^p^ final states (see Section 6.3 below), we do not discuss R3 further 
here. The experimental result is again in good agreement with factorization. 
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Finally, the quantities and R5 are constructed such that the dominant tree am- 
plitude cancels out in the numerator. This leaves as the leading term 

where the dots denote the interference terms with the penguin amplitude and other 
smaller contributions. For a similar approximation (with somewhat larger correc- 
tions) holds with TT <-> p exchanged everywhere. Hence i?4 5 are expected to be nearly 
independent of 7, and in first approximation they access the real part of the color- 
suppressed tree amplitude. These observables may therefore be interesting for assessing 
the viability of the large-Q;2 scenarios S2 and S4. However, the current experimental 
values of R4 and R5 are not sufficiently precise to support or disfavor these scenarios. 



6.3 Time-dependent rates in the n^p^ system 

We now analyze the asymmetries in time-dependent measurements of B'^ and 5° decays 
into the vr^p^ final states. We begin with setting up conventions and notation for a 
general final state / and its CP conjugate / [59]. 

We define Af, Af, Aj, and to be the amplitudes for the four decay modes, where 
the bar on A refers to the decay of the S° meson. In the standard approximation, which 
neglects CP violation in the B^-B^ mixing matrix and the width difference of the two 
mass eigenstates, the decay amplitude squared at time t of the state that was a pure B^ 
at time t = can be parameterized by 

\Am' ^ \{fm))? = ^(l^/r + {l + C^cos(Am^i) - 5^sin(Am^i)} , 

(100) 

where Ams > denotes the mass difference, and P the common total width of the 
5-meson eigenstates. For an initial B^ the signs of the cos and sin terms are reversed. 
For decays to the CP-conjugate final state one replaces / by /. 

For the following discussion we adopt the phase convention CP|i?°) = —\B^) and 
define the amplitude ratios 

(101) 

In terms of these 

which are phase-convention independent. (The CKM angle (3 is defined according to 
the convention of the Particle Data Group.) The system of four decay modes defines 
five asymmetries, C/, 5*/, Cj, Sj - alternatively parametrized as C = | {Cf + Cj), 
S = \{Sf + Sf), AC = ^ {Cf - Cf) and A^ = | {Sf - Sf) - together with the global 
charge asymmetry related to the overall normalization of (100): 

1+Acp _ \Af\'+\Af\' 



1 - Acp \Aff + l^/l^ 



(103) 
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or 

\A,\- + \A,\--\Aj\'^-\Aj\- 
- \A,\^+\A,\^ + \Aj? + \Aj? ■ 

Under a CP transformation pf goes into l/p/; hence AC and A^" are CP-even, but ^cp, 
C, and S arc CP-odd. For the special case that / is a CP eigenstate there are only two 
different amplitudes since / = /, and Aqp, AC and AS* vanish. 

For the np system we identify / = 7i~p'^. We write the amplitudes as 

•^f ~ Ati-p+ — Ap^ I'^i I ^pTT ( 6 ~^ Tj, j ) 

(105) 

Aj = = - |AW| T., (^e^^ + . 



For / = 7r+p~ one must exchange tt p on ^4, T, and P. These equations define the 
"tree" amplitudes T and the "penguin-to-tree" ratios, while the factors A^^p and Apj^ 
are given by (11). From (95) wc sec that An-p+ is proportional to Apj^ up to the small 
weak annihilation contribution dpuPi + /Sf, which suggests to extract this factor as we 
have done. The sign in the last amplitude is related to the CP convention for the B- 
meson state and the convention that all decay constants and form factors are taken to 
be positive.^ 

^The phases of the S-meson decay amphtucles are obtained as foUows. Lot C \n^) = £,it\t^^) and 
C \p~) = ^p\p'^)- Our CP convention for the B meson imphes C\B°) = |-B°), i.e., = 1- Isospin invari- 
ance suggests to choose Cir = 1 and = —1, since tt" is C-even and is C-odd, but in the following we 
shall keep the charge-conjugation phases for the charged pions and p mesons arbitrary. It follows from 
this that |7r-7r+) is CP-even, but CP |7r-(p)p+(-p)) = ^^Cp\^+ i~p)p- (p)) = -iA;k+{p)p- i^p)}- 
The second equality assumes that we form a wave-packet state by integrating over the relative momen- 
tum p and uses that the tt~ p'^ is in a P-wave state. Using CPO (CP)^ = for the operators in the 
effective weak Hamiltonian and dropping the momentum labels that are understood to be equal on the 
left and right-hand sides, we obtain 

{n-7r+\O^B°} = -(tt-tt+IOI^O) , {7T-p+\O^B°) = ^^Cp {7r+p-\0\B°) . 

This can be used to show that /9^+7r- = — e~^*'^, and hence 3^^^^ = sin2Q: in the absence of the penguin 
contribution. Similarly, one finds S^OpO = sin 2a for the neutral np final state. 

The conclusion is more subtle in the case of a non-CP eigenstate such as tt^p^, since after applying the 
CP transformation p^^- p+ involves the B matrix elements of two different final states. To clarify the issue 
we evaluate the matrix elements in naive factorization before and after applying the CP transformation, 
neglecting in addition the penguin contribution. This gives 



-2i7 _ -2i7 t J ^0 



The standard C and P transformations of the vector and axial-vector currents together with the CP 
convention for the B meson imply f„+ = fp+ = -^pfp-, F^^'^ = -^*F^^'^*, and Aq^'' = 

£,pAq^^ . The previous equation is consistent with this. However, the relative phase (actually, only a 
sign) of /b times the form factor to the light-meson decay constant cannot be inferred from symmetry 
considerations alone. Neither can it be determined from the semileptonic decay rates, which provide 
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Table 12: Magnitudes of the amplitude parameters in B ^ tttt and B ^ np 
decays. 





Theory 


S2 


S3 


S4 


IT 1 

1 TTTT 1 


n qi+0.05 

'-'•^-'--0.07 


0.75 


0.92 


0.81 


1 -Pnn / -^TTTT 1 




0.49 


0.37 


0.48 


IT 1 
1 1 




0.88 


0.99 


0.88 


1 -^np/^np 1 


"-'•-'-"-'-0.04 


0.12 


0.12 


0.20 


IT 1 

1 P""" 1 


-"-•"J'-o.o? 


1.03 


1.11 


1.07 


1 -PpTrf^pn 1 


"-'•-'-"-'-0.05 


0.19 


0.14 


0.15 



The magnitudes of the "tree" and "penguin-to-tree" parameters in QCD factorization 
are given in Table 12, where we also show the corresponding quantities for the tttt case. 
Note that the sign of Ppn/Tpj^ is negative in naive factorization, where all amplitudes are 
real, while the other two penguin-to-trcc ratios arc positive. The largest uncertainty on 
the penguin-to-tree ratios is caused by \Vuh\-i iTig, and weak annihilation. The important 
point is that the penguin amplitudes are smaller for the 7rp final states compared to 
TTTT. Furthermore, I-Pttp/T^pI and \PpTr/TpT^\ are similar in magnitude and interfere with 
an opposite sign with the tree amplitude. We should note that the near equality of the 
two penguin-to-tree ratios is the result of intricate dynamics specific to factorization. 
\P-irp/T.,rp\ is explained by the fact that al{7ip) + r^ag(7rp) ~ aKirp), since the scalar 
penguins have a small effect on PV amplitudes. On the other hand, for \PpK/TpT^\ the 
scalar penguin amplitude OglP"^) large and determines the sign of 04(^77) — rjag(p7r), 
such that the PV and VP amplitudes have similar magnitude but opposite sign. 

The results for the time-dependent asymmetry parameters are given in Table 13 
including again our default prediction with errors and the four standard scenarios to 
exhibit possible correlations. Since the penguin-to-tree ratios are small, it is instructive 
to compare the complete results with an expansion of the asymmetries in these ratios. 

access only to the magnitudes of these quantities. The expression after the second equality makes it 

clear that wc cannot simply assume all decay constants and form factors to be positive numbers. 

To fix the sign we have to resort to theoretical input such as lattice calculations or QCD sum rules. 
We then find that the (re-phasing invariant) ratios fsF^^^^ I fTt+ ^ JbAq^'' / fp+ are positive (see, for 
instance, [60]), hence fsF^^'^ /fn- and the ratio of form factors and decay constants in PT^-p+ is 
negative. Ho-wever, in this paper we have ignored the possible phase and sign conventions for the decay 
constants and form factors, always assuming them positive. Hence, we must write 

o - + = -e-'^''' ^ 

in the naive factorization approximation. In QCD factorization these considerations carry over to the 
full amplitude. This explains the minus sign in the B-meson decay amplitude of (105). 
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Table 13: Parameters of the time-dependent B 
tries as defined in the text. S and are computed for (3 
spending to sin(2/5) = 0.734. 



TT^p^ decay rate asymme- 



23.6°, corre- 



Theory 



SI 



S2 



S3 



S4 



Experiment 



^CP 


0.01 


C 


0.00 


s 


0.13 


AC 


0.16 


AS 


-0.02 



+0.00+0.01- 
-0.00-0.01- 
+0.00+0.01- 
-0.00-0.01- 
+0.60+0.04- 
-0.65-0.03- 
+0.06+0.23- 
-0.07-0.26- 
+0.01+0.00- 
-0.00-0.01- 



-0.00+0.10 

0.00-0.10 
-0.00+0.02 

0.00-0.02 
-0.02+0.02 

0.01-0.01 
-0.01+0.01 

0.02-0.02 
-0.00+0.01 
■0.00-0.01 



0.01 
0.01 
1.00 
0.03 
0.00 



0.02 
0.00 
0.20 
0.01 
-0.02 



-0.08 
0.02 
0.14 
0.15 
0.00 



-0.08 
0.05 
0.09 
0.00 

-0.03 



-0.21 ±0.08 
0.36 ±0.18 
0.19 ±0.24 
0.28 ±0.19 
0.15 ±0.25 



Defining P,rp/^7rp — ae^^", Pfm/Tf^ — —be^^^ (such that 6a, Si, are small and a, b are 
positive), and ApT^Tfy,^/ (At^pTt^p) — Re^^"^ , and treating a, 6, and 5t as small, the leading 
terms read 

2 

^cp = , , „o (asin(5a + R^bsm.5,,) sin7 + . . . , 
4i?2 

C = ^ ^^^^ (asin^a - 6sin(5b)sin7 + . . . , 

^ 1 + j^2 + (l + i^2-)2 (QCOS^a + &COsgfc)cOS7 + ... , 

^ 2R 2R ( ^ /2sin2a , ^ A 

- 6 cos bb I -^j^^2 ^ + sin(2/:; + 7) 1 ^ + ■ ■ ■ , 

2i? ^ . , 2R \ . ^ /2sin2a . 
l + i^2 cos2Q;sin6T- _^ ^2 1 " ^« I _^ ^2 siJ^7 + cos(2/j + 7) 

^ . ^ f2R^sin2a . .^^ A 1 

+ 6sin(5bl ^^^2 sin7 + cos(2/? + 7) 1 > + ..., 

with a = TT — /? — 7. The numerical values of a and 6 are given in Table 12. With our 
input parameters the tree-to-tree ratio is given by = 0.91 21 with phase 5t = (1±3)°, 
where the sizable error on R is entirely due to the form factors. 

The asymmetries Aq^, C, and AS* are suppressed by the pengTiin-to-tree ratios and 
the sin of a strong phase, hence they are always small in QCD factorization. This can be 
seen explicitly in Table 13. Note that A^ has a potentially large coefficient in front of 
sin 5t, however a large relative phase of the PV and VP tree amplitudes would constitute 
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Figure 12: Dependence of AC on the ratio of i?-meson form factors (for fixed 
F^-^(O) = 0.28) for 7 = 70° (center solid), 7 = 40° (upper dashed) and 
7 = 100° (lower dashed). The experimental value with its la error is also 
shown. Our central input value for the form-factor ratio is 1.32. 



a rather unexpected failure of QCD factorization. The data show a large value of C, 
which is related to the large direct CP asymmetry in B ^ 7r~p+, since 

Acp (1 -AC)-C , . + AcP (1 + AC) + C 

It seems impossible to accommodate a penguin-tree interference as large as needed to 
reproduce the central experimental values of Acp and C within QCD factorization. 

In first approximation AC is determined by the ratio R alone, with the largest uncer- 
tainty due to the form factor ratio Aq^'' (0) / F^^'" (0) . Since R is close to 1, the second 
term in the expansion in the penguin-to-tree ratios is also important and introduces some 
dependence on 7. This ratio is graphically displayed in Figure 12. 

The most interesting asymmetry is S. The numerical values of the parameters R, a, b, 
6a, and 5b are such that the first correction to S in (106) nearly cancels, hence S is nearly 
proportional to sin(2a). Furthermore, the expression 2R/ (1 + i?^) is not very sensitive to 
R near R = 1, leaving little uncertainty from the 5-meson form factors. This suggests 
that a or, more precisely, 7 (given the B-B mixing phase) can be accurately determined 
from this asymmetry. In the upper panel of Figure 13 we overlay the experimental la 
band for 5* to the theoretical prediction given as a function of 7 for (3 = 23.6°, which 
(up to a discrete ambiguity) corresponds to the current measurement of the B B mixing 
phase. The figure also contains two further curves, which define the theory error band, 
including the error ±2.4° on f3. Since this observable might be used to put a constraint 
on 7, we adopt an estimate of the theoretical error where the parameters qa that model 
the size of weak annihilation are allowed to vary independently for the PV and VP 
annihilation coefficients. This is more conservative than the standard error estimate, 
since it allows for a significant difference in the magnitudes of Pttp/T^^p and PpT,/TpT,, 
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Figure 13: Dependence of S on the angle 7 for j3 = (23.6 ±2.4)°. Upper panel: 
S from the vr^p^ final states. Lower panel: St^t^ from the 7r"7r+ final state. 
The experimental values of 5* and St^t^ with their la errors are overlaid. The 
dashed curves specify the theoretical error. 

which breaks the near cancellation in the correction term to S in (106). The figure 
shows that the theoretical error remains small even with this conservative estimate. 
Using Sexp = 0.19 ± 0.24 we obtain the allowed ranges 

7 = (72 ±11)° or 7 = (151 ±10)° (from^), (108) 

where the limiting values are defined by the intersection of the theory error band with 
the la experimental error band. The error given is dominated by the experimental error. 
The first range is in reasonable agreement with the standard unitarity-triangle fit. Note 
that if we allow for the second solution for /?, given by /? — 7r/2 — /5, we obtain the above 
ranges with 7 180° — 7 to a very good approximation. 

It is interesting to compare the constraint on 7 from S with the constraint obtained 
from the analogous asymmetry in the decay to the tt^tt^ final state. Here the correction 
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due to the penguin amplitude is larger. The current average experimental value St^t^ = 
-0.49 ± 0.27 yields 

7=(66^}^)° or 7 = (174^9)° (from 5.,) , (109) 

as shown in the lower panel of Figure 13. Despite the rather different values of the 
asymmetry, the allowed values of 7 are consistent with each other, in particular for the 
first range that is also compatible with the standard unitarity-triangle fit. This provides a 
strong argument for the validity of the theoretical framework underlying the computation 
of the hadronic amplitudes. While this is reassuring, one should keep in mind that there 
are several aspects of the current nn and np data that cannot be accommodated in QCD 
factorization, in particular the large direct CP asymmetries for Tr^Tr"*" and t:~ p'^ . Other 
pieces of data, such as the measurements of S'tttt, are not compatible between the BaBar 
and Belle experiments. It remains to be seen whether these anomalies disappear without 
affecting the apparent consistency of the remaining data. 

The analysis of the up final states has been performed under the assumption that the 
p-meson width is negligible. The experimental error from the quasi two-body assumption 
on the extraction of the parameters of the time-evolution of — > 'K^p'^ decays has 
been estimated as 0.08 for the CP-violating parameters ^cP) C*, and S", and as 0.03 for 
AC and AS* [61]. This leaves enough room for an improvement of the experimental 
error, so that we can look forward to more comprehensive unitarity-triangle constraints 
from the up system once higher statistics data samples become available. 

7 Final states containing r/ or r)' 

In this section we give results for the 23 B~ and decays into final states containing 
the mesons rj or 77'. A dedicated investigation of the AS* = 1 decays in QCD factorization 
has already been performed in [18]. For comparison of the different scenarios that we 
defined in the present paper we repeat the results for these decays in Tables 14 and 15. 
Small differences in the "Theory" column relative to the "Default" in [18] result from 
slight changes in the input parameters as mentioned in Section 3. 

We shall not repeat here the discussion of the complex dynamics of decays to mesons 
with flavor-singlet components that is particularly important for the penguin-dominated 
A^" = 1 decays. Let us briefly summarize the conclusions reached in [18]. The singlet 
penguin amplitude is presumably small and plays no important role in the enhancement 
of the B r]'K branching fractions. Rather, the particular pattern of the B — > rf^'^K^*^ 
decay rates is caused by the interference of different non-singlet penguin amplitudes. 
In particular, the large B — > vj'K branching fractions are obtained naturally in QCD 
factorization. There are, however, large uncertainties in applying QCD factorization to 
r^^') flnal states related to 77-77' mixing, a possible singlet contribution to the B — > 77 ^'^ 
form factors, and a novel, soft spectator-scattering term. The numerical results given in 
the tables and their errors corroborate these findings. 

The AD = 1 decay modes all have small branching fractions (and, perhaps, large CP 
asymmetries) unless they involve the color- favored tree coefficient a\. The corresponding 



61 



Table 14: CP-avcraged branching ratios (in units of 10^^) of penguin- 
dominated B — > PP decays (top) and B — > PV decays (bottom) with = 1. 
The errors and scenarios have the same meaning as explained in Section 5. 



Mode 


Theory 


SI S2 S3 S4 


Experiment 


B- rjK- 
B- ri'K- 


-1 Q +0.5 +2.4 +0.5 +1.6 

-■-■^-0.5-1.6-0.6-0.7 
1 -1 +0.1+2.0+0.4+1.3 
-■-■-■--0.1-1.3-0.5-0.5 

AQ 1 +5.1+26.5+13.6+33.6 
^y--"- -4.9-16.3-7.4-14.6 
AC cr +4.7+24.9+12.3+31.0 
"lU-O -4.4-15. 4 -6.8-13. 5 


1.1 1.6 2.4 1.6 
1.0 1.1 1.6 1.1 
52.1 78.3 64.6 76.1 
46.0 72.8 60.7 70.3 


3.1 ±0.7 

< 4.6 
77.6 ± 4.6 
60.6 ±7.0 


B- r]K*- 
B° riK*° 
B- ri'K*- 
B° rj'K*^ 


in 0+1. 9+8. 1+1. 8+16. 5 
-"-"-■■o -1.7 -4.4 -1.3 -5.5 
-1 A 7+1.1+7.8+1.4+16.2 
-'-'-'■' -1.0-4.3-1.2-5.5 
r 1 +0.9+7.5+2.1+6.7 
^■-■- -1.0-3.8-3.0-3.3 
Q+0.4 +6.6+1. 8+6.2 
'^■^ -0.4 -3.3 -2.5 -2.9 


14.0 19.4 19.1 19.9 
10.7 18.2 18.6 18.6 
3.5 7.6 7.1 2.2 
3.7 6.7 6.0 1.9 


25.4 ±5.3 
16.4 ±3.0 

< 35 

< 13 



results are summarized in Tables 16 and 17. We discuss specifically here only the AD — 1 
decays to final states with one pion or p meson. The two independent decay amplitudes, 
simplified according to Section 4.2, are given by 



-2 A 



A. 



-■nq 



+ A 



riq-K 



A. 



+ A 



^pu 



Spu ("2 - Pi) + 2al + a 



+ V2A^^^al + V2A^,,J, 



0(2 



-'pu 



-OL2 - Pi) + "4 



(110) 



where = a^+Pg^. The amphtudes for B — > nrj' and B — > prj'^'^ are obtained from these 
results by replacing (tt, 77) (tt, V) and (tt, 77) {p,V^'^)^ respectively. We find that the 
singlet coefficient and the "charm content" in the 77*^'^ are small. The amplitudes can 
then be further approximated by setting and A^^ ^(/) to zero. 

The neutral 5-meson decays to 7r°77*-''* and p'^77'^') have small branching fractions, be- 
cause the color-suppressed tree amplitudes proportional to 02 tend to cancel each other. 
A consequence of this is that the B° 
smaller than the B^ 

tudes are smaller. The charged decays have branching fractions of a few times 10^^. We 
note from (110) that the penguin coefficient enters both parts of the amplitude with 
equal sign. We should therefore expect larger penguin-tree interference in B~ — > 7r~77 
than in B^ — > tt^tt"^. Data are available for B~ 7r~77, whose branching fraction is in 



p°77^'^ decay rates are predicted to be much 

Ti'^rj^'^ rates, because the residual PV and VP penguin ampli- 
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Table 15: Direct CP asymmetries (in units of 10 ^) of penguin-dominated 
B — > PP decays (top) and B — > PV decays (bottom) with AS — 1. 



Mode 


Theory 


SI S2 S3 S4 


Experiment 


B" r]K~ 

B- r]'K- 


inn +6 4+11 7 +4 8 +25 3 
— ±0.y _6.9_i7.5 -8.5-21.8 
Q n+2.8+ 5.4+2.8+8.2 
y-*-* -2.1 -12.6 -6.2 -7.8 
9 /I +0.6+0.6+0.3+3.4 
^■^-0.7-0.8-0.4-3.5 
-1 Q +0.4+0.3+0.1+0.8 
-•-■^ -0.5 -0.3 -0.2 -0.8 


-33.3 -19.6 2.7 9.6 
-10.2 -10.3 -3.5 0.5 
2.3 1.6 -0.6 -0.8 
1.9 1.4 1.2 0.7 


-32 ± 20 

2±4 
8± 18 


B- rjK*- 
B- 7]'K*- 


q cr +0.9+1.9+0.8+20.7 
'J-'-' -0.9 -2.7 -0.8 -20.5 
Q Q +0.9 +1.1 +0.2 +3.8 
"^-O-l. 1-0.8-0.2-3.5 
-lA 9+4.7+8.5+4.9+27.5 
i + .Z _4.2 -13.8 -14.6 -26.1 
c cr +1.6+3.1+1.8+6.2 
^■^ -1.3-5.1-5.9-7.0 


2.7 3.0 -9.2 -5.7 

3.8 2.9 1.5 0.8 
-20.8 -9.8 4.5 22.1 
-5.8 -4.1 -2.4 1.7 


-5 ±28 
17 ±27 



Table 16: CP-averaged branching ratios (in units of 10 ^) of tree-dominated 
B — > PP decays (top) and B — > PV decays (bottom) with AD — 1. 



Mode 



TlioorA" 



SI 



S2 



S3 



Si 



)oriiiieiit 



B- 

B- 
5° 



7r°77 
► 7r~r]' 



A 7+1.9+1.8+0.6+0.4 
-1.7-1.5-0.3-0.3 



0.28 



+0.09+0.43 +0.02 +0.19 
-0.08-0.26-0.02-0.08 



q 1 +1.3+1.2+0.6+0.3 
"J--"- -1.2-1.0-0.3-0.3 
A -17+0.05+0.28+0.10+0.14 
-0.05-0.16-0.05-0.06 



2.7 

0.45 
1.8 
0.27 



4.1 
0.31 
3.1 
0.36 



4.7 
0.35 
3.1 
0.24 



3.8 
0.30 
2.9 
0.35 



3.9 ±0.9 

< 2.9 
< 7 

< 5.7 



5° 



7777 
7777' 



n 1 « +0.03 +0.43 +0.09 +0.10 

^■^'^ -0.03 -0.18 -0.03 -0.05 
A 1^+0.04+0.59+0.14+0.07 
^■^^ -0.04 -0.16 -0.06 -0.05 

u.uu _o.()i -0.05 -0.0:-! -0.0:-! 



0.13 
0.10 

0.03 



0.24 
0.30 

0.13 



0.23 
0.20 

0.10 



0.27 
0.31 

0.16 



B- 

B- 
5° 



► 77P 
77P° 

► 77'p~ 
77V° 



A 4+4.6+3.6+0.7+0.7 
^■^-3.7-3.0-0.4-0.7 
A AO +0.02 +0.16 +0.02 +0.05 
U.UO _o.oi -0.10 -0.01 -0.02 

^ q +3.1 +2.4 +0.5 +0.5 
"■•^-2.5-2.0-0.3-0.5 

A A-i +0.01+0.11+0.02+0.03 
"-'■'-'-'- -0.00-0.06-0.00-0.01 



8.3 
0.05 
6.0 
0.01 



6.3 
0.04 
4.3 
0.03 



9.1 
0.06 
6.1 
0.03 



6.3 
0.09 
4.2 
0.06 



< 6.2 

< 5.5 

< 33 

< 12 



B^ 
B' 
B^ 



r]uj 
rj'cu 

ri(j) 
r]'(j) 



A qi +0.14+0.16+0.35+0.22 
^■"^^ -0.12-0.11-0.14-0.16 
A 9A +0.10+0.15+0.25+0.15 
^■'^^ -0.08 -0.05 -0.10 -0.11 



0.001 
0.001 



0.23 
0.17 
0.001 
0.002 



0.53 
0.36 
0.001 
0.002 



0.44 
0.29 
0.000 
0.002 



0.65 
0.44 
0.001 
0.003 
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Tabic 17: Direct CP asymmetries (in units of 10~^) of tree-dominated B PP 
decays (top) and B PV decays (bottom) with AD = 1. We only consider 
modes with branching fractions larger than 10^''. 



Mode 


Theory 


SI S2 S3 S4 


Experiment 


B" — > 7r~?7 

B- TT-T)' 


14 Q +4.9 +8.3 +1.3 +17.4 
■-^—5 4— 74—08—173 
1 7 Q+5.2+ 7.9+1.2+33.4 
-■-'•y -4.1 -14.1 -1.4 -32.9 
o +2.8+10.5+0.7+20.4 
"•"-3.1- 9.0-0.7-20.4 
-iQ 0+5.5+7.7+4.1+35.7 
-•-^•^ -4.3-7.8-3.3-35.8 


-25.9 -15.4 0.7 5.6 
-11.0 -15.8 3.8 8.5 
-14.9 -10.5 8.6 11.1 
-11.9 -12.1 1.9 6.1 


-51± 19 

— 


B^ — > rjT] 
B^ r]r]' 


ao o +18.1 +65.6 +19.0 +23.4 
-13.1 -22.2- 9.5-16.2 
rr- O+17.1+141.0+20.7+13.7 
OD.O_iQQ_ 16.1-15.8-16.3 

4fi n +14-6 +138.4 +22.0 +40.2 
'±U.U_^4g_ 11.0-17.5-34.0 


-79.1 -41.3 -45.1 -20.1 
-87.2 -34.3 -46.8 -28.0 
-75.3 -27.8 -18.4 -3.5 




B" — > rip~ 
B- ^ r)'p- 


9 4+0.7+6.3+0.4+0.2 
^■^-0.7-6.3-0.4-0.2 
/I -1 +1.2+7.9+0.5+7.0 
^■^ -1.1-6.9-0.8-7.0 


-2.7 -1.9 -2.4 -6.0 
4.3 5.2 4.1 1.1 




S° r)u 
B^ rj'u 


oo A +10.0+65.3+20.9+19.2 
9.5-55.8-21.4-20.8 
A 9+0.1+53.0+11.6+20.4 
'-'■^ -0.1 -76.5 -11.5 -20.1 


-44.9 -10.5 -31.4 -27.5 
0.2 3.7 -1.9 -15.2 





good agreement with our result. The upper limit on the r]p~ final state is already below 
or near our prediction, and hence we expect this decay to be discovered soon. 

It should be noted that some of the decay rates for final states containing an rj' meson, 
in particular the modes B~ — * n'r]', p~f]', are very sensitive to the singlet contribution 
to the B ^ rj' form factor in (62), which we simply put to zero. For instance, setting 
F2 = 0.1 increases Bt{B- n-rj') to about 7 x 10"^ and Bt{B- p'ri) to 16 x 10"^. 

The direct CP asymmetries in Table 17 confirm the possibility of significant penguin- 
tree interference in the decays B^ n^r]^'\ However, the measurement of Acp(7r~?7) = 
(—51 ±19)%, together with a large CP asymmetry Acp{71~tt~^) = (51 ±23)% of opposite 
sign and a large asymmetry Acp{riK~) = (—32 ± 20)% of the same sign, is difficult to 
understand in QCD factorization and presumably in any framework (see, for instance, 
the comments in [62]), unless one can accommodate very large SU(3) flavor-symmetry 
breaking in the amplitudes pertaining to the vrvr, 7rr]^'\ and Kr]^'^ final states. Indeed, if 
the dominant interference is between and ai, then the two CP asymmetries must have 
the same sign, barring very large final-state dependence of these coefficients. Possible 
corrections could come from the singlet penguin amplitude ctg, but for this amplitude 
to have an effect as significant as indicated by the data, the estimates in [18] or [62] 
would have to be grossly in error. This leaves 0:2, which if large could affect the phase of 
the penguin-to-tree ratio by a noticeable amount. However, for this to reverse the sign 
of a large CP asymmetry the color-suppressed tree coefficient would have to exceed the 
color-allowed one (cti), which would affect all color-suppressed B decays. 
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8 Results for decays 

8.1 Simplified expressions for the decay amplitudes 

So far we have focused our attention on the decays of B~ and mesons, which are 
currently under investigation at the B factories. Our formahsm apphes equally to the 
decays of Bs mesons, which cannot be studied at these facilities. These decay modes will 
soon become accessible at hadronic B factories operated at the Tevatron Run-11, and in 
the longer term at LHC-b and BTcV. In this section wc present results for the branching 
ratios and CP asymmetries of Bg decays into PP and PV final states. Our predictions 
are collected in Tables 18-24. 



8.1.1 Decays with AS = 1 

Tables 18 and 19 contain results for penguin-dominated Bg decays with AS — 1. The 
modes with two pseudoscalars in the final state have large branching fractions, which in 
fact are among the largest of all rare B decays studied in this work. The expressions 
for the Bg — > KK amplitudes, simplified according to the approximations described in 
Section 4.2, are given by 

(111) 

(s) 

The expressions on the right-hand side must be multiplied with Xp and summed over p = 
u, c. The amplitudes for Bg — > KK* and Bg — > K*K are obtained from these expressions 
by interchanging K <-> K* or X <-> K* everywhere. Apart from the small annihilation 
terms parameterized by /^l these modes are characterized by a simple pattern of tree- 
penguin interference, which resembles that in the decays B~ — > 7r~K^ and B^ 7i^K~ 
in (73). This means that many of the analysis strategies discussed in Section 5.1.1 
can be applied also here, once experimental data on the decays Bg ^ KK (and the 
corresponding PV modes) will become available. The expressions for modes involving 77 
or T)' are more complicated and will not be given here. The exact decay amplitudes for 
these modes can be found in Appendix A. Note from Table 18 that the branching fractions 
for the decays Bg — > ri^'^uj and Bg ri^'^(j) are several orders of magnitude smaller than 
the corresponding branching fractions for the decays Bg — > r]^''^r]^'\ For the rj^'^u final 
states this occurs because the cu meson is assumed to have no strange component. For 
the case of the r]^'^(f) final states there is a strong cancellation between the PV and VP 
penguin amplitudes a2(^*-'V) ^^(^ al{(f)r]^'^) . The corresponding branching fractions can 
be enhanced by an order of magnitude by choosing a different value for rUg or giving up 
the assumption of universal annihilation. The direct CP asymmetries of the penguin- 
dominated Bg decays with AS — 1 are predicted to be small except for the modes with 
very small branching fractions. 
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Table 18: CP-averaged branching ratios (in units of 10^®) of penguin- 
dominated Bs PP decays (top) and Bs PV decays (bottom) with 
A^* = 1. The errors and scenarios have the same meaning as in Section 5. 



Mode 
K+K' 



B, ^ ^o^o 



Theory 

99 7 -1-3.5 -1-12.7 -1-2.0 -1-24.1 
-3.2- 8.4-2.0- 9.1 
r,A 7 +2.5 4-13.7 -)-2.6 +25.6 
-2.4- 9.2-2.9- 9.8 



1 ^ fi-hl-6+9.9+2.2+13.5 
-•-^■"-1.5-6.8-2.5- 5.5 
rA A +5.5 +32.4 +8.3 +40.5 
-5.2 -22.4 -6.4 -16.7 



SI 
28.0 
24.0 



S2 
33.4 
35.6 



S3 
34.3 
36.7 



S4 

36.1 
38.3 



Bs 



rjrj 

vv' 



11 



.- + 1.2^26,3^15,2^36.6 
' —1.0-17.2- 8.5 -15.-1 



14.9 
52.4 

11.1 



22.7 
79.5 

63.1 



21.4 
72.3 

60.2 



22.5 
77.7 

65.5 



Bs 
Bs 
Bs 

B.. 



K+K*- 

K°K*° 

K-K*+ 



A 1 +1.7+1.5+1.0+9.2 
^■^ -1.5-1.3-0.9-2.3 

Q +0.4 +1.5 +1.3 +10.4 
"J-^ -0.4-1.4-1.4- 2.8 
r r +1.3+5.0+0.8+14.2 
^■^-1.4-2.6-0.7- 3.6 

1 .9 +0.1+1.6+1.1 ^13.2 

-0.1-2.2-0.9- 3.2 



7.4 
3.8 
3.3 

1.1 



4.3 
4.2 
9.9 

8.7 



9.0 
9.1 
13.0 

11.1 



13.7 
14.3 
9.0 

7.9 



Bs 
Bs 
Bs 
B. 



r]ijj 
r]'uj 

77V 



n ni 9 +0-005 +0-010 +0-028 +0.025 
U.UIZ _o.oo4 -0.003 -0.006 -0.006 
n nO/l H-O.Oll +0.028 +0.077 +0.042 
-0.009 -0.006 -0.010 -0.015 



+0.02- 
-0.02- 



-0.95+0.54 +0.32 
0.14-0.12-0.13 



0.12 



n nc; +0-01 -t-l.lO +0.18 +0.40 
U.UO _Q 01 -0.17-0.08 -0.04 



0.017 
0.024 
0.15 
0.05 



0.011 
0.024 
1.02 
1.08 



0.010 
0.039 
0.24 
0.07 



0.009 
0.033 
1.47 
2.10 



Table 19: Direct CP asymmetries (in units of 10~^) of penguin- dominated 
Bs PP decays (top) and Bs PV decays (bottom) with — 1. We only 
consider modes with branching fractions larger than 10~^. 



Bs 
B. 



Mode 

K+K- 



Theory 



4.0: 



-1.0+2.0+0.5+10.4 
-1.0-2.3-0.5-11.3 

A A +0.2 +0.2 +0.1 +0.2 
'-'■^-0.2-0.2-0.1-0.3 



SI 
3.2 
0.9 



S2 

3.0 
0.7 



S3 



S4 



-4.5 -4.7 
0.6 0.6 



Bs 
Bs 
B. 



7]7] 

rjr)' 



_i « -(-0.5 +0.6 +0.4 +2.2 

-■-■"-0.4-0.6-0.7-2.2 
A 4+0.1+0.3+0.1+0.4 
"■^-0.1-0.3-0.1-0.3 
9 I +0.5+0.4+0.2+1.1 
^■^ -0.6-0.4-0.3-1.2 



-1.7 -1.1 -0.4 -0.1 

0.4 0.4 0.3 0.2 
2.1 1.7 1.3 1.0 



Bs 
Bs 
Bs 
Bs 



K+K*- 
K-K*+ 



9 9+0.6+8.4+5.1+68.6 

^■^-0.7-8.0-5.9-71.0 
-1 7+0.4+0.6+0.5+1.4 
-■-■ ' -0.5-0.5-0.4-0.8 
_o 1 +1.0+3.8+1.6+47.5 
'^■^ -1.1-2.6-1.3-45.0 
A 9+0.0+0.2+0.1+0.2 
"■^-0.1-0.3-0.1-0.1 



1.2 
1.7 
-5.2 
0.2 



1.8 
1.5 
-2.4 
0.3 



-34.8 

1.0 
18.8 

0.1 



-10.0 

0.8 
26.6 

0.1 



Bs 
Bs 



r](f) 
V'<f> 



-8.411?: 



-30.1+14.6+36.3 
-71.2-44.7-59.7 
(?r, 9+15.9+132.3+80.8+122.4 
""^■^-10.2- 84.2-46.8- 49.9 



-6.5 
-61.1 



3.8 
-8.9 



-8.7 
-34.0 



-9.5 
7.5 
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Tabic 20: CP-averaged branching ratios (in units of 10 ^) of tree-dominated 
Bg — > PP decays (top) and Bg — > PV decays (bottom) with = 1. 



Mode 


Theory 


SI S2 S3 S4 


Bg ^ 7r0^7 

Bg TT^T]' 


n nyc; +0-013 +0.030 +0.008 +0.010 
V.VIO _o.012 -0.025 -0.010 -0.007 
n 1 1 +0.02+0.04+0.01+0.01 
-0.02 -0.04 -0.01 -0.01 


0.097 0.069 0.077 0.073 
0.15 0.10 0.10 0.10 


Bg 71% 

Bg ^ ttV 


a; 0.0005 

r, 1 9+0.03+0.04+0.01+0.02 
^ -0.02 -0.04 -0.01 -0.01 


0.0004 0.0004 0.0033 0.0024 
0.16 0.12 0.12 0.12 


Bg ^ p^T] 
Bg ^ p^Tj' 


A 1 7+0.03+0.07+0.02+0.02 
^•J^' -0.03-0.06-0.02-0.01 
n 9 c +0.06+0.10+0.02+0.02 
•J-^^ -0.05 -0.08 -0.02 -0.02 


0.23 0.17 0.18 0.17 
0.35 0.24 0.24 0.24 



Table 21: Direct CP asymmetries (in units of 10"^) of tree-dominated Bg 
PP decays (top) and Bg PV decays (bottom) with A^* = 1. We only 
consider modes with branching fractions larger than 10~^. 



Mode 


Theory 


SI S2 S3 S4 


Bg ^ TfO^' 


97 o +6.0 +9.6 +2.0 +24.7 
^ < .o _5 7 _2 Q -27.2 


19.5 26.1 36.6 35.3 


Bg ^ ttV 


97 9+6.1+9.8+2.7+32.0 
^'•^-6.8-5.6-2.4-37.1 


20.0 25.2 27.2 24.8 


Bg ^ p°^7 
Bg ^ pOr/' 


97 0+6.4+9.1+2. 6+25.9 
^ ' -0-6.7-5.7-2.2-28.4 
OS Q+6. 1+10.3+1. 5+24.8 
^°-^-7.5- 6.3-1.8-27.5 


21.0 25.1 16.4 15.6 
20.0 26.1 36.7 32.2 



Tables 20 and 21 contain results for tree-dominated Bg decays with AS — 1. The 
corresponding branching fractions are very small, typically of order few times 10^^. 
Accordingly, the direct CP asymmetries can be large in all cases but will hardly be ob- 
servable in the near future. As an example we quote the simplified amplitude expressions 
^Bs^ttOco = and 



A, 



pu 



Ci2 



(112) 



The amplitudes for modes with t] or rj' arc given in Appendix A. 

Some of the Bg decays with A^ = 1 receive only annihilation contributions. As 
shown in Table 22 these modes have tiny branching fractions of order few times 10~^ to 
few times 10"^. The simplified expressions for the corresponding decay amplitudes are 



Bs 



5. 



pc 



(113) 



The amplitude for Bg n^p^ is obtained from the first expression by interchanging 
7T ^ p everywhere. The expressions for the Bg — > tttt amplitudes are obtained by setting 

p — > TT. 
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Table 22: CP-averaged branching ratios (in units of 10 ®) of annihilation-do- 
minated Bs — > PP decays (top) and Bg — > PV decays (bottom) with A^* = 1. 



Mode 


Theory 


SI S2 S3 S4 


Bs n^n^ 

Bs TT^TT^ 


n 09/1 -1-0.003+0.025+0.000+0.163 
^■^'^^ -0.003 -0.012 -0.000 -0.021 
n ni 9 +0-001 +0-013 +0.000 +0.082 
U.UIZ _o.ooi -0.006 -0.000 -0.011 


0.027 0.032 0.149 0.155 
0.014 0.016 0.075 0.078 


Bs ^ TT+p- 

Bs 7rV° 


^ 0.003 
^ 0.003 


0.002 0.003 0.019 0.014 
0.002 0.003 0.019 0.017 


Bs ^ 7T-p+ 


!^ 0.003 


0.002 0.003 0.019 0.015 



8.1.2 Decays with AD = 1 



Most Bs decays with AD = 1 are dominated by tree topologies. The branching fractions 
and direct CP asymmetries of these decays are given in Tables 23 and 24. The decays 
Bs — > 7r~K'^, 7r~ir*+, and p~K'^ have large branching fractions of order (l-2)-10~^. The 
corresponding neutral modes have much smaller rates. The direct CP asymmetries are 
predicted to be of moderate, sometimes even large magnitude, ranging from order 10% 
for the charged modes to significantly larger values for the neutral modes. The simplified 
expressions for the corresponding decay amplitudes are 



Spu "2 - "4 



(114) 



The right-hand sides of the expressions must be multiplied with X^^ and summed over 
p = u,c. The amplitudes for Bs — > nK* and Bs pK axe obtained by interchanging 
(tt, K) ^ (tt, K*) or (tt, K) ^ (p, K) everywhere. Note that these decays are governed 
by a relatively simple pattern of tree-penguin interference, which allows extractions of 
the penguin coefficients a.1 and gives sensitivity to 7. 

The branching fractions for the remaining Bg decays with AD = 1 are smaller, 
typically of order few times 10~^ (except for Bs — > K^r)') . The simplified expressions for 
the Bs —>■ K^u;, K^cj) decay amphtudes are 



V2A 



Kui 



(115) 



The amplitudes for modes with rj or rj' are given in Appendix A. In the first case the tree 
contribution is color suppressed, while the second process is a pure penguin decay. This 
explains the small branching fractions. Correspondingly, we predict generically large 
direct CP asymmetries for these modes. 
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Table 23: CP-averaged branching ratios (in units of 10 ^) of Bg 
(top) and Bs PV decays (bottom) with AD = 1. 



PP decays 



Mode 



Theory 



SI S2 S3 S4 



Bs 
Bs 



-rn 9+4.5+3.8+0.7+0.8 
-'-'-'•^-3.9-3.2-1.2-0.7 



0.49 



hO.28- 
-0.24- 



-0.22+0.40 +0.33 
■0.14-0.14-0.17 



6.8 

0.95 



8.1 
0.68 



10.4 
0.60 



8.3 
0.61 



Bs 
Bs 



n 0/1+0.19+0.64+0.21+0.16 
"-'•'-'^ -0.16 -0.27 -0.07 -0.08 
9 ri +0.3 +1.5 +0.6 +1.5 

^•'-'-0.3-1.1-0.3-0.6 



0.65 

2.6 



0.42 
3.0 



0.38 
2.7 



0.37 
2.9 



Bs 
Bs 



'+4.6+3.5 +0.7+0.8 



n 9C-+0.08- 
"-'•^'-'-0.08- 



-0.10- 
■0.06- 



-0.32 +0.30 
■0.14-0.14 



10.1 
0.15 



6.7 
0.39 



9.0 
0.36 



6.8 
0.33 



Bs 
Bs 



K+p- 



9/I cr +11.9+9.2+1.8+1.6 
^^•■J- 9.7-7.8-3.0-1.6 
A (-1 +0.33+0.21+1.06+0.56 
'-'•"-L -0.26-0.15-0.38-0.36 



21.3 
0.82 



19.2 
0.58 



24.6 
0.70 



191 
0.6J 



Bs 
Bs 



n p-1 +0.20+0.15+0.68+0.40 
'-'•'-'-'- -0.18-0.11-0.23-0.25 
n, 97+0.09+0.28+0.09+0.67 
'-'•^' -0.08-0.14-0.06-0.18 



0.33 
0.43 



0.50 
0.54 



0.58 
0.64 



0.63 
0.46 



Bs 
Bs 



r. 9(^+0.15+0.49+0.15+0.57 
^■'^^ -0.13 -0.22 -0.05 -0.15 
n 98+0.04+0.46+0.23+0.29 
"J-^" -0.04 -0.24 -0.10 -0.15 



0.52 

0.23 



0.29 
0.24 



0.54 
0.39 



0.57 
0.67 



Table 24: Direct CP asymmetries (in units of 10 ^) of Bs PP decays (top) 
and Bs — > PV decays (bottom) with AD = 1. 



Mode 



Theory 



SI 



S2 



S3 



S4 



Bs 
Bs 



77- K+ 



-6.7 



+2.1+3.1+0.2+15.5 
-2.2-2.9-0.4-15.2 



A-] c +16.6 +14.3+ 7.8+40.9 
^±.U _12.0 -13.3 -14.5 -51.0 



-10.0 
21.4 



-8.5 
33.6 



7.7 
14.3 



10.9 
4.6 



Bs 
Bs 



18.5+28.6+ 5.2+34.6 
13.2-32.2-12.5-45.6 



46.8 

_qc +8.6+6.0+3.8+19.3 
OU.U _g^2 -7.4 -2.5 -17.3 



24.1 
-28.8 



39.3 
-29.5 



28.6 
-22.2 



24.2 
-18.2 



Bs 
Bs 



7r-K*+ 



n ,-+0.2+1.4+0.1+19.9 
'-'•"-0.1-1.7-0.1-20.1 



/ic 7+14.3+13.0- 
-16.0-11.6- 



-28.4+80.0 
■28.0-59.7 



0.5 
-79.3 



1.6 
-41.9 



-18.5 
0.3 



_1 r +0.4+1.2+0.2+12.1 
i- O -0.4-1.4-0.3-12.1 
r>A 7+7.1+14.0+22.8+51.3 
-5.2-12.4-17.7-52.3 



-22.0 
15.4 



Bs 
Bs 



K+p- 



-1.7 
18.3 



-1.7 
24.5 



10.1 
-11.8 



6.2 
11.6 



Bs 
Bs 



AO Q +13.6+18.0+30.6+57.7 
^O.y _i3.4 _18.2 -30.2 -49.3 

in q +3.0+4.7+3.7+5.0 
lU.O_2.4_3.o-4.1-7.5 



-67.5 

-6.4 



-40.9 
-10.5 



-9.6 
-6.3 



An. 9+17.0+24.6+7.8+65.9 
^U.Z 5 _3o 8 _i4 -96.3 

_ro +16.9+41.4+19.9+44.9 
-11.9-11.7-13.9-35.7 



-30.1 
-7.4 



Bs 
Bs 



*o 



20.4 
-70.8 



36.3 
-54.5 



-11.7 
-24.9 



0.6 
-32.7 
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9 Conclusions 



In this paper we extended our previous analysis [10] of 5-meson decays to nn and 
ttK final states and to final states with r]^'^ mesons [18] within the QCD factorization 
approach to all two-body final states with pseudoscalar mesons or one pseudoscalar 
and one vector meson, including also decays of Bg mesons. The main motivation for 
performing a comprehensive analysis of all 96 final states is to obtain a global assessment 
of the phenomenology of QCD factorization with a consistent common input to all decay 
modes, and to display possible correlations between the various modes. 

The theoretical analysis follows [10]. The computation is performed at next-to- 
leading order in ctg for the hard-scattering kernels at leading order in l/m^, and simi- 
larly for the kernels that multiply the subleading twist-3 quark-antiquark distribution 
amplitudes. Spectator scattering effects at subleading power in l/nib (including weak 
annihilation), which do not generally factorize, are estimated by a phcnomenological 
model and assigned a 100% uncertainty (including an arbitrary strong- interaction phase). 
On the technical side, the generalization of the decay amplitudes computed in [10] to 
pseudoscalar-vector final states is for most parts straightforward, involving only a few 
sign changes in the decay amplitudes and a few new hard-scattering kernels. For the 
analysis of this paper only the kernels for the twist-3 quark-antiquark amplitudes of vec- 
tor mesons needed to be computed anew. We also discussed a new electroweak penguin 
effect that contributes only to neutral vector mesons. 

We have taken this comprehensive analysis as the occasion to summarize in a unified 
notation all results available at next-to-leading order in QCD factorization. The matrix 
elements of the effective weak Hamiltonian are decomposed according to their flavor 
structure. A complete list of all decay amplitudes expressed in terms of the flavor 
coefficients (which generalizes and simplifies the more familiar Oj and 6j notation) is 
given in Appendix A. The next-to-leading order results for the coefficients are given in 
Section 2, including the results from [10] for completeness. 

The comparison of theory with data shows many interesting effects, which we sum- 
marize here. We should note, however, that no experimental information exists to date 
for the majority of the 96 decay modes considered in this paper. This will allow further 
tests of the theory in the future. When passed successfully, this implies a rich source 
of information on fiavor-changing transitions in the quark sector from purely hadronic 
decays. 

Results related to the tztz, tzK final states 

A detailed discussion of these modes can be found in [10]. Since 2001 the experimental 
errors have been reduced by almost a factor of two. QCD factorization continues to 
provide a natural explanation for the magnitudes of the tree and penguin amplitudes 
relevant to these decays. The theory with default input parameters docs not fare well 
on the 7r~7r"'" decay mode, which it predicts too large, and on the n'^K^ mode, which 
it predicts too small. The former discrepancy is often interpreted as evidence of a 
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large value of 7 or large rescattering. We find that an alternative "hadronic physics" 
explanation is possible if the S — > tt form factor is about 15% smaller than usually 
assumed, the strange-quark mass is at the lower end of the currently favored range, and 
if the color- suppressed tree amplitude is enhanced, for example by a sizable spectator- 
scattering effect. While the evidence for this scenario is not conclusive now, crucial 
input will be provided by a measurement of the semileptonic decay rate near = 
0. Interestingly, this scenario is also favored by measurements of the decay rates for 

— > n'^p^. With regard to the "7r°^° anomaly", we find that the large experimental 
ir^K^ rate cannot be explained by a different value of 7, and that "hadronic 
physics" explanations appear extremely unlikely. While a significant modification of the 
electroweak penguin amplitude due to "New Physics" may explain the effect, we consider 
it to be more likely a statistical fluctuation. 

It will require more data to arrive at a conclusive picture for the direct CP asymme- 
tries. While it appears now certain that the asymmetries are small for the nK modes, in 
agreement with the predictions of QCD factorization, a quantitative comparison needs 
better statistics and a more accurate theoretical computation. The experimental situa- 
tion for the 7r~7r+ flnal state is still unsettled. If a direct CP asymmetry of order 50% 
is conflrmed in the future, the factorization framework would be in trouble. Similar 
comments apply to other observations of large direct CP asymmetries, which at present 
all have large experimental errors. 

Results related to the Trp final states 

The TTp system discussed in Section 6 exhibits some advantages for studies of CP viola- 
tion that render it highly interesting even within the limitations of the quasi two-body 
assumption. We deflned and discussed several ratios of irp branching fractions that 
should shed light on the magnitudes of the hadronic amplitudes underlying this class 
of decays. The theoretical predictions for these ratios are in good agreement with the 
available data. In particular, QCD factorization predicts that the two distinct penguin- 
to-tree ratios in the irp system are about a factor of three smaller than the corresponding 
ratio for tttt, and have smaller errors. If this can be conflrmed, it implies that 7 can be 
determined relatively accurately from time- dependent CP violation measurements. We 
considered the five quantities that parameterize the time-dependent asymmetries in the 
decays B^, B^ — » n^p^ in detail. If data and theory are taken at face value, wc deter- 
mine 7 ~ 70° with an error of about 10° from the asymmetry S. This value is consistent 
with a less accurate result obtained from the corresponding quantity in the tttt system (if 
one takes the average of the two experiments, which are mutually incompatible). This 
appears intriguing but should be taken with some reservation, since the central values of 
the direct CP asymmetries Aqp and C are once more larger than theoretically expected. 
We anticipate that, with more data available soon, the np system will play an important 
role in the understanding of hadronic decays and CP violation. 
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Results related to the penguin-dominated PV final states 

The comparison of the penguin- dominated final states nK, nK* and pK allows for 
crucial tests of the factorization framework, since the PP, PV, and VP penguin am- 
pHtudes take very different values for reasons specific to factorization. The nK* data 
arc indeed consistent with a smaller penguin amplitude; however, as has been noted 
before [20, 22, 23], the reduction is not as large as predicted by theory. With default 
parameters one underestimates the amplitude by about 40%, but the theoretical error 
on the penguin amplitude for vector-meson final states is large, in particular the one 
from weak annihilation. This taken into account, the prediction may be in agreement 
with the measured branching fractions within their errors. However, the situation is 
not satisfactory, since one would wish to have an explanation of the data that does not 
invoke weak annihilation. 

The estimate of weak annihilation is necessarily model-dependent in the QCD factor- 
ization approach. In particular, the error range (technically implemented by requiring 
qa < i) has to be specified as an "educated guess". It is therefore desirable to derive 
experimental constraints on weak annihilation. It was already found in [10] that the 
7iK branching fractions do not favor a sizable annihilation amplitude, since this would 
require a fine-tuning of its strong phase to keep the branching fraction small enough. 
On the other hand, such a coincidence cannot be excluded. Assuming this to happen 
and further assuming universal annihilation amplitudes, we find that several branching 
fractions of pseudoscalar-vcctor final states, in particular K(j), are much above the data 
for Qa > 2. This, together with an estimate of the pure annihilation mode D^K~ 
[63], seem to imply that weak annihilation cannot be much larger than the upper limit 
defined by the phenomenological treatment adopted in the present analysis and in [10]. 

As a by-product of the analysis of penguin-dominated PV modes, we also obtain an 
estimate of the difference of the time-dependent CP asymmetries in J /ipKs, <t>Ks, and 
rj'Ks decays. As may have been expected, the factorization approach does not contain 
any mechanism that could enhance the CKM-suppressed amplitudes with a different 
weak phase, limiting the CP-asymmetry differences to a few percent. 
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Note added in proof 



After submission of this paper the following experimental results relevant to our analysis 
have been published: (1) The branching fraction for the decay — tt^tt" now reads 
(1.9±0.5)-10~^ [98, 99], which is significantly larger than expectations. (2) The branching 
fraction for the decay K~ has been measured by Belle to be (3.9 ± 0.8) • 10~^ 

[100], in good agreement with the predicted suppression of the VP penguin amplitude 
discussed in Section 5.1. (3) The branching fraction for the decay t[~ K*^ from 

Belle now reads (8.5 ± 1.3) ■ 10^^ [100], less than half as large as before but in good 
agreement with the previous CLEO measurement. If this central value were confirmed 
by BaBar, this would remove one of the major discrepancies between data and the results 
of QCD factorization, see Section 5.1. 
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Appendix A: Explicit results for the decay amplitudes 



The results for the large set of Bq PP, PV decay amplitudes can be expressed most 
concisely in terms of traces over flavor matrices. We collect the three i?-meson states 
into a row vector B = {B^ , B^, Bs) and represent the final-state pseudoscalar and vector 
mesons by matrices 



V2^ V2^ V2 



n 



K- 



\ 



v 



V 



Vs + V'sJ 

K*- 



(116) 



respectively. (One should also add a cc component for rj and r]', whose contribution is 
given by (51).) In addition, we define a column vector 

/ \ 



A 



■p 

is) 



;ii7) 



VA-7 

containing CKM matrix elements, as well as matrices 

/ Sr,u \ / 1 \ 



-'pu 


\ 0/ 



Q = -^Q 



-| 

\o -y 



(118) 



Finally, we use the definitions of the quantities AMiM2i '^1(^1^2) and j3i{MiM2) given 
in (11), (9) and (18). Then the entity of all decay amplitudes Ab^MxM2 reproduced 
by evaluating the master expression 



Am,M2 jsMi (ai Up + al + < ew Q) ^2 A,, 



p=u,c 



+ BMiAp ■ Tr 



+ B (/?2 t/p + PI + /3f_Ew Q) MiM2Ap 
+ BAp ■ Tr Uj, + Pl + 6^^ew Q) ^1^2 

+ B (/352 U, + (3^s3 + /3|3,EW Q) M^A, ■ TVM2 

+ BAp • TV [(psi Up + + b's4,Ew Q) Ml] ■ TVMaj , (119) 



74 



where ctj = ai{MiM2) and Pi = Pi{MiM2). We recall that the ccj terms include ver- 
tex, penguin and hard spectator contributions, whereas the /Si terms result from weak 
annihilation. Contributions proportional to Up are from the current-current operators, 
those proportional to the unit matrix are from QCD penguins, and those proportional 
to the charge matrix are from electroweak penguins. 

To determine the amplitude for a particular decay, for instance S° — > Tr~^K~, one 
inserts P for Mi and M2 and extracts the terms corresponding to 7r'^K~. When 7r+ 
comes from Mi the prefactor is A^^k, otherwise it is Akw For pscudoscalar-vector final 
states both possibilities, Mi = P, M2 = V and Mi = V, M2 = P, must be summed 
and the Am^m-z prefactor is determined as above. For some pure annihilation amplitudes 
the form factor B — ^ Mi does not exist, and AmiM2 is not defined. In this case the 
expressions AM-^M2 Pi must be replaced by Bm^m^ h- See Section 2.2 for the definition of 
the relevant quantities. 

We now list our results for the various decay amplitudes expressed in terms of the 
and j3i parameters. Several sets of amplitudes have the same representation in terms of 
flavor parameters, apart from obvious substitutions of labels. 

A.l Decays with A5 = 1 

is) 

In this section the expressions for decay amplitudes must be multiplied with Ap and 
summed over p = u,c. Throughout, the order of the arguments of the af (M1M2) and 
/3f (M1M2) coefficients is determined by the order of the arguments of the Am-^M2 pref- 
actors. There is a total of 17 S — > PP and 31 B ^ PV amplitudes, which split up as 
(4,4,9) and (8,8,15) into the flavor states {B-,B^,Bs). 

B — > tiK'^*'> and B pK decay amplitudes 
There are four independent amplitudes, given by 

■^B-^'K-Ro — ^ttK [^pu /32 + 0:4 — |q;4 Ew + /^f + /^3,Ew] ' 

V2Ab-^7tOK- = Kk [^pu («1 + 132) + «4 + "4,EW + + /^3,Ew] 

+ ^^TT \Ppu Ot2 + |«3,Ew] ' 

AbO-^-k+K- — ■^■kK \Ppu Oil + + Q^^EW + ~ i/^3,Ew] ' 

Abo^t^oRo = A^R [—(^4 + |<^4,EW ~ /^3 + l/^3,Ew] 

+ [(5p„q;2 + |4ew] ■ (120) 
Isospin symmetry implies that 

The expressions for the B ttK* and B pK amplitudes are obtained by setting 
{ttK) — > {tiK*) and {tiK) — ^ {pK), respectively. 
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B K^*^ri^''> and B Kuj/(j) decay amplitudes 
There are two independent amplitudes, given by 

\/2AB--.K-r, = ^Kr^g [^pu {(^-2 + Wst) + 2q| + l^f^EW + 2/^53 + W^^^-Ew] 

+ ^ARtis Spu {P2 + PS2) +al + al- l^a^EW " i«4,EW + P3+ /^3,EW 
+ Ps3 + /^53,EW 

EWj ' 

\/2AB0^K^n = ^Krig [Spu ^2 + 2af + ^ag^EW + 2^93 - /^ss.Ew] 

+ -\/2A^^^ + — |q^^eW ~ i'^4,EW + /^f ~ i/^3,EW + /^53 ~ i/^53,EW 

+ [al - i<EW + /5f - i/^Cew] • (122) 

The amphtudes for B^ Krj', B^ K*r]^'\ and 5° Ku}/(f) are obtained from this 
result by replacing {Kr]) [Kt]'), {Krf) {K*ri^'^), and {Krj) — > {Kuj/cf)), respectively. 
When ideal mixing for uj and (/> is assumed, set Ap^,^^ and ^^^^ to zero. Furthermore, 
with our approximations A^^^ = Aj^^^ = 0. 

Bg —>■ TTTT and Bg — > vrp (ieca?/ amplitudes 
There are two independent amplitudes, given by 

+ ^pTT [5pu fel + &4 + fe^^Ew] ' 
2 ^Bs-»7rOpO = -^Trp [(^pu ^1 + 264 + |^_ew] 

+ S^[5p.&i + 26^ + |6^,Ew]- (123) 

The amplitudes for Bg 7r^p+ is obtained from the first expression by interchanging 
77 p everywhere. In the limit of isospin symmetry the following relation holds: 

2 As.-^^OpO = AB,-,n+p- + ^B^-»7r-p+ ■ (124) 

The expressions for the Bg — > tttt amphtudes are obtained by setting p — > tt. 



76 



Bg K^*)k^*) decay amplitudes 

There are two independent amplitudes, given by 

^Bs-^K-K+ = -S^K [Spu 61 + 64 + b^ Ew] 

+ A^a [5,^ a, + al + a^^^w + l^l + l^l - - ¥U^] ■ (125) 

The ampUtudes for i?s KK* and K*K are obtained from these expressions by 

replacing {KK) [KK*] and [KK] {K*K), respectively. 

Bg — > nr]^'\ Bs — > nuj/(f), and Bg — > pjy^'^ c?eca?/ amplitudes 
There is only one independent amplitude, given by 

2 AB.-^n^n ^ ^-^vq [^py. (^1 + 2^51) + i^^,EW + 3^5,Ew] 

+ ^A^^^ [5p^ a2 + Kew] ■ (126) 

The amplitudes for Bg ttt]', Bg 7ruj/(p, and Bg pi]^'^ arc obtained from these 
results by replacing (tt?]) — > {tii]'), {Tirj) {nu/cp), and {tit]) (p^*-'-*), respectively. 
When ideal mixing for cu and is assumed, set B^^i^^^ and Stt^^ to zero. 

Bg ■q^')-qi') and Bg — > 'r]^^ijjj<^ decay amplitudes 
There is only one independent amplitude, given by 

2 ■AB.^nn' = Brj^rj'^ [6pu (61 + 2651) + 26^ + ^If^^^w + + ^4,ew] 
+ ^Brj^v's [^PU ^Sl + 2i'54 + |&54,Ew] 

+ \/2Ajj^^>^ [6pu 0L2 + 2af + l^g^E^ + 2/3^3 + 2/3^4 - /^^g - (^^s^^^nA 

+ '^^r]sV's 0^ + 0(4^ — i'^3,EW ~ |'^4,EW + + /^4 ~ i/^3,EW ~ i/^4,EW 

+ /^S3 + f^Si ~ |/^53,EW ~ i/^54,EW 
+ '^^VsV'c [^PC + Oil] 
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+ B^'.v, [^vu {bi + 2bsi) +2lf, + ife^ EW + + ^4,Ew] 

+ V2Br,'^r,s [Spu bsi + 21)^34. + |^4,Ew] 

"I" "3 "T "4 2"3,EW 2"4,EW "T /-'3 "T /-'4 2'^3,EW 2'^4,EW 

+ /^SS + /^54 ~ |/^S3,EW ^ i/^54,EW 

+ 2A^,^J5pea2 + a|]. (127) 

The amplitudes for Bg r]r], Bg T]'r]', and Bg ri^'^uj/(f) are obtained from this 
result by replacing {rir]') (^7^7), (^^') (W)) ^'^^ ivv') ~^ respectively. 
When ideal mixing for uj and is assumed, set A (o , A (j) , B (/) , B (j) to zero. 

Furthermore, with our approximations A {,) = A (/> = 0. 



A.2 Decays with AD = 1 

In this subsection the expressions for decay amplitudes must be multiplied with Ap^^ 
and summed over p — u,c. Throughout, the order of the arguments of the af (M1M2) 
and /3f (M1M2) coefficients is determined by the order of the arguments of the Am-^M2 
prefactors. There is a total of 17 S ^ PP and 31 S PV amplitudes, which split up 
as (4, 9, 4) and (8, 15, 8) into the flavor states (S", Bg). 

B — >■ TTTT and B ^ np decay amplitudes 

There are three independent amplitudes, given by 

V2 AB--^n-pO — ^Trp [^pu {oi2 — P2) — OC^^ + |q;3,EW + |'^4,EW ~ ^3 ~ /^3,Ew] 

+ ^pTT [Spu {ai + (32)+al + <^Ew + ^3+ ^3,ew] > 



a 



4,EW 



p _ Irp _ 1 /qp 

2/^3,EW 2 



/^4,Ew] 



-2 A 



A 



np 



+ A 



PK 



^pu {<^2 — Pi) — Q!4 + 2'^3,EW + 2'^4,EW ~ ^3 ~ 2(5^ 



^ 2f-'3,EW 2f-'4:,EW 



Spu {0i2 — — a\-\- 2Qi3,EW + 2'^4,EW ~ ^3 ~ 2(5\ 



^ 2'^3,EW 2/^4,EW 



(128) 



The amplitudes for B and — > TT are obtained from the first two expres- 

sions by interchanging n ^ p everywhere. In the limit of isospin symmetry the following 
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relation holds: 

2ABO^T:OpO — AbO-^t^+p- + ^B0^7r-p+ ~ (^AB--.Tr-pO + ^S-^7rOp- j • (129) 

The expressions for the S — > tttt amplitudes are obtained by setting p — > tt. 

B K^*)k^*) decay amplitudes 

There are three independent amplitudes, given by 

•A.B-^K-K° = Ap^j^ \5pu P2 + " |Q^4,EW + + /^3,Ew] ' 
AbO-^K-K+ — ^KK [SpuPl + + /^4,Ew] 

+ [K - m,Ew] ■ (130) 

The amplitudes for B KK* and B — > are obtained from these by replacing 

{KK) {KK*) and {KK) {K*K), respectively. 

B — i> TTri'''\ B — > prj^'^ and B Tiuo/cj) decay amplitudes 
There are two independent amplitudes, given by 



\/2A. 



Spu {0:2 + P2 + Ws2) + Saf + + Ictg^EW ^ 2"4,ew 

+ /^f + /^3,EW + 2/353 + 2/3|3^EW 

+ ^/2A^^^ [5pu Ps2 + "3 - i"I,Ew + P'k + /5s3,Ew] 

+ V^^TTJJc \5pc OL2 + Q!3] 

+ AgTT [(^pw (ttl + /32) + + «4,EW + /^f + /^Cew] ' 

= A. 

•• -III 

+ /^3 ~ 2''^3,EW ~ i/^4,EW + 2/3^3 — Ps3,EW ~ ^(^S4,EW 
+ \/2At^jj^ [Spu Psi + <^3 ~ i<^3,EW + /^S3 ~ |/^S3,EW ~ i/^54,EW, 

+ V2A„ri^ [6pc a2 + af ] 
+ A^TT (-0:2 - Pi) + o^ 



1 p 



-2 A 



5pu (^2 - /5i - 2/351) + 2af + + ^tta^EW " i"4,EW 



(131) 



2'-^3,EW 2'^4,EW 



' 2'^3,EW 2'^4,EWj • 



The amplitudes for B tit]' , B pi]^'\ and B tilo / (j) are obtained from these results 
by replacing {nr]) — * (tt?]'), (tt?]) — > {pv^'^), and (vrr/) — > {nuj/cf)), respectively. When 
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ideal mixing for uj and is assumed, set A^^-i^^^ and A^^p^ to zero. Furthermore, with our 
approximations A^r^^^ — A^^^^ — 0. 

B — > Tj^'^Tj^'^ and B — > rj^'^uj/ct) decay amplitudes 
There is only one independent amplitude, given by 

2 Aso-^rm' = "^vW, ^pu («2 + + 2/?5i) + Saf + al+ |aif ew " i'^4,EW 

+ /3f + 2f3l — ^/5f^EW + |/^4,EW + 2/^53 + 4/5^4 — /^^s^eW + /^S'4,EW 
+ •\/274^^^/ \5pu Psi + Oi^— I Qis.EW + ^S3 + 2/354 — |/353^EW + i/^54,Ew] 

,EWj 

+ '^^VsV's ~ i^,EW + ^4 ~ i^4,Ew] 

+ Arj'^^^ 5pu {a2 + (3i + 2psi) + 2af + + |af,EW " ^"4,ew 

+ /3f + 2/3f — |/3f^EW + i/^4,EW + 2/353 + ^/^54 ~ /^53,EW + /^54,Ew] 
+ \/2Ajjf^jj^ [Spuf^Sl + 0^3 ~ |'^3,EW + I^S3 + 2/?54 — |/?53^eW + i/^54,Ew] 
+ 

+ V2Brj'^ri^ [2654 — &54^Ew] 

+ '^^ri'^Vs ~ |^,EW + ^4 ~ i^4,Ew] • (132) 

The amplitudes for B'^ ^ rjrj, B^ ^ rj'r]', and B^ — > rj'^'^Lo/cf) are obtained from this 
result by replacing {rjrj') — > (r^r^), (r???') — > {v'v'): ^^'^ ivv') ~^ {v^'^^/4>): respectively. 
When ideal mixing for u and is assumed, set A (,) , A (i) , B (,) , B (,) to zero. 

Furthermore, with our approximations A (,) — A (,) — 0. 

Bg —>■ ttK^*^ and Bg pK decay amplitudes 
There are two independent amplitudes, given by 

•^Bs^-K-K+ = Aktt \5pu + Q;4 + tt4^EW + /^3 ~ i/^3,Ew] ' 
•\/2 ^B^^TrO/fO — ^Kn [^pu a2 — + foif^EW + i'^4,EW ~ /^f + |/^3,Ew] " (133) 

The expressions for the Bg — > nK* and Bg pK amplitudes are obtained by setting 
(ttK) — > (ttK*) and (ttK) — > (pK), respectively. 
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Bg K^*\^'^ and Bg Kuj/cj) decay amplitudes 
There is only one independent amplitude, given by 



^/2A 



A 



+ V2A,^K H - Kew + /3f - i/3f,Ew] ■ 



(134) 



The amplitudes for Bg Ki]', Bg K*r]^'\ and Bg Kcu/cf) arc obtained from this 
result by replacing {Kr]) [Kt]'), {Kt]) {K*r]^'^), and {Kr]) (Ku/cf)), respectively. 
When ideal mixing for cu and (p is assumed, set Akuj^ -^k^^ to zero. Furthermore, 
with our approximations Akoj^ — Ak^^ — 0. 
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Appendix B: Convolution integrals 

The convolution integrals of the hard-scattering kernels with meson light-cone distri- 
bution amplitudes can be evaluated using expansions of the distribution amplitudes in 
terms of Gegenbauer polynomials. In [10] the corresponding expressions for final states 
containing two pseudoscalar mesons were given including the first three terms in this 
expansion. Here we list the corresponding results for the convolution integrals involving 
the twist-3 distribution amplitude which are needed for pseudoscalar-vector final 

states. We include the first three terms in the Gegenbauer expansion (27). 

The convolution integral entering the vertex corrections for the coefficients ag § is 

J dx ^y{x) h{x) = 9 - 6i7r - i7i^ + . . . , (135) 

with h{x) as given in (38). The hard spectator contributions involve the divergent 
integral 

[ dxp^^^^{l)XH-{Q + 9a^^ + na^^ + ...), (136) 
Jo ^-x 

where $„(1) = (^n,±- The penguin contributions involve the convolution 

Gvis) ^ [ dx ^y{x) G{s -ie,l-x), (137) 
Jo 

where G{s,x) is the penguin function defined in (40). We obtain 

Gv{s) = 1(6 + 2aX^ + al^) - 4s (9 + 12aX^ + 14a^^) 

- {SaX^^ + 35q;J^) + 360s^ 

+ 12sVl - 4s [l + (1 + 4s) aX^^ + (1 + 15s - 30s^) aj^] 
X (2 arctanhVl — 4s — in) 

- I2s^ [1 + (3 - 4s) aX^^ + 2(3 - 10s + ISs^) aj^] 

X (2 arctanhVl - 4s - ivr) ^ + . . . . (138) 
For the special cases s = and s = 1, this expression reduces to 

Gv(0) = 1 + ^ + ^ + . . . , 

Gv{l) = -35 + 4v^7r + ^ + (-^ + 20v^7r - ^) 

+ 56^.+ ^) + (139) 
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Appendix C: Summary of experimental results 



In Tables 25-28 we compile the available experimental data on the CP-averaged branch- 
ing fractions and CP asymmetries in B ^ PP and B PV decays, distinguishing the 
two classes of decays of the type = 1 and AD = 1. We also present our weighted 
averages of the data (ignoring correlated errors, which are small). Where measurements 
are inconsistent, the combined error is inflated by a factor of 5" = s/x^li^ ~ !)> which 
is shown in parenthesis. 



Table 25: CP-averaged branching ratios (top, in units of 10 ^) and CP asym- 
metries (bottom, in %) for B — > PP decays with AS — 1. Upper limits are at 
90% confidence level. We show S = VxV(^ - 1) 

in cases where 5" > 1. 



Mode 


BaBar Belle CLEO 


Average 


Penguin-dominated decays 


B- TT^K- 
& TT+K- 

B° ^ 7r"X" 


20.0 ± 1.6 ± 1.0 [65] 22.0 ± 1.9 ± 1.1 [84] 18.8;^:^ti;8 [91] 
12.8tJ:? ± 1-0 [73] 12.8 ± lAt\i [84] I2.<dtlit\\ [91] 
17.9 ±0.9 ±0.7 [73] 18.5 ± 1.0 ± 0.7 [84] 18.0t2:ito.9 [91] 
10.4 ± 1.5 ±0.8 [75] 12.6 ±2.4 ±1.4 [84] 12.8^^;°!^;^ [91] 


20.6 ± 1.3 
12.8 ± 1.1 
18.2 ±0.8 
11.2 ± 1.4 


B- rjK- 
fiO ^ rjK° 
B- rj'K- 
§0 ^ r]'K^ 


2.8 ± 0.2 [71] 5.3 t\l ± 0.6 [84] 2.2 (< 6.9) [96] 
2.6 1^]! ± 0.2 (< 4.6) [71] < 12 [84] < 9.3 [96] 
76.9 ± 3.5 ± 4.4 [69] 78 ± 6 ± 9 [84] 80 t^l ± 7 [96] 
55.4 ± 5.2 ± 4.0 [69] 68 ± 10 1^ [84] 89 ± 9 [96] 


3.1 ±0.7 

< 4.6 

77.6 ±4.6 

60.6 ±7.0 
{S = 1.3) 



Mode 


BaBar Belle CLEO 


Average 


Penguin- dominated decays 


B- w-K^ 
B- n^K- 


-17 ± 10 ± 2 [77] [84] 18 ± 24 ± 2 [95] 

-9±9±1[73] 23 ±11^4 [84] -29 ±23 ±2 [95] 

-10 ± 5 ± 2 [76] -7 ± 6 ± 1 [84] -4 ± 16 ± 2 [95] 
3 ± 36 ± 9 [75] — — 


-2 ±9 
(5 = 1.4) 

1± 12 
(5=1.8) 
-9 ±4 

3 ±37 


B- r]K- 
B- r]'K- 

Sn'Ks '■ 


-32+11 ± 1 [71] — 
3.7 ± 4.5 ± 1.1 [69] -1±7±1[83] 3 ±12 ±2 [95] 
-10 ±22 ±3 [69] 26 ±22 ±4 [86] — 

2 ±34 ±3 [69] 71±37i^[86] — 


-32± 20 

2±4 

8± 18 
(5=1.1) 

33 ±34 
(5 = 1.4) 
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Table 26: CP-averaged branching ratios (top, in units of 10~^) and CP asym- 
metries (bottom, in %) for B — > PV decays with AS — 1. Upper hmits are at 
90% confidence level. 



Mode 


BaBar Belle CLEO 


Average 


Penguin- dominated decays 


B- 'K-K*'^ 

B- TT^K*- 

B" TT+K*- 


15.5 ± [74] 19.4 t^Jl^i til [89] ± 1.6 (< 16) [94] 

— — < 31 [94] 
14.8tl:lil|to.9 [79] 16t^±2[92] 

— — < 3.6 [94] 


13.0 ±3.0 

( Q 1 A\ 

W - ■'-4) 

< 31 

15.3 ±3.8 

< 3.6 


B- -> K^p- 
B- K-p° 
B° ^ K-p+ 


— — < 48 [65] 
< 6.2 [74] < 12 [89] < 17 [94] 

7.2,t\i ± 1.3 [66] l?>-ltlit\it^A [79] 16.0 ± 2.8 (< 32) [94] 

— < 12 [88] < 39 [65] 


< 48 

< 6.2 

8.9 ±2.2 
{S = 1.4) 

< 12 


B- r}K*- 

B" rjK*" 
B~ -q'K*- 
B° r]'K*° 


22.1^1^:^ ± 3.3 [83] 26.5 1^:^ ± 3.0 [83] 26.4 1|:^ ± 3.3 [96] 
19.8 ± 1.7 [83] 16.5^4:^ ± 1.2 [83] 13.8 t^j ± 1.6 [96] 
— < 90 [83] < 35 [96] 
< 13 [83] < 20 [83] < 24 [96] 


25.4 ±5.3 
16.4 ±3.0 

< 35 

< 13 


B- K-uj 
B° ^ K'^LJ 
B- K-(l) 

50 ^ 


5.0 ± 1.0 ±0.4 [70] 6.7tJ-j±0.6 [80] 3.212-4 ± 0.8 (< 7.9) [94] 
5.3lJ:2 ± 0.5 [70] ± 0.5 (< 7.6) [80] 10.0 t^ls ± 1.4 (< 21) [94] 
10.0 1^:^ ± 0.5 [72] 9.4 ± 1.1 ± 0.7 [82] 5.5 j ± 0.6 [93] 

7.&l\i ± 0.5 [72] 9.0ti;8 ± 0-7 [82] 5.41^;^ ± 0.7 (< 12.3) [93] 


5.3 ±0.8 

5.1 ±1.1 

9.2 ±1.0 
{S = 1.4) 

7.7 ± 1.1 



Mode 


BaBar 


Belle 


CLEO 


Average 






Penguin 


■dominated decays 






50 








2&tf,t}l [90] 


26 ±35 


BO 




28± 17±8 [66] 


22^22+6 [79] 




26± 15 


B- 


K-uj 


-5 ± 16 ± 1 [70] 


6118 ± 1 [80] 




0± 12 


B- 


K-(j) 


3.9 ± 8.6 ± 1.1 [72] 


1 ± 12 ± 5 [83] 




3±7 


50 


S<l,Ks'- 


80 ± 38 ± 12 [67] 
-18 ±51 ±7 [67] 


-56 ± 41 ± 16 [86] 
-73 ±64 ±22 [86] 




19 ±68 
{S = 2.3) 

-38 ± 41 
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Table 27: CP-averaged branching ratios (top, in units of 10~^) and CP asym- 
metries (bottom, in %) for B — > PP decays with AD — 1. Upper hmits are 
at 90% confidence level. 



Mode 



BaBar Belle 

Tree-dominated decays 



CLEO 



B- tt-ttO 

^ TT+TT- 

B° ^ 7r°7rO 


5.5 tj;^ ± 0.6 [73] 5.3 ± 1.3 ± 0.5 [84] 4.6tJ:^lg:^ [91] 
4.7 ± 0.6 ± 0.2 [76] 4.4 ± 0.6 ± 0.3 [84] 4.5+J-^+g-| [91] 

< 3.6 [73] < 4.4 [84] < 4.4 [91] 


5.3 ±0.8 

4.6 ±0.4 

1.6 ±0.7 
< 3.6 


B~ — > 7r~ry 
B° ^ 7T°rj 

B- TT-f]' 

B° ^ n'^rj' 


4.2 IJ:^ ± 0.3 [71] 5.2 1?:? ± 0.6 [84] 1.2ti-J (< 5.7) [96] 

— — < 2.9 [96] 
< 12 [78] < 7 [83] < 12 [96] 

— — < 5.7 [96] 


3.9 ±0.9 
(5=1.1) 

< 2.9 
< 7 

< 5.7 



Penguin-dominated decays 



B- R-K^ 
B" ^ A'" A'" 



< 2.2 [65] 

< 1.6 [65] 



< 3.4 [84] 

< 3.2 [84] 



< 3.3 [91] 

< 3.3 [91] 



Pure annihilation decays 
B'^^K-K+ I < 0.6 [76] < 0.7 [84] < 0.8 [91] 



Mode 


BaBar Belle 


Average 


Tree-dominated decays 


B- tt-ttO 

B° TT+TT- 


-3+]?±2[73] -14±24l^[84] 
30 ± 25 ± 4 [76] 77 ± 27 ± 8 [85] 

2 ± 34 ± [76] - 123 ± 41 +5 


-7± 14 

51± 23 
(5 = 1.2) 

-49 ± 61 
(S = 2.3) 


B~ — > ■jT~ri 


-5lt?°±l[71] - 


-51± 19 
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Table 28: CP-averaged branching ratios (top, in units of 10~®) and CP asym- 
metries (bottom, in %) for B — > PV decays with AD = 1. Upper hmits are 
at 90% confidence level. 



Mode 



BaBar Belle 

Tree-dominated decays 



CLEO 



Average 



10.4+^;|±2.1 [94] 



B- 

B- 
B° 

BO 



* IT p 



9.3 ± 1.0 ± 0.8 [64] 
11.0 ± 1.9 ± 1.9 [64] 
13.9 ±2.7 [61] 
8.9 ±2.5 [61] 
22.6 ± 1.8 ±2.2 [66] 
< 2.5 [64] 



8.0t|^±0.7 [87] 



29.1 1^:^± 4.0 [81] 
6.0t||±1.2 [81] 



< 



43 [94] 



27.6 i?:4± 4.2 [94] 
< 5.5 [94] 



9.1 ±1.1 
11.0±2.7 
13.9 ±2.7 
8.9 ±2.5 
24.0 ±2.5 
< 2.5 



B — > TT a; 



5.4 ± 1.0 ±0.5 [70] 

< 3 [78] 



5.7l:\i±0.6 [80] 

< 1.9 [80] 



11.3 i^:^ ± 1.4 [94] 



"941 



5.9 ± 1.0 
{S = 1.2) 

< 1.9 



B- 
B° 
B- 
B^ 



> rip 

Tjp'^ 



< 6.8 [83] 



< G.2 [83] 

< 5.5 [83] 

< 14 [83] 



< 15 [96] 

< 10 [96] 

< 33 [96] 

< 12 [96] 



< 6.2 

< 5.5 

< 33 

< 12 



Penguin- dominated decays 



B- 
~B- 



^ TT (p 

>■ K-K 



*o 



< 0.41 [72] 



< 5 [97] 

< 5 [97] 
< 5.3 [94] 



< 0.4 

< 5 

< 5.3 



Mode 


BaBar 


Belle 


CLEO 


Average 


Tree-dominated decays 


B- 




-17 ±11 ±2 [64] 






-17± 11 


B- 


^ttV 


23 ± 16 ± 6 [64] 






23± 17 




TT'^ p~ 


-11 +^^±4 [68] 






-11± 17 


5° 




-62^24 ±6 [68] 






-62 ±27 


50 














Acp 


-18±8±3 [66] 


-38t^?l^ [81] 




-21 ±8 






36 ± 18 ± 4 [66] 






36± 18 






19 ± 24 ± 3 [66] 






19± 24 






28 ±4 [66] 






28± 19 






15 ± 25 ± 3 [66] 






15 ±25 


B- 




4± 17± 1 [70] 


48t^3±2 [80] 


-34 ± 25 ± 2 [95] 


9 ±21 
(5=1.8) 
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